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Preface 


At the request of the Defense Atomic Support 
Agency, Oak Ridge National Laboratory has under- 
taken the preparation of a handbook to aid engineers 
charged with the responsibility of designing shields to 
protect military equipment and personnel in the vicin- 
ity of a nuclear weapons burst. This document con- 
stitutes the fifth chapter of the Handbook issued thus 
far, the earlier chapters being Chapter 3, “Methods for 
Calculating Neutron and Gamma-Ray Attenuation,” 
Chapter 4, “Neutron and Gamma-Ray Albedos,” 
Chapter 5, “Methods for Calculating Effects of Ducts, 
Access Ways, and Holes in Shields,” and Chapter 6, 
“Methods for Predicting Radiation Fields Produced by 
Nuclear Weapons.” 

At the time Chapters 3, 4, and 5 were published it 
was planned that the Handbook would eventually be 
bound in two volumes, with Volume I serving as a 
textbook and ready reference and Volume II serving as 
an engineering guide. Volume I would include an 
introductory chapter, a chapter describing the radiation 
fields produced by nuclear detonations, and Chapters 3 
through 5. Volume II would present two or more 
chapters describing engineering methods based on the 
more sophisticated techniques described in Volume I. 
During the preparation of Chapter 6, however, the plan 
to combine the individual chapters into volumes was 
abandoned — at least for the present — the primary 
reason being that it would complicate the later publica- 
tion of updated chapters. We also changed the order of 
the chapters, with the result that the chapter originally 
planned as Chapter 2 was published as Chapter 6. One 
reason for shifting the chapter numbers was to improve 
the continuity of the Handbook by putting Chapter 6, 
which describes the various sources of radiation 
produced by a nuclear weapons burst and presents 
techniques for calculating the transport of the radia- 
tions from the point of burst to the surface of the 
shield, immediately preceding the engineering chapters. 
Another reason was that an additional chapter on basic 
shielding concepts, the present one, had been prepared 
and it would logically follow the introductory first 
chapter. With this reordering the basic information 
generally applicable to all neutron and gamma-ray 
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shielding problems is presented in Chapters 2 through 5, 
and the information specifically applicable to nuclear 
weapons radiation shielding problems is presented in 
the chapters beginning with Chapter 6. The engineering 
methods dealing with the design of shields for protec- 
tion against initial weapons radiation will be divided 
into two chapters: Chapter 7, entitled “Engineering 
Method for Designing Initial Radiation Shields for 
Blast-Hardened Underground Structures,” which will be 
published very shortly, and Chapter 8, “Engineering 
Method for Designing Initial Radiation Shields for 
Above-Ground Structures.” (Like Chapter 6, both 
Chapters 7 and 8 will be classified.) It now appears that 
methods for calculating dose rates in structures may be 
included as an additional chapter in the Handbook. 
Whether or not the Handbook will include a chapter 
dealing with shields to protect against fallout radiation 
is at present uncertain. If a fallout chapter is published, 
it will rely heavily on a method already available and 
will be integrated in the Handbook only so that all 
aspects of the radiation shielding problems will have 
been covered. 

In order to prepare this Handbook, it has been 
necessary for Oak Ridge National Laboratory to obtain 
the assistance of several consultants and subcontractors. 
For this chapter, for example, Paul N. Stevens, a 
consultant from the University of Tennessee, together 
with H. Clyde Claiborne of the Laboratory, prepared 
the first draft with which the editors worked. Other 
chapters similarly represent a cooperative effort of 
ORNL staff members and those of other organizations. 

As is always the case for handbooks, the authors and 
editors must rely heavily on the reviews of others as an 
aid in the development of the various chapters. The list 
of individuals who have contributed in this manner has 
already grown very large, and it would be almost 
impossible to acknowledge each person here. However, 
there are always reviewers who we feel have made such 
significant contributions as to warrant individual 
acknowledgement. With respect to this chapter, we 
particularly wish to acknowledge the help given by D. 
K. Trubey of the Laboratory, who spent considerable 
time reviewing and commenting on the entire chapter. 


Appreciation is also expressed to Maj. F. A. Verser 
and to Lt. Cols. Charles D. Daniel and William A. 
Alfonte, who as past DASA Shielding Project Officers 
handled the early administration of the contract and 
assisted in establishing the scope of the Handbook. 
The work they began is currently being ably performed 
by Maj. R. W. Enz. 
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Finally, we wish to thank Mrs. Virginia M. Hamrick, 
who assisted in the editing of the chapter, and Mrs. 
Virginia Glidewell, who typed and helped proofread the 
many drafts which are always necessary precursors to 
such a publication. 
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2.0. Introduction 


Weapons radiation shielding analysis usually involves 
three major steps: describing the field of radiation 
against which protection is required; determining the 
extent to which shielding materials will attenuate 
particles from that field; and converting the number of 
particles that penetrate the shield into a dose that can 
be related to a physical effect or to a biological 
response. This chapter describes the basic concepts 
underlying these steps. For those who plan to perform 
detailed shielding calculations by the methods described 
in Chapters 3 through 6, a thorough knowledge of these 
concepts is essential. Such knowledge is of less im- 
portance to those who plan to design shields merely by 
applying the engineering data provided in Chapter 7 and 
subsequent chapters, but it would be extremely helpful 
in all cases. 

The first step in the analysis — describing the 
radiation field — begins, of course, with an examination 
of the data available on the design and burst conditions 
of the weapon to be considered.* The resulting descrip- 
tion can vary from a simple statement of the total 
number of particles of a particular type at the “point” 
of detonation to a detailed specification of the energies 
of the particles and their locations and directions at all 
times over some specified volume. It is obvious that in a 
practical situation the consideration of volume 
introduces some type of shielding material (for ex- 
ample, the elements in the weapon or the atmosphere), 
and the effect that these materials have on the radiation 
field must be determined. As a result, the first two steps 
of the analysis become inseparable, the two together 
comprising a radiation transport problem. The custom 
is to divide the transport problem into two parts, one in 
which the attenuation by the atmosphere is calculated 
and another in which the attenuation by some barrier (a 
building or a specially designed shield) is determined. 

The various methods that have been developed for 
attacking the transport problem are described in Chap- 
ter 3 of this Handbook, and the application of some of 


*The interrelation of these parameters is discussed in detail 
in Chapter 6. 


the methods is demonstrated in subsequent chapters. 
The present chapter is devoted primarily to describing 
the quantities required as input to these methods and to 
discussing the types and behavior of the particular 
radiations to which the methods are applied in weapons 
radiation shielding. The chapter also discusses how the 
results of the transport calculations can be related to a 
physical or a biological effect, that is, how the third 
step of the analysis can be accomplished. 

Discussions of the quantities used to describe radia- 
tion fields are presented in Sections 2.1 and 2.2. The 
most detailed descriptions are given in  seven- 
dimensional phase space and are required in very 
sophisticated transport calculations. However, many 
transport methods cannot handle such detail and 
because of this, plus the fact that the detailed descrip- 
tions are so difficult to obtain, descriptions based on 
fewer parameters have come into common use. Also, 
some special descriptions of radiation fields that are 
basic to many of the mathematical models have 
emerged and these forms will often suffice. 

The extent to which shielding materials will attenuate 
particles from a radiation field depends both on the 
type and energy of the radiation and on the kind of 
materials comprising the shield. Some particles will pass 
through the shield unimpeded, and others will collide 
with the nuclei of the medium. Particle-nucleus colli- 
sions can result in a change in the direction and energy 
of the incident particle or in its absorption, the 
attenuation within the shield being determined by the 
probabilities for the various types of interactions. The 
quantities representing the probabilities, called cross 
sections, are discussed in Section 2.3. 

The formal definitions and explanations given in 
Sections 2.1-2.3 are presented in terms of idealized 
particles and can represent any type of neutral radia- 
tion. Insofar as shielding against nuclear weapons is 
concerned, the only radiations that need be considered 
are neutrons and gamma rays. While both these radia- 
tions can be regarded as neutral particles, there are 
certain obvious differences in their behavior. For 
example, the speed of a gamma-ray “particle” is 
constant and does not depend on its energy, whereas 


the speed of a neutron is determined by its energy. 
Thus some caution must be exercised when specializing 
particle concepts to gamma rays. 

Neutrons and gamma rays can undergo many different 
interactions with materials, and calculations of their 
transport through a shield can be quite difficult. The 
situation is complicated by the fact that the interaction 
of a neutron can produce a gamma ray. The converse is 
also possible, but is much less likely. The types of 
interactions that neutrons and gamma rays can undergo 
are discussed in Sections 2.5 and 2.7 respectively, and 
the interactions that give birth to these radiations are 
discussed in Sections 2.6 and 2.8 respectively. 

Not all quantities that can be determined by transport 
calculations (for example, differential flux) are easily 
measured or easily related to a physical or biological 
effect. However, a quantity that has come to be known 
as dose and that is determined from the energy 
deposited by the radiation in the receiver (for example, 
in an instrument or in a human body) is measurable, 
and calculated results can be related to the same 
quantity. As a result, whenever a physical or a 


biological effect is being considered, dose is the 
quantity that is used to specify the design limit on the 
intensity of the radiation field. 

Unfortunately, several concepts of dose exist, and 
considerable confusion has arisen with respect to their 
definitions. Because the problem has not yet been fully 
resolved, all viewpoints are presented in Section 2.9 of 
this chapter and the definitions that are employed in 
this Handbook are specified. Section 2.10 discusses 
allowable dose limits and the probable effects of such 
acute doses as might be encountered in a nuclear 
accident or in warfare. 

Finally, it is pointed out that although this chapter 
was prepared primarily to describe basic concepts and 
particle interactions of especial interest to the designer 
of weapons radiation shields, its contents will be 
equally useful to those engaged in reactor shield design 
since the radiations involved are the same. The informa- 
tion will also be useful to those who design radiation 
shields for accelerators or spacecraft, but for these 
purposes it is incomplete since many more different 
types of radiations must be considered. 
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2.1. Quantities Used to Describe Radiation Fields 


A complete description of a radiation field would 
specify the number of particles of a given type and 
energy which at a given time exist at a given position 
and travel in a given direction, with all particle types, 
energies, positions, and directions being considered at 
all times. However, such complete descriptions are 
seldom required — except in the most sophisticated 
particle transport calculations — and a radiation field 
cannot be measured in such detail. Therefore particle 
populations are usually described in terms of a fewer 
number of parameters. The quantities used most 
frequently in shielding analysis are listed in Table 2.1 
and discussed in the following paragraphs. The names 
given to the differential quantities have not been 
standardized, and occasionally one finds that a term 
specifically identified as being differential is actually 
less differential than a companion term not so 
designated. While in the first column of Table 2.1 an 
attempt is made to specify the degree of differentiation 
in each case, the names used in the text are those in 
keeping with general usage. 


2.1.1. PARTICLE DENSITIES 


A knowledge of the particle density over all phase 
space is equivalent to a complete solution of a particle 
transport problem and comprises more information 
than is available from most calculational schemes now 
in use. When given in seven-dimensional phase space, 
which consists of three spatial coordinates, the par- 
ticle’s energy, two direction-defining angles, and tie, 
this particle density is defined by 


n(r,£,0,t) dE dQ=the number of particles per unit 
volume at space point 7 and time ¢ 
having energies in dF about energy 
E and directions in dQ about the 
unit direction vector 2. 


The particle density so defined is doubly differential — 
in energy and in direction — and less detailed forms will 
often suffice. For example, one may use the steady- 


state* particle density differential in energy only, 
commonly referred to as differential particle density 
and defined by 


n(r,E) dE =the number of particles per unit volume at 
space point 7 having energies in dE about F 


= f. n(7,E,) dQ dE . (2.1) 


Or one may use the steady-state total particle density, 
defined as the number of particles per unit volume at 
space point F and given by 


n(r) = f ; f ” n(7,E,9) dE dQ = f, “n(7,E) dE. (2.2) 


2.1.2. FLUX DENSITIES 


Even though the concept of particle density is 
basically simple and has a unique interpretation, experi- 
ence has shown that the flux density’ or, as it is more 
commonly referred to, the flux,’ serves better as the 
dependent variable in solutions of the transport equa- 
tion (refer to Chapter 3, Section 3.1). The flux is 
related to the particle density through the particle’s 
speed, and when described in terms of seven- 
dimensional phase space the flux is given by 


€(7,E,0,1) =v n(,£,Q,0) , (2.3) 


where pv is the particle’s speed and corresponds to the 
energy E. (The speed is the scalar magnitude of the 
particle’s velocity vector v.) 


*The steady-state, or time-independent, condition will be 
denoted in phase space notation by dropping the time symbol f; 
for example, seven-dimensional phase space becomes six di- 
mensional and in phase space notation is given by (7,£,). 


tT Aithough this quantity is truly a density and the Inter- 
national Commission on Radiation Units and Measurements 
recommends the use of the term flux density, the simpler term 
flux is ingrained in shielding terminology and will be used in 
this and other chapters of this Handbook. 


Name 


Table 2.1. Quantities Used to Describe Particle Populations 


tit 1 i 
Quantity Cominonly Weed Symbol Usual Units Comments 
Particle density Particle density nr,E, Q,t) particles nv,E,Q,t) dE dQ =number of particles per unit volume at the space point 7 and time ¢ 
differential in cm?-MeV'steradian which have energies in dE about energy E and directions in dQ about the unit vector 2. 
energy and 
angle 
Particle density Differential n(r,E) particles n(,E) dE =f n(r,E,2) dQ dE. 
differential in particle cm3-MeV an 
energy density n(,E) dE = number of particles per unit volume at the space point r that have energies 
in dE about £. 
Total particle Particle density nr) particles/em3 n(r) = f n(,E) dE. 
density 0 
n(r) = number of particles per unit volume at the space point r. 
Particle flux Angular flux o£, 2,0 particles 0(7,E,2,t) =y n(r,E,9,t) = particle speed (vy) times the particle density; that is, the 
. : _ = SS ee . . . 3 z ie : 
feta! differ cm? -sec‘MeV'steradian ies nk the angler flux derive from (cm/sec) (particles/cm? "MeV-steradian). 
ae a — 0(7,E,Q,t) dE dQ = flow per unit area and time at space point r and time ¢ of particles 
eae having energies in dE about E and directions in dQ about 2 
or 
O(7,E,2,t) dE dQ = total tracklength traced by particles per unit volume and time at 
the space point r and time t which have energies in dF about E and directions in dQ 
about 9; that is, the angular flux is viewed as (particles:cm/sec)/cm?-MeV‘steradian. 
Particle flux Differential OE) particles (FE) dE = f 0(E,Q) dE dQ. 
: P SS 4n 
density differ- flux cm? -sec‘MeV - 
ential in @(r,E) dE = number of particles having energies in dE about E which enter a sphere of 
energy unit cross section per unit time at the space point F 
or 
o(r,£) dE = total tracklength traced by particles per unit volume and time at the space 
point r which have energies in dE about E. 
Total particle Flux or) particles o(r) = f @(r,£) dE. 
flux density 2 0 
em**sec 


(Fr) = number of particles which enter a sphere of unit cross section per unit time at 
the space point 7 
or 
(r) = total tracklength traced by particles per unit volume and time at the space 
pointr. 


Table 2.1. (continued) 


Quantity 


Group flux den- 
sity differential 
in angle 


Group flux 
density 


Discrete ordinates 
flux density 


Energy flux 
differential in 
energy and 
angle 


Energy flux 
differential in 


energy 


Total energy 
flux 


Fluence 


Name 


Commonly Used Symbol Usual Units Comments 
—— . pies Eg ie 
Group angular @o(7,2), particles Otr, 2) =f o(,£,Q) dE. 
flux 2. ne A Ege 
cm*-sec’steradian 
Oo(r,n) particles @o(7,2) = angular flux within the energy group AEF ¢. 
cm?-sec-unit cosine ®o(7,4) du= Oo(7,2) dQ = total tracklength per unit volume and time at the space 
point r of particles having energies within the energy group AE @ and directions 
defined by direction cosines in du about p, where p = 27r/|r}. 
= . ss Eg. a 
Group flux ® Gg) particles Oo) = f @(r,£) dE. 
2. Egtl 
cm**sec _ 
®,(7) = total flux within the energy group AF ¢. 
: . : _ oil Ug = 
Discrete ordi- ®;GD particles %;¢ p=—— Pou) du. 
nates flux 2. : : "Abn 4 Bar) 
cm*‘sec"unit cosine 
®; ¢ p =sroup angular flux averaged over specified segments of u-space, Aup, and 
corresponding to discrete spatial mesh points ri. 
Angular energy IGE,9,t) MeV 1,E,2,0dE dQ=E 0(,E,9,t) dE dQ= energy flow per unit area and time at the 
flux em?-sec-MeV'steradian space point r and time t due to particles which have energies in dE about E and 
directions in dQ about Q; here the units of EF ®(7,£,9,t), which is equivalent to 
1(r,E,9,t), are (MeV/particle) (particles/cm? -sec-MeVsteradian). 
Differential IGE) MeV I¢,E) dE = E &(7,E) dE = energy flow per unit time into a sphere of unit cross 
energy flux em2-sec‘MeV section at the space point 7 due to particles which have energies in dE about E. 
Total energy Ir) MeV IQ) = f I(7,E) dE: 
= a. 0 
flux cm? sec as : ; : ; 
I(r) = total energy flow per unit time into a sphere of unit cross section at the 
space point Fr. 
: é ty+At 
Fluence wy Appropriate units for = f (ft) dt. 
ty 


flux integrated over 
a specified time 
interval At 


y = time-integrated flux, where ®(f) is any of the fluxes listed above. 


e_ 


Table 2.1. (continued) 


Name - 
i 1 Usual Unit 
Quantity Comintonly Used Symbo sual Units Comments 
Current differential Angular current J(r,v,2,t) , particles J(,v,2,0) dv dQ=Qy n(r,v,2,0) dv dQ = directed flow per unit area (normal to the 2 
in energy and Tem? "sec" ni /sec)-steradian direction) and time at the space point r and time ¢ of particles having speeds in dy 
angle about v and directions dQ about Q. 
TE, 2,0) particles JVE,2,0 dE dQ = 2 v n(E,9,1) dE dQ = 20(,E,,1) dE dQ = direction vector N 
cm?2-sec'MeV'steradian times the angular flux. 
Current differential Total angular 70,9) particles T(r,Q) = f J(,E,Q) dE=20(7,Q). 
in angle current cm?-sec’steradian fo a age = 
J(r,2)d2 = directed flow per unit area (normal to the Q direction) and time at the 
space point r and time ¢ of particles having directions in dQ about Q. 
Net current differ- Differential net J(,E,0) particles J(,E,t) dE = f J(V,E,2,t) dQ = f 2OG¢E£,,t) dQ. 
ential in energy current cmt sen Me 4n 4n 
J(E,t) dE = net flow per unit area and time at the space point r and time ¢ of particles 
having energies in dE about E where the unit area is normal to the direction of the 
resultant vector J(r,£,t) (which is simply the angular flux-weighted vector summation 
of the unit vectors 2 over 4m solid angle). 
Total net current Net current J) particles JM) =f T7,2) da -f 2,2) aN. 
———— 4n 4n 
cm? “sec Fr. . . . . a . 
J(r) = net flow of particles per unit area and time at the space point r where the unit 
area is normal to the direction of the resultant vector. 
Cot a hi E, aes = = 
Group current Group angular I G(r, 2) particles Jolr, 2) -f 8 J(r,E,Q) dE = NOG(,2). 
ecu in current oni see-steradian os oe Ett = 
nee Jglr, 2) dQ = the directed flow per unit area (normal to the Q direction) and time at 
the space point and time ¢ of particles having energies within the energy group 
AE g and having directions in dQ about &. 
Group net current Group net IG& particles IG” = f Ig ® dQ= f 2OG,Q) dQ. 
density current int ssee an an 


ie) = net flow of particles per unit area and time at the space point r having energies 
within the energy group AE, where the unit area is normal to the direction of the 
resultant vector. 


a 
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The flux defined in the preceding paragraph is doubly 
differential and is usually referred to as the angular 
flux. Greater insight into the use of the angular flux as 
the dependent variable in mathematical descriptions of 
particle transport is provided by its interpretation either 
as the tracklengths traversed per unit volume and time 
or as the flow of particles per unit area and time. 

The tracklength interpretation of angular flux follows 
from the observation that the speed of an individual 
particle can be considered as its scalar tracklength per 
unit time. The product of particle density and speed is 
then the sum of the tracklengths traced by all the 
particles within a unit volume per unit time,* in which 
case the definition of the angular flux would be 


€(7,£,Q,t) dE dQ = the total tracklengths traversed per 
unit volume and time at space 
point 7 and time ¢ by particles 
having energies in dE about energy 
E and directions in dQ about 22. 


The interpretation of the angular flux as a flow of 
particles per unit area and time is closely related to the 
concept of angular current (see discussion of angular 
current in Section 2.1.4). It will be shown in the 
discussion on current that the angular flux is identical 
with the magnitude of the current vector J and thus can 
be interpreted as 


(7,E,9,1) dE dQ = the flow per unit area and time at 
space point 7 and time t¢ of particles 
having energies in dE about EF and 
directions in dQ about 2. 


When the transport and deposition of the particle’s 
kinetic energy are of concern, the energy flux differ- 
ential in energy and angle is often used. This quantity, 
referred to as the angular energy flux, is defined by 


1(7,£,0,0) dE dQ = the energy flow per unit area and 
time at space point 7 and time ¢ due 
to particles having energies in dE 
about E and directions in dQ about 
Q, 


*Note that the per-unit-time units of the flux are associated 
with the particle’s speed, which is a function only of the energy. 
However, the time dependence of the flux is a consequence of 
the time behavior of the particle density, which does not have 
time units even though the time symbol ¢ is included in the 
phase space notation to denote a dependence on time; that is, 
n(7,E,8,t). 


and is related to the angular flux by 
I(7,E,9,0) = E ®(,E,Q,t) . (2.4) 


Many calculational models employ less detailed de- 
scriptions of the flux; for example, they may use a 
description in which steady state is assumed and the 
angular dependencies are removed by the appropriate 
solid-angle integration. This reduces the description to 
four-dimensionat phase space and results in a dependent 
variable representing the particle flux differential in 
energy. Commonly known as the differential flux, this 
quantity is given by 


®(7,£) = S, _ OED) dQ. (2.5) 


Like the angular flux, the differential flux can be 
interpreted in terms of tracklength per unit volume and 
time or in terms of the number of particles that enter a 
unit sphere per unit time. In the latter case, the 
solid-angle integration can be regarded as a summing of 
particles that enter a sphere of unit cross section, 
regardless of their directions of motion. The sphere is in 
effect generated by the rotation of a circular unit area 
during the integration over 47 solid angle (see Fig. 2.1). 
In this context the definition for the differential flux 
can be restated as 


&(r,E) dE = number of particles having energies in dE 
about £ which enter a sphere of unit cross 
section per unit time at space point 7. 


While this definition of the differential flux is descrip- 
tive, it is not exact. The mathematically rigorous 
definition is 


lim N(E) dE 


@(7,E) dE = 
(rE) Aas a 8 


(2.6) 


which implies the limit process AA > 0, with M(E) dE 
denoting the number of particles having energies in dE 
about £ which enter an “incremental” sphere of cross 
section AA per unit time. 

- The concept of the incremental sphere is the best way 
to visualize the energy flux differential in energy. 
Referred to as the differential energy flux, this quantity 
may be defined by 


I(r,E) dE = the energy flow per unit time into a sphere 
of unit cross section at space point 7 due to 
particles having energies in dE about E. 
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Fig. 2.1. Incremental Sphere Concept of Flux. 
The differential energy flux is given by 
I(7,E) dE = f. 1(V,E,9) di dE=E ®(,E)dE. (2.7) 
Tv 


Other quantities used are the fotal flux, the total 
energy flux, the group flux, the group angular flux, and 
the discrete ordinates flux. The total flux, defined al- 
ternatively as the total particle tracklength per unit 
volume and time at space point 7 or as the number of 
particles that enter a sphere of unit cross section per 
unit time at space point 7, is obtained by integrating 
the differential flux over all energies: 


(7) = J, ” (RE) dE. (2.8) 


Similarly, the total energy flux, defined as the total 
energy flow per unit time into a sphere of unit cross 
section at space point 7, is obtained by an integration of 
the differential energy flux over all energies: 


I(F) = S "1(F.E) dE . (2.9) 


The total flux has only limited application to practi- 
cal shielding problems because of the strongly energy- 
dependent nature of the particle’s behavior. A more 
useful approach is to divide the total energy range into 
NOG energy intervals, called energy groups, 


AE Gg =E,- Eguy* G=1,2...NOG, 


*The subscripts g and gt+1 refer to the upper and lower limits 
respectively of the Gth energy group, and G=1 corresponds to 
the highest energy group. An alternate convention would as- 
sociate G=1 with the lowest energy group and the subscripts 
gt+l and g would then correspond to the upper and lower 
energy limits. 


and to define the group flux as the integral of the 
differential flux over the corresponding energy group: 


E 
(N= 1% BFE\dE, (2.10) 
e i 1 
with the obvious constraint that 
=P Oe). (2.11) 


NOG 


The group angular flux (group flux differential in 
angle) has a similar definition and is obtained by 
integrating the angular flux over a specific energy 


group: 


og67,0)= f°% oe Dae. (2.12) 
Eg+i 


For problems that involve directional symmetry, the 
group angular flux can be rewritten in terms of a new 


angular variable uw = Q2-r/|r|, the direction cosine: 


Bo (ry) du = ¥ (7,0) dQ. (2.13) 


The group angular flux can then be defined as the total 
particle tracklength per unit volume and time at space 
point 7 of particles with energies within energy group 
AE and directions defined by direction cosines which 
lie in du about yu. 

A mathematical treatment which evaluates the angu- 
lar flux at discrete angles, called discrete ordinates 
forms the bases for many detailed numerical solutions 
of the transport equation. A form of discrete ordinates 
could be simply group-angular-flux values that corres- 
pond to given spatial positions and particle directions. 
However, a more useful form is obtained by integrating 
the group angular flux over specified segments of py 
space, Aun = Mg ~Mgy,,D=1,2...NOA, 


= af el a Z 
® = 7 _ Pou) du. 2.14 
60> Ra fi" ScGudu. 2.14) 


The “mean value” approximation to Eq. 2.14 is given 
by 


&g pM =FoFutp), (2.15) 


© 


© 


ge 


where Lp is a mean value* for the direction cosine over 
the increment Aypy. The discrete ordinate fluxes 
®; ¢,p! are those values of @g p(F) which correspond 
to the spatial meshpoints 7;: 


®; 6,0 =%,v%) - (2.16) 


2.1.3. FLUENCES 


A useful measure of total exposure to a flux for 
applications involving energy deposition or dose is the 
integral quantity called fluence. Fluence is defined by 
the International Commission on Radiation Units and 
Measurements as the quotient of AN divided by AA, 
where AN is the number of particles that enter a sphere 
of cross-sectional area AA and the A’s imply a limiting 
process.’ This definition is equivalent to regarding 
fluence as a time-integrated flux over some specified 
time interval. As such, the fluence can be written as 


v =f" wate, (2.17) 


where At corresponds to some specified time interval 
and ®(t) can be any one of the several kinds of flux 
described in Section 2.1.2. For example, the energy 
fluence is 


F=f TED dt. (2.18) 


(NOTE: The symbolism for particle fluences and 
fluxes has been confusing in the past, with the symbols 
nvt and ® being used for both. In order to avoid 
confusion in this chapter, ® is reserved for particle 
fluxes and when particle fluences are intended, the 
symbol yw is used.) 


*up may be adjusted to give the equality; for well-behaved 
functions the closer up is to the true mean, the better the 
approximation. 

tThe following subscript notation is used throughout this 
section and in Chapter 3: subscripts J, G, and D denote 
functions whose values are associated with the /th space inter- 
val, the Gth energy group, and the Dth angular interval, respec- 
tively; i and i+1 refer to a function evaluated at the lower and 
upper limits of the Jth space interval; g+1 and g refer to a 
function evaluated at the lower and upper limits of the Gth 
energy group; d+1 and d refer to a function evaluated at the 
lower and upper limits of the Dth angular interval. 


2.1.4. CURRENT DENSITIES 


The characteristic property of the current variable is 
its close relation to the convective (leakage) effects in 
the theoretical description of particle transport. The 
most general form of the current variable is differential 
both in energy (or speed) and in angle. Called the 
angular current density or, more frequently the angular 
current, this quantity is symbolized by J(r,v,Q,1) and is 
defined as the directed flow per unit area (normal to 
the Q direction) and time at the space point 7 and time 
t of particles having speeds in dy about v and directions 
in dQ. about 9. 

The relation between current and particle density can 
be established by considering that (1) the product of » 
and the particle density can be regarded as a vector sum 
of the individual co-directional velocity vectors (¥), 
yielding the resultant vector J =n v, or (2) the (v X dt 
X dA X n) particles contained within the volume 
element shown in Fig. 2.2 will all exit through the 
differential area dA within differential time dt if v = 
vQ. These models can be expressed mathematically in 
terms of the angular current as 


Jr, v,2,0) dv dQ. = 2 v n(7,v,9,0) dv dQ , (2.19) 


and since the particle’s kinetic energy is a function of 
its speed, Eq. 2.19 can be rewritten as 


WE,2,t) dE dh = Dv n(7,£,8,0) dE dB , (2.20) 


where J(7,£,0,t) dE = J(r,v,Q,t) dv and n(7,E,Q,t) dE = 
n(r,v,Q,t) dv. Then when v n(7,£,Q,f) is identified as 
the angular flux, the relation between the angular flux 
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Fig. 2.2. Particle Flow Concept of Current. The particles 
contained within the volume element will all exit through the 
differential area dA within dt if vy = vQ, 


and the angular current pointed out in Section 2.1.2 is 
obtained: 
IC,.E,9,1) dE dQ=2.6(,£,9,) dEdQ. (2.21) 


Other less detailed forms of the angular current are 


the group angular current Ig(,Q) and the total angular , 


current J(r,Q), which are obtained by integrations over 
an energy group and all energies, respectively: 


I(,E,Q) dE = 2. O6(7,Q) 5 (2.22) 


ant “see: E 
To@.B)= f g 
gt 


j¢,0)= Sf "FGED) dE =0. 47,2). (2.23) 


It is apparent that the angular current variables all have 
essentially the same simple relation with the corres- 
ponding angular flux variables, because the energy 
integrations are performed directly on the flux vari- 
ables. For example, in the case of the group angular 
current, 


Jo(7,0) = i 2 0(7,£,2) dE 
Eg+i 


=2f "8 06,6,2) dE=2466,H). (2.24) 
Eos) 

The integral of the angular current over all directions 
(47 solid angle) constitutes a vector summation, and the 
resultant vector is regarded as the net current, often 
called simply current. 


The net current differential in energy only, referred 
to as the differential net current, is defined by 


J(r,£,t) dE =the net flow per unit area and time at 
space point r and time ¢ of particles having 
energies in dE about E, where the unit 
area is normal to the direction of the 
resultant vector J(r,E,t) 


. f IE, O,t) dQ dE 
4n 


= f. 2 O(,6,0,1) dQ dE , (2.25) 
WT 


which is the angular flux-weighted vector summation of . 


the unit vectors {2 over 47 solid angle. The group net 
current J(r) and total net current J(r) ate 47 solid- 
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angle integrations of the group angular current and total 
angular current, respectively: 


Ig) = J IeG@D d= f 2&_69)dQ; (2.26) 


I(r) f, IGM) aQ= f, 2 OGD) dO. (2.27) 

The flow of Q-directed particles across an arbitrarily 
oriented differential area is a necessary concept in the 
description of the directed flow of particles in terms of 
a specific coordinate system and can be related to the 
angular current J(Q)* by consideration of Fig. 2.3, 
where the direction vector n is normal to the differ- 
ential area. The number of {-directed particles crossing 
the differential area dA per unit time is equal to 
J(Q)-(n dA). A scalar current J,,() that describes the 
flow of the {2-directed particles per unit area normal to 
the direction a is defined as 


J, (Q) dA = J(Q)(n dA) ; (2.28) 
it follows that 
J,(Q) = aT) = Bn &@) = cov 6 BG), (2.29) 


where 


n= the unit vector corresponding to an arbitrary 
direction, 


n=a coordinate-identifying subscript, for ex- 
ample, n =x when n =i, 


J,,(Q) = the flow of Q-directed particles per unit area 
(normal to the direction 7) and time, 


B(Q) = = the angular flux variable corresponding to the 
angular current J(Q). 


Note that J,,(Q) is a scalar quantity but is uniquely 
related to the coordinate system through the direction 
unit vector n. The corresponding vector current is given 


by 


J,(Q) =n J, (Q) =n Q)n , (2.30) 


where the vector J,() is the component of the vector 
J(Q) with respect to the n direction and J,,(2) is the 


*In this further discussion of current, the notation IQ) will 
be used to denote any of the angular currents and J the 
corresponding net current: for example, J(Q) may represent 
IGE, 2 vt), IgG 2,t), or J(7,2) and J may represent J(AE,b), 
Ig, orJ(n). 


 Y 


ge: 


Cc 
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Fig. 2.3. Schematic of Particle Flow Across an Arbitrarily 
Oriented Surface, 


dA (DIFFERENTIAL AREA 
NORMAL TO 7 ) 


projection of the vector J(Q) on the 7 axis and also the 
magnitude of the component J,,(®). 
The three components of J(Q) in cartesian coordi- 
nates are given by 
J,(Q) = 7 ,Q) = 7 OQ) FQ) , 
7,(Q) =F JQ) =F @) FD), (2.31) 
T,(2) =F ,(Q) =k &@) (ED), 


and J(Q) is equal to the vector sum of the three 
components: 


IQ) = Fy(Q) + FQ) + FQ) 
=i5,(Q) + jJ,Q) +k J,(Q). (2.32) 
The net flow of particles across a unit area normal to 


the 7 direction is obtained by integrating J,,(Q) over all 
directions: 


: f _InQ) dB = f _AIQ) d= nT. (2.33) 


Further, J, can be represented in terms of the partial 
currents* J and J, such that 


J,2dt Jn, (2.34) 


*The partial currents are usually defined to be positive 
numbers, with their difference Up -J n) determining the sign of 


J: 
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where 


J; =the flow per unit area (normal to the 7 direc- 
tion) and time of particles having positive 
direction cosines relative to the n direction 


=fa-ai>oin® dQ ? 


J, =the flow per unit area (normal to the 7 direc- 
tion) and time of particles having negative 
direction cosines relative to the n direction 

= fan <oin(DaM. 

The partial currents do not occur explicitly in most 
theoretical formulations of particle transport; however, 
the boundary conditions of many problems of interest 
are best expressed in terms of the partial current 
variable. 


2.1.5. ADJOINT FLUX 


The adjoint flux does not correspond in the usual 
sense to a real physical quantity. Instead it is the flux 
associated with the adjoint equation when that equa- 
tion is interpreted to’ describe the transport of hypo- 
thetical radiation particles, sometimes called “adjunc- 
tons.” The rationale of the concept of adjunctons 
follows at least in part from the similarity of the adjoint 
equation to the Boltzmann transport equation.t Inte- 
gral forms of these two equations are given by Eqs. 2.35 
and 2.36 respectively, and a coupling between the two 
is provided by Eq. 2.37: 


&*(P) = [K(P>P') OP) aP +f), (2.35) 

x(P) = a x(P’) K(P’ > P) dP’ + S(P) , (2.36) 

A= fs) o*@) aP = fx(P) f(P) aP , (2.37) 
where 


®*(P) =the adjoint flux, defined as the ex- 
pected value of the present and future 
contributions to the effect of interest 
by a radiation particle leaving a col- 
lision at a point whose coordinates in 
phase space are denoted by P, 


tT Derivations, relations, and other useful and interesting 
properties of the various forms of the Boltzmann equation are 
given in Chapter 3 of this Handbook. 
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K(P' > P) dP =the probability that a particle experi- - 


encing an event at a point P’ in phase 
space will experience another event at 
a point whose phase space coordinates 
lie within dP about P, 


f@)=the payoff function, which is the 
probable, or expected, direct contri- 
bution to the effect of interest by a 
particle leaving a collision at point P in 
phase space, 


x(P) dP = the expected number of particles leav- 
ing collisions at points whose phase 
space coordinates lie within dP about 
P, x(P) being called the “emergent 
particle density,” 


S(P)dP=the expected number of _ initial 
(source) events which occur within the 
phase space coordinates dP about P, 


A=the effect of interest which can be 
calculated if either (the adjoint flux 
and source function) or (the emergent 
particle density and payoff function) 
are known, 


The adjoint equation, which mathematically defines 
the characteristics of the adjunctons, has often been 
referred to as the “backward” equation because the 
movement of adjunctons is backward to that of real 
particles. That is, they originate at a point in phase 
space associated with the detector and move in the 
direction of phase space that corresponds to the real 
source. This behavior requires that in a collision the 
adjuncton gain energy, which contradicts the usual 
energy conservation principles unless prior to the 
collision the target nucleus has the kinetic energy it 


would have acquired in a real collision process and the 
particle has the final energy of the adjuncton. This 
precludes the assignment of physical properties to the 
adjunctons that would resemble those of real particles. 

The lack of physical properties is not a serious 
shortcoming, however, since adjunctons are not useful 
as radiation particles per se. Rather, the adjoint flux 
which may be mathematically evaluated for some point 
in phase space has the precise interpretation of “value 
function.” The value function has been defined to be a 
measure of the present and future contributions of a 
real particle at some point in phase space to an effect of 
interest and as such has been used to bias the Monte 
Carlo analysis of the behavior of real particles (neutrons 
and gamma rays). This is accomplished by modifying 
the basic stochastic processes so that event densities 
that contribute most strongly to the effect of interest 
are sampled most frequently. The transmission proba- 
bility used in the Monte Carlo analysis is modified by 
using the value function as an importance function, 
with the particle’s “weight” being adjusted to preserve 
parity with the true analog calculation. 

Calculation of the answer of interest through Monte 
Carlo solution of the adjoint problem avoids the 
difficulties associated with estimation to a point detec- 
tor, provided, of course, that the source has a more 
agreeable configuration. Also, the adjoint concept 
provides a fundamental basis for perturbation analyses. 

Since the adjoint equation is very similar to the usual 
Boltzmann transport equation, it is not surprising that 
methods developed to solve the latter will also apply to 
the former. For example, adjoint fluxes have been 
calculated using the techniques of diffusion theory, 
Monte Carlo, and transport theory. Consequently, there 
should be no serious practical or fundamental diffi- 
culties in calculating the adjoint flux. 


© 


c 
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2.2. Special Radiation Fields 


Certain radiation fields are basic to many of the 
mathematical models that describe particle transport, 
and they possess functional simplicities that sometimes 
permit analytical solutions. In this section four such 
fields are presented in terms of the various flux and 
current variables. They are a monodirectional beam, an 
isotropic flux field, a point source, and an isotropic 
plane source. 


2.2.1. MONODIRECTIONAL BEAM 


The monodirectional beam, shown schematically in 
Fig. 2.4, consists of a one-dimensional system of 
particles that are all moving in the same direction Qo. It 
is described in most detail in terms of the angular 
current: 


I(r,E,Q) dE dQ= QIy FE) 6(2,29) dEdQ, (2.38) 


where Jo is the magnitude of the current (particles 
cm? sec!) for a unit area normal to the Qp direction, 
F(E) is the normalized energy spectrum with the 


property 

ie F(E) dE =1, (2.39) 
and 8(Q,Q9) is an angular delta function* with the 
property 
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Fig. 2.4. Monodirectional Beam Incident on the Shield. 


J, (2) 88,2) dQ = 2) , (2.40) 


a mathematical operation that aligns all direction 
vectors with Qo. 

When the monodirectional beam is specified to be 
monoenergetic, which it frequently is for calculational 
convenience, the Dirac delta function with the property 


J RE) 8E - Eo) dE = flE) (2.41) 


may be used to describe the energy “spectrum” of the 
beam. That is, 


J(r,E,Q) dE dQ 
= QJ 6(E — Eo) (2,29) dEdQ. (2.42) 


Integration of Eq. 2.42 over all directions (47 solid 
angle) yields the differential net current 


I B)= J Hr #,D a 
=JoF(E) S 8(2,29) dQ 
= Qo Jo FE), (2.43) 


and a subsequent integration over all energies gives the 
total net current 


=f" Wee) ak 
= Qo Jo J ” FUE) dE = 29 Jo (2.44) 


*A delta function is simply a mathematical notation to indi- 
cate that the function is non-zero at only one particular 
value of the argument. 
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The angular flux* for the beam is equal to the scalar 
magnitude of the angular current and is given by 
(7,E,Q) =Jo 6(E ry Eo) (QQ) . (2.45) 


The differential flux is obtained by integrating the 
angular flux over 47 steradians, 


HGEE)=f WEED) ad 


 =Jy B(E — Eo), (2.46) 
and the total flux is obtained by integrating the 
differential flux over all energies, 


OP) = Jy OEE) dE = Jy. (2.47) 


The beam concept is primarily used to specify a 
source configuration at a surface boundary and as such 
is usually regarded as a boundary condition. The partial 
current with respect to the positive x direction is 

Tg = (Qoi) Io (2.48) 
The partial current in the negative direction, J,, cannot 
be used as a boundary condition in this situation 
because the reflected particles would also contribute to 
the partial current in the negative direction, which is 
unknown prior to the solution of a problem with an 
incident beam source configuration. 


2.2.2, ISOTROPIC FLUX FIELD 


The radiation field for many physical situations can 
be regarded as being isotropic with respect to the 
angular distribution of the particle velocity vectors; that 
is, motion in any direction is equally probable. The 
angular flux for this condition is simply the differential 
flux ®o(r,E) divided by 47 steradians, 


Po(r,£) 
7 ry 


(rE) = 
O7,E, 2) A 


(2.49) 


*The flux, usually referred to as the “incident flux,” would 
not be the flux existing in a real situation since backscattering 
would contribute to the flux at a surface. 
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and the total flux is the differential flux integrated over 
all energies, 


Bo(r) = f “ o(r,£) dE . (2.50) 
0 
The angular current is 
eaten ee VS ees Dar 
We. 6,D) = 0 0G, 6,9) =e) (2.51) 
T 


and the differential net current is obtained by inte- 
grating the angular current over all directions. For an 
isotropic flux field the differential net current is 
characteristically zero: 


Ki,E) = JG,B,Q) dQ 
®o(r,E) OHO 
= NdQ=0. 2.52 

4n fi. ( ) 

Under the concept of an isotropic flux field the 
positive partial current is sometimes used to specify an 
equivalent surface source configuration. Consider, for 
example, a unit area normal to the direction vector 7 as 
shown in Fig. 2.5. The scalar angular current that 
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Fig. 2.5. Orientation of Q-Directed Particles with Respect to 
a Unit Area Normal to the Direction Vector 7. 
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describes the flow of 9-directed particles in terms of a 
unit area normal to the 7 direction is given by 


I(t EQ) = IR E,Q) = 0-Q OED), — (2.53) 


and the positive partial current that corresponds to the 
n direction is obtained by appropriate integrations: 


TGO=S a. S. In@#D ae aD 


= bey a na PO ap a0 
. 


_ £07) 


ne cos 0(27 sin @ dé) 


= Po(7) 

4° _ 

It can be seen that the total flux (7), which 

corresponds to a 47 solid-angle integration of the 

angular flux, is exactly four times as large as the 
corresponding positive partial current. 

In many real situations the angular flux, although 
anisotropic, can be regarded as being isotropic in the 
forward half-space (or in 2a angular directions). For 
this case the positive partial current would have 
essentially the same magnitude as a positive partial 
current in an isotropic flux field would have. As a result 
it has become standard practice in shielding calculations 
to define a new total flux, referred to as the 27 total 
flux* and given by 


(2.54) 


©27(7) = J #(7,9) dB. (2.55) 
where the 27 solid-angle integration is accomplished in 
the forward half-space. Since for the isotropic flux field 
the 47 total flux is exactly two times the 27 total flux, 
the positive partial current can be expressed as 
I") = 93 7(9/2. (2.56) 
From the above it is obvious that care must be 
exercised in applying this formula, which was derived 
on the basis: of a J, MO) source condition, since either the 
2n total flux 2*(r) or the 47 total flux ® 07) may be 
specified for a typical problem, The corresponding 
partial current would be either ©27(r)/2 or ®,(r)/4, 
and a factor-of-2 error could result unless it is clear 
which total flux is specified. 


*The 27 total flux is also called the incident flux. 


2.2.3. POINT SOURCE 


The general volumetric source term in seven-dimen- 
sional phase space is defined as 


S(r,E,Q,t) dr dE dQ =the number of particles having 
energies in dE about E am direc- 
tions in dQ about Q that are 
emitted per unit time at the time 
t from the differential volume dr 
at the space point 7. 


A useful specialization of the general source term is the 
angular point source, which is given by 


So(ro,E,2,t) dE dQ =the number of particles having 
energies in dE about £ and direc- 
tions in dQ about Q that are 
emitted per unit time at time f 
from a point source located at 
the space point Fo. 


The angular point source is related to the general source 
term by the identity 


So(Po.£,2,t) 
=f f S(E,Q,1) 6(Q,2’) dr dQ’, (2.57) 
4n “all space 
where the general source term is specialized for the 


point-source configuration by inclusion of the Dirac 
delta function 


S(7,E,Q',t) = 8(r — 7) So(E,2',1) , (2.58) 


and So (E,0',1) is an energy-, angle-, and time-dependent 
source intensity factor. For the usual problem the 
assumption of energy-angle separability is permissible, 
and the source intensity factor can be expressed as 


So(E,Q,t) = So(t) FE) (2) , (2.59) 


where So(ft), the source intensity factor, is 


sone f [> So0E8,0 ak ad, (2.60) 
4n “0 
and the other terms are defined as 


5(rF — ro) = the Dirac delta function which locates the 
point source at 7o, 


F(E) = normalized energy spectrum, 


(2) = normalized directional probability density 

function having the property f an G2) dQ 
=] . 

g(Q) dQ = fraction of particles having directions in 
aQ, about Q, 

§(9,9') = angular delta function which aligns the 
directions of the emitted particles with 
the direction 2, 


Q= direction vector which lies along the 
vector ¥ — Fp as shown in Fig. 2.6. 


Evaluation of the integrals in Eq. 2.57 shows that 


So(Fo,E,Q,t) = So(t) F(E) g(Q) . (2.61) 


For isotropic emission g(Q) dQ = dQ/4r, in which case 


se 


So(Fo, E, es t= re F(E) . (2.62) 
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The differential point source is obtained by integrating — 


Eq. 2.61 over all directions. Thus 


So(ro,£,t) = S(t) F(E) : (2.63) 
Similarly, the total number of particles emitted per unit 
time at time ¢, S(t), is obtained by integrating the 
differential point source over all energies. Finally, if 
total emission of particles is desired, the time- 
dependent source intensity factor So(z) is integrated 
over some specific interval of time. It should be noted 
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Fig. 2.6. Point Source Specified in General 7 Space. 
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that the various mathematical descriptions of the point 
source presented here apply to both steady-state and 
time-dependent source conditions. For the steady-state 
problem the phase space notation becomes (7,£,), and 
an emission rate constant in time is presumed. 

Many source configurations can be represented by a 
point source or a simple arrangement of point sources; 


for example, isotropic line and plane sources consist of 


appropriate distributions of isotropic point sources. 
Also, the kernel technique, known as the method of 
Green’s functions, has for its basis the point kernel, the 


point kernel being the solution to a problem with a unit - 


point source. In summary, the point source is a most 
elementary and versatile source configuration and has 
considerable application to complex source configura- 
tions such as the detonation of a nuclear weapon. 


2.2.4. ISOTROPIC PLANE SOURCE 


An isotropic plane source, depicted in Fig. 2.7, can be 
regarded as an arbitrary distribution of unit isotropic 
point sources, and in six-dimensional phase space is 
defined as 


S(xo y,z,E,Q) dE dQ =the number of particles having 
energies in dE about EF and 
directions in dQ about Q that 
are emitted per unit time from 
a unit area located at xo,y,z and 
normal to the direction vector n 


x 


= S(y,2) dE dO , (2.64) 
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Fig. 2.7, Orientation of Isotropic Plane Source with Respect 


to the Rectangular Coordinate System. 
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where 
Xo = location of the source plane along the x axis, 


S(y,z) = planar distribution of unit isotropic point 
sources, 


FE)_ ; aries 
de = unit angular point source energy distribution 
for isotropic emission. 


For an isotropic plane source having a uniform 
distribution of point sources, S(y,z) = S (particles cm ~ 
sec '), and the expression for the source becomes 


S(xo,E,Q) dE dQ. = s75) dE dQ. (2.65) 
T 


This source is of considerable importance in shielding 
work; for example, it is used to represent radiation 
fields such as those due to radioactive fallout. 

Note that S(xo,E,Q) is based on a unit area that is 
normal to the x axis. The angular flux that corresponds 
to the isotropic plane source has for its basis a unit area 
normal to the Q direction and therefore differs from 
S(xo,E,Q) by a factor of (Qn)! or cos 6 with 


azimuthal symmetry; that is, 


B(x, £,Q) =(Q:n)! S(xo,E,Q) 


_ of) 1 
8 Ep 268) 


The utility of ®(x9,£,Q) is sometimes limited by the 
singularity experienced as (Q:n) approaches zero. The 
term S(xo,£,9) has proved to be the more useful form, 
and integrations over direction Q and energy E lead to 
physically meaningful quantities. 

The differential form for the plane source is obtained 
by integrating Eq. 2.64 over all directions, 


S(x,E) = f _ Slo £2) dB =5 FUE) , (2.67) 


and the total emission (particles cm~* sec™!) is 
obtained by the subsequent integration over all en- 
ergies, 


Siro) = f é S(o,E) dE =Sp. (2.68) 


2.35 


In describing the transport of radiation through 
matter, the interactions that the various particles may 
have with the nuclei of the medium must be considered. 
The probability of the occurrence of a given inter- 
action, which will vary with the type of particle, its 
energy, and the material it traverses, is usually given in 
terms of microscopic or macroscopic cross sections. 
These cross sections are described in the following 
paragraphs, along with their relation to reaction proba- 
bilities and reaction rates. 


2.3.1. MICROSCOPIC CROSS SECTIONS 


The microscopic cross section, designated by the 
symbol o, can be considered to be the effective area that 
a nucleus presents to an oncoming particle, regardless of 
the particle direction. That is, a particle will interact 
with a nucleus with a probability of unity when the 
particle is incident on the area and with a probability of 
zero when it misses the area. The units usually 
associated with o are cm?/atom or barns/atom, where 1 
barn* = 10°?* cm?. Because of the wave-like behavior 
of the particles and the complex nature of the nucleus, 
the magnitude of o in general bears no direct relation 
with the physical size of the nucleus, but is strictly a 
function of the properties of the atom and the energy 
and type of the interacting particle. 


2.3.2, REACTION PROBABILITIES 


If a nucleus having a cross section o is contained 
within some characteristic volume of material, V,, 
which has a projected area A, then the probability that 
a particle incident on the volume will interact with that 
particular nucleus is the ratio of the effective area o to 
the projected area A, (see Fig. 2.8a). It follows then 
that o/A, is the reaction probability per nucleus 
associated with the area A, and that it can be 
interpreted as the average number of reactions oc- 
curring in V, per incident particle. 


*The term barn is attributed to have resulted from a humor- 
ous remark to the effect that 10-24 cm? looked as big as a barn 
to a neutron. 
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Quantities Used to Describe Radiation Interactions 


If the nuclide density NM (nuclides/cm*) is uniform 
throughout V,, the total number of nuclei contained 
within V,, is given by NV, and the reaction probability 
associated with the total number of nuclides is equal to 
the linear sum of the reaction probabilities associated 
with the individual nuclei provided that the character- 
istic volume V, is optically thin (no self shielding). The 
reaction probability now associated with the character- 
istic volume V, is NV (0/A,). 

This summation of probabilities is shown schemat- 
ically in Fig. 2.8b. The optically thin condition is 
assured if the characteristic volume V, is differential in 
size. (The characteristic area A, would also be differ- 
ential in size.) Otherwise, the condition shown in Fig. 
2.8¢c might result, in which case the effective cross 
section per nucleus would be decreased because of the 
self-shielding effect. This follows since the reaction 
probability would no longer be a simple linear sum- 
mation of the individual probabilities. This can easily be 
seen for the example of an opaque volume; increasing 
the thickness has no further effect. 


2.3.3. MACROSCOPIC CROSS SECTIONS 


The reaction probability associated with the optically 
thin characteristic volume V, can be rearranged and 
written as (ON)(V,/A,). If V,/A, is recognized as the 
characteristic differential distance dS, traveled by 
particles within V, and if oN is identified as the 
macroscopic cross section X,t then ZdS c is the reaction 
probability associated with the differential tracklength 
dS,. It follows that a formal definition for the 
macroscopic cross section 2 is the reaction probability 
per differential tracklength, or, presuming that differ- 
ential tracklength is implicit in 2 and interpreting 
reaction probability as reactions per incident particle, 
the definition for 2 becomes the number of reactions 
per particle tracklength. However, the units for 2 are 
usually given as reciprocal centimeters (cm~'). 


TThe symbol u is used for gamma rays. 


© 


c 


There is another approach in-defining the macro- 
scopic cross section. In particle transport problems large 
numbers of particles are involved, and so it is sometimes 
convenient to describe the process in terms of averages 
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Fig, 2.8. Relationship Between Microscopic Cross Sections - 


and Reaction Probabilities. 

(2) A nucleus having a microscopic cross section o is 
contained within a characteristic volume of material V, which 
has a projected area A,. A particle is entering V, and the 
probability that it will interact with the nucleus is given by the 
ratio o/A,. This ratio can be interpreted as the number of 
reactions per incident particle per nucleus. 

(b) Several nuclei, each having a microscopic cross section a 
are contained within the volume V,, which is sufficiently thin 
that the o’s will not overlap. The probability that the incident 
particle will interact with a nucleus in V, is given by 
(a/A oNV¢, where N is the number of nuclides per em? in Vig: 

(c) The characteristic volume V, is sufficiently thick for the 
a’s to overlap. In this case the self-shielding effect would 
prevent the reaction probability from being a simple linear 
summation of the individual probabilities. 


or expected values. It can be assumed that the particle 
has a certain probability of undergoing an interaction 
with a nuclide within a distance dx. This average or 
expected value, called the macroscopic cross section 
Z in units of inverse length, depends on the particle 
energy and the density and nuclear properties of the 
medium. Thus in penetrating a distance dx in a 
medium, a particle will suffer Xdx average number of 
interactions per unit time. 

In order to examine the concept more closely, 
consider the simple model of a beam of particles 
normally incident on an infinite slab of some material. 
Let ®(x) be the uncollided flux at thickness x. The 
number being removed in dx is 


—d®(x) = &(x)Edx , (2.69) 


which integrates to 
@(x) = Bye =* , (2.70) 


where ®g is the incident flux. The probability that a 
particle will suffer an interaction in a distance less than 
Lis 


> 


L YP(x)d. L 
Pa <L)= f ef, De =*dx=1 —e- FL 


(2.71) 


where Le =~ is the probability density function for an 
interaction. The mean distance a particle will travel (the 
mean-free path )) before it suffers an interaction is 


1 
= —=x = 5S 
r f xZe dx =: (2.72) 


Af the foregoing discussion is accepted as establishing 
the definition of a macroscopic cross section, which is a 
function of both the properties of the nuclide (nucleus 
only for the neutron) and the particle energy and type, 
the microscopic cross section o can be defined as simply 
=/N, where N is the nuclear density and o represents 
the cross-sectional area per nucleus seen by the incident 
particles. 

It should be clear on the basis of the foregoing dis- 
cussion that 2 should not be used as a probability 
density function to determine reaction probabilities 
except in special cases. Even though the reaction proba- 
bility for a particle normally incident on a differential 
volume dx within the slab is given by Zdx, the 
probability that a reaction will occur within the 


optically thick slab is not simply f ” Sdx = ZL. As was 


shown above, the reaction probability for the slab is 
given by (1 — e— 2), which approaches ZL only for 
the special case of an optically thin slab since (1 — 
e-=L) > SL for SL <1* 


2.3.4. REACTION RATES 


The concept of a macroscopic cross section as 
described in the preceding section suggests that particles 
which move through matter will undergo interactions at 
a rate proportional to the total particle tracklength, a 
quantity that is closely related to the particle flux &(r) 
(see Table 2.1). The macroscopic cross section 2 is in 
fact the required proportionality constant between the 
reaction rate and the particle flux; the reaction rate per 
unit volume at the space point 7 is given by 

R(*) = TH) (2.73) 
the dimensional correctness of which is verified by 
noting that 


reactions reactions particle-cm 


cm?-sec_ particlecm cm?-sec ~ 


The differential reaction rate is related to the differ- 
ential particle flux and is given by 


RGE) = X(E) Hr,E) , (2.74) 


where 


RE) dE = number of reactions at space point 7 per 
unit volume and time which involve 
particles having energies in dE about the 
energy £, 


*Another and more definitive description of the optically 
thin condition. 
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X(£) = macroscopic cross section for the reac- 
tion, which depends on the particle’s 
energy as well as the nuclear character- 
istics of the medium, 


@(r,£) = differential particle flux density. 


The reaction rate for an energy interval, or group 
(AEg = Ey — Eg+1), is simply related to the differential 
reaction rate by 


Ro fn _ ROE AE 
= if 2) O(7,E) dE. (2.75) 


This leads to a definition of the reaction rate in terms 
of an energy-averaged macroscopic cross section 2, 
(group cross section) and the group flux: 


Ro) = 26 G(r), (2.76) 
where 
Spf XE) OE) dE 
lige - (2.77) 
i F,E) dE 
and 
(2.78) 


PA=fe%  UEB)AE. 


The concept of the reaction rate in its many forms 
is fundamental to the formulation of the various 
equations that analytically describe particle transport 
(for example, the Boltzmann transport equation) and to 
the determination of useful information, such as energy 
deposition and biological dose from known or calcu- 
lated values of the particle flux density. 
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2.4. Radiations Produced in Nuclear Weapons 


As pointed out in the introduction to this chapter, 
several types of radiations are produced by nuclear 
detonations, but generally the only radiations that the 
shield designer must be concerned with are neutrons 
and gamma rays. It is the nature and classification of 
these radiations that are discussed in the following 
paragraphs. The processes whereby neutrons and 
gamma rays interact with materials are discussed in 
Sections 2.5 and 2.7 respectively, and the processes 
whereby they are produced are described in Sections 
2.6 and 2.8 respectively. 


2.4.1, NEUTRONS 


The presence of a highly penetrating electrically 
neutral radiation had already been experimentally 
demonstrated* by several groups of investigators when 
Chadwick? in 1932 established the existence of a 
neutral particle having a mass nearly equal to that of a 
proton. Since then, this neutral particle, which has 
become known as the neutron, has proved to be not 
only a key to nuclear weaponry but also a bane to those 
persons charged with providing protection against 
nuclear radiation fields. 

Neutrons, like protons, are constituents of the 
nucleus, and both are referred to as nucleons. The 
ability of a neutron to transform into a proton — and 
the converse — suggests that protons and neutrons may 
be simply different states of a single particle, and a full 
understanding of the relation between the two is one of 
the current problems in nuclear physics. Neutrons can 
be regarded as fundamental particles since they are 
readily identifiable by the following physical properties: 


Charge 0 

Rest mass 1.00898 atomic mass units 
Spin Y% 

Magnetic dipole moment —1.913 nuclear magnetons 
Statistics Fermi 

Intrinsic parity Even 


*This radiation was produced when beryllium was bombarded 
by alpha particles, 


Outside a nucleus the neutron is not a stable particle 
and decays with a half life of 12.8 min. However, this 
effect is always neglected in neutron transport problems 
because of the extremely short transport lifetimes (in 
general much less than 1 sec). 

The theory of wave mechanics suggests that all 
moving particles, including neutrons, have wave proper- 
ties. The deBroglie wavelength //p can be assigned to a 
neutron, and if its kinetic energy is less than the energy 
equivalent of its rest mass,t that is, if & <939.5 MeV, 
the nonrelativistic wavelength is given by 


2.8600 X 10712 


TIT : (2.79) 


eee 
Pp 


where 


d = deBroglie wavelength (cm), 
£ = neutron energy (MeV), 


h=Planck’s constant (6.624 X 10°?” erg:sec, or 
4.132 X 10715 eV‘sec), 


p = momentum of the neutron. 


The neutron’s wavelength is inversely proportional to 
its velocity and provides valuable insight into the 
general nature of neutron interactions with matter. 
Thermal neutrons have wavelengths of the same order 
as the atomic spacings in crystals. Therefore, when 
thermal neutrons interact with a crystal, diffraction 
patterns in accordance with Bragg’s law are produced. 
(The main difference between the diffraction of neu- 
trons and of gamma rays is that neutrons interact with 
the nuclei whereas gamma rays interact primarily with 
the orbital electrons.) Neutrons with energies from 1 to 
100 MeV have wavelengths of the order of nuclear 
dimensions, and neutrons with relativistic energies 
greater than 10 GeV have wavelengths comparable with 


tThe conversion factor that relates a particle’s rest mass in 
atomic mass units (amu) to the energy equivalent of its rest 
mass is 931 MeV/amu. The mass of a neutron is 1.00898 amu, 
which when multiplied by 931 MeV/amu yields 939.5 MeV for 
its energy equivalent. 


the spacing of nucleons within the nucleus, which 
provides another tool for exploring the structure of the 
nucleus.* 

Neutrons are classified according to energy, with the 
divisions occurring somewhat naturally because of 
behavior peculiar to particular energy ranges. In some 
cases the divisions are somewhat arbitrary, and they can 
be ambiguous since the energy structures overlap and 
differences in opinion exist as to the proper termi- 
nology. However, classification is helpful in discussing 
neutron behavior. Some of the classifications used are 
given below. 

Slow Neutrons. Neutrons with energies from zero to 
about 1 or 2 keV are sometimes referred to as slow 
neutrons. This is a broad classification, however, and is 
not in general use since it does not serve as a very 
precise description of the energy state of low-energy 
neutrons. In reactor physics and technology jargon, 
the term “slow neutrons” refers to thermal neutrons. 

Thermal Neutrons. Thermal neutrons are in thermal 
equilibrium with the atoms (or molecules) of the 
medium in which they are present. A particular thermal 
neutron may gain or lose energy in any one collision 
with the nuclei of the medium, but there is no net 
energy change for a large number of neutrons diffusing 
in a nonabsorbing medium. For this condition the 
kinetic energy of the neutrons is distributed statistically 
according to the Maxwell-Boltzmann equation as de- 
rived from the kinetic theory of gases. The equation for 
the spectral distribution is 


dn 2mm E 
= = eg E/RT 71/2 | 2.80 
mE) =e CakTyl? © G20) 


where 
dn = n(£E) dE = number of neutrons having energies in 
the energy interval dE about £, 
n= total number of neutrons in the system, 
F = neutron energy (eV), 
k = Boltzmann constant (8.61 X 10~* eV per °K), 
T = absolute temperature (°K). 


It can be shown that, when the Maxwell-Boltzmann 
distribution applies, the energy corresponding to the 
most probable neutron velocity (2200 m/sec at 20°C) is 
given by kT (0.0252 eV), the average kinetic energy of 


*Additional details on the properties of neutrons are available 
in most texts on nuclear theory (see, for example, a book by 
Evans’). 
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thermal neutrons is *%4kT, and the most probable 
neutron energy is 44 KT. 

For absorbing media the true energy distribution of 
the neutrons deviates from the Maxwell-Boltzmann 
distribution because neutron absorptions tend to in- 
crease with decreasing energy, which causes an increase 
in the average energy.* This tendency is illustrated in 
Fig. 2.9. The spectral distribution shown in Fig. 2.9a 
exhibits very little distortion, which is characteristic of 
poorly absorbing media, while the distribution shown in 
Fig. 2.9b clearly demonstrates the effects of the 
preferential absorptions at the lower neutron energies. 
This phenomenon is called hardening and can in part be 
compensated for by the use of an “effective neutron 
temperature” 7,,, which is defined as a number that, 
when used in the Maxwell-Boltzmann distribution, gives 
the best least-squares fit to the true distribution as 
determined by a Monte Carlo calculation. The effective 
neutron temperature 7,, has been found to be approxi- 
mately related to the environmental temperature T,, in 
the following way:* 


1.84 2,(k7) 
Tet, [ + ae =T, [1+0.46A] , (2.81) 


where 


s 


2,(kT) = macroscopic absorption cross section evalu- 
ated at the energy KT, 


2, = macroscopic scattering cross section for 
thermal neutrons, 


£= average logarithmic energy decrement per 
collision, 


42 ,(KT) 
Ex ; 


s 


Cold Neutrons. — Neutrons with an average energy 
less than that of thermal neutrons are sometimes 
referred to as cold neutrons. They are of no concern in 
shielding calculations. 

Epithermal Neutrons. — Epithermal neutrons are the 
neutrons in a system whose energies exceed those of a 
Maxwell-Boltzmann distribution for the effective tem- 
perature of the neutron-supporting medium. 

Resonance Neutrons. — Neutrons having energies in 
the range between 1 and 100 eV are frequently called 
resonance neutrons. Many nuclei exhibit strong absorp- 
tion of neutrons at well-defined energies in this energy 
range, and these absorptions are commonly referred to 
as resonance absorptions. 
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Fig. 2.9. Comparison of the Neutron Energy Disbributions with Maxwell-Boltzmann Distribution in (¢) Poorly Absorbing 
Medium and (+) Highly Absorbing Medium, (From Coveyou et al, ref. 4.) 


Intermediate Neutrons. — Neutrons in the energy 
range 1 keV to about 0.5 MeV are often referred to as 
intermediate neutrons. The slowing down and removal 
aspects of neutron transport primarily involve inter- 
mediate neutrons, and their behavior strongly affects 
the design of high-performance shields. 

Fast Neutrons. — Neutrons having energies between 
0.5 MeV and about 15 MeV are regarded as fast 
neutrons. Most neutrons produced by fission and fusion 
reactions are fast neutrons. 

Very High Energy and Ultra High Energy Neutrons. — 
Neutrons referred to as having very high energy usually 
fall in the 15- to 50-MeV range and those referred to as 


having ultra high energy fall in the range > 50 MeV. 
These neutrons are either cosmic radiations or by- 
product particles produced during accelerator opera- 
tions. They are not of interest in weapons radiation 
shielding, but are included here to complete the 
neutron classifications. 


2.4.2 GAMMA RAYS 


During early investigations of radioactive materials 
three types of radiation were observed that were 
arbitrarily identified as alpha, beta, and gamma rays. 


Gamma rays, unlike the alpha and beta rays, which 
were found to carry positive and negative charges 
respectively and to be easily stopped, were found to be 
electrically neutral and highly penetrating. It was 
postulated that gamma rays were actually electro- 
magnetic radiations similar to X rays, and this was 
confirmed experimentally in 1914 by E. Rutherford. 

According to classical theory, electromagnetic radia- 
tion consists of transverse electric and magnetic fields 
that propagate as waves at the speed of light. They may 
exist in regions that are void of matter and contain no 
electric charges or currents. Within the framework of 
classical theory a continuous flow of energy is associ- 
ated with the propagation of the electromagnetic wave, 
and there is no upper or lower limit to the amount of 
energy that can be transported by electromagnetic 
radiations. 

While wave theory gives a good description of the 
diffraction, reflection, and interference of electro- 
magnetic radiations, it is quite inadequate for describing 
phenomena such as nuclear emission, the photoelectric 
effect, Compton scattering, and bremsstrahlung.* How- 
ever, the quantum theory, which was formulated by 
Max Planck in 1900, does provide a satisfactory 
explanation for such phenomena. This theory suggests 
that electromagnetic radiations are emitted as discrete 
bundles of energy of size hv, where Planck’s constant h 
is equal to 6.624 X 107?” ergs/sec and vp is the 
frequency. These pseudo-particles are commonly re- 
ferred to as photons. They have energies of E = hv and 
possess the properties of momentum (p = hp/c) and 
relativistic mass (M = hv/c’). 

Wave theory and quantum theory seemingly con- 
tradict each other in many ways, but in 1928 Max Born 
argued that electromagnetic radiation was indeed 
quantum in nature and that the classical wave theory 
could be interpreted as the probability of finding a 
photon at a particular point in space and time. It is 


*These phenomena are discussed in Section 2.7. 
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within the context of quantum theory that electro- 
magnetic radiations are treated in this chapter. 

In the years since 1914 it has been established that 
electromagnetic radiations originate whenever transi- 
tions take place between two energy levels of nuclei or 
between two electronic levels of an atom. The fre- 
quencies and wavelengths of the radiations are de- 
termined by 


T=hv=he/, (2.82) 


where the level width [" is the difference in the | 


excitation energy of the two levels. Since the level 
widths are greater for nuclear transitions than for 
electronic transitions, the radiations accompanying nu- 
clear transitions are characteristically of higher energy 
than those accompanying electronic transitions. For 
historical reasons the radiations produced by nuclear 
transitions are usually referred to as gamma rays and 
those produced by atomic transitions are called charac- 
teristic X rays. Typical of the latter are the soft X rays 
(fluorescence radiation) produced as a result of the 
photoelectric effect discussed in Section 2.7. Two other 
types of electromagnetic radiations are bremsstrahlung 
and annihlation radiation. Bremsstrahlung radiation 
arises as a result of the acceleration or deceleration of 
free charged particles, and annihilation radiation is 
produced when a positron and an electron combine, 
both particles being annihilated in the process. (Posi- 
trons are produced by the pair-production process 
described in Section 2.7.2). 

From the foregoing it is seen that electromagnetic 
radiations are classified by origin; that is, (1) gamma 
rays accompany nuclear transitions, (2) X rays are 
emitted in atomic transitions, (3) bremsstrahlung radia- 
tion results from the acceleration or deceleration of free 
charged particles, and (4) annihilation radiation is 
emitted when a positron combines with an electron. As 
it turns out, and as is discussed in Section 2.7, 
bremsstrahlung and annihilation radiations, as well as 
soft X rays, are not usually highly penetrating and the 
transport of these radiations is usually ignored in 
shielding design. 
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2.5. Interactions of Neutrons 


Neutrons are attenuated by interactions with the elastic spheres (sometimes referred to as billiard-ball 
nuclei of the elements comprising the medium that they collisions) and all other interactions. The total cross 
are traversing. The neutron interactions most important — section can also be thought of as the sum of the 
to shielding are elastic scattering, inelastic scattering, total scattering cross section and the total absorption 
radiative capture, and charged-particle reactions. The cross section. In the latter interpretation the scattering 
relative importance of these interactions is, at least in cross section includes both elastic and nonelastic 
part, determined by their respective cross sections, the scatterings, and the absorption cross section includes 
magnitude of which varies with the neutron energy and _. all processes ‘that cause the product nucleus to 
the element (see Section 2.3). differ either in atomic number or in atomic weight from 

The relation of the various cross sections is shown in the target nucleus. Because of the change in the target 
Fig. 2.10. From this figure it is apparent that the total nucleus, absorption processes are usually referred to as 
cross section for a neutron interaction can be thought _ nuclear reactions rather than interactions, although the 
of as being the sum of the elastic-scattering cross two terms are often used interchangeably. 
section and the nonelastic cross section. The two The accuracy (or reliability) of a given calculation can 
classifications distinguish between those interactions in be strongly influenced by the quality of the cross- 
which the neutron and nucleus collide as two perfectly section data used. Despite their obvious importance, 
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detailed neutron cross sections are not presented in this 
document because of the vastness and complexity of 
the information. Much cross-section data are readily 
available in a number of publications.*° Also available 
is a comprehensive and periodically updated index'® to 
the latest literature on microscopic neutron data. 


2.5.1. ELASTIC SCATTERING 


As mentioned above, the elastic scattering of neutrons 
with nuclei can be viewed as the collision of perfectly 
elastic spheres.* This process is characterized by the 
conservation of kinetic energy and momentum and is 
exactly described by the methods of classical me- 
chanics. The various theoretical analyses are usually 
performed in the center-of-mass (c.m.) frame of ref- 
erence, which assumes that the center of mass of the 
neutron-nucleus system remains at rest. An external 
observer (e.g., one performing a physical measurement) 
would view the collision process in the laboratory frame 
of reference. Theoretical results are more useful when 


expressed in terms of the laboratory system coordi- 
nates, but the theoretical relations are more easily 
derived in the c.m. system. 

If it is assumed that the target nucleus in the labora- 
_ tory system is stationary, the change in kinetic energy 
experienced by the neutron is uniquely related to its 
' change in direction by the following equation: 


E= Eo ain oes 


(A+1 G22) 


where 


Eo = initial neutron energy in the laboratory system, 
£ = final neutron energy in the laboratory system, 


A =ratio of the atomic mass of the nuclide to the 
mass of the neutron, 


n=cos 0, the cosine of the scattering angle in the 
c.m. system. 


The minimum value of £/F, which corresponds to the 
maximum loss in energy experienced by a neutron ina 
single collision, occurs for a head-on collision, @ = 7 or 
n = —1. The minimum final neutron energy is given by 


*This description applies only to “potential” elastic scatter- 
ing, which is the type discussed in this section. Another kind of 
elastic scattering, called capture scattering, is very much like 
inelastic scattering in that it involves the formation of a com- 
pound nucleus and the subsequent emission of a neutron with 
nearly its original kinetic energy (see later discussion). 
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(2.84) 


=[(A —1)/(A+1)]?. (2.85) 
The cosine of the scattering angle in the laboratory 
system, symbolized by y, is simply related to the cosine 
of the scattering angle in the c.m. system by 


1+7A 


Gioaeie (2.86) 


B= 


If it is further assumed that the scattering in the c.m. 
system is isotropic (an excellent assumption for hy- 
drogen), then the following relations hold for the 
elastic-scattering process: 

1. The probability that a neutron will undergo a 
change of direction into a differential solid angle dQ 
corresponding to a conical element defined by the 
scattering angles @ and @ + dé in the c.m. system is 


dQ 1 
p(@) dé ae =5 sin 6dé. (2.87) 


2. The fraction of all elastic-scattering collisions that 
involve neutrons of energy Eo scattered into the energy 
interval dE about the final energy £ is given by 


dE 


P(E;Eo) dE "FG a) : 


(2.88) 


3. The average final energy of all neutrons of initial 
energy Eo is given by 
= E 
E= 
e 


° E p(E-Eo) dE = [(1 + 0)/2] Ey. (2.89) 
i?) 


4. The average value of the cosine of the scattering 
angle in the laboratory system, an important parameter 
in neutron transport calculations, is given by 


aes, (2.90) 


5. The average logarithmic energy decrement per 
collision is defined by 


In = p(EsEo) de = 1 SE 


(2.91) 
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An interesting property of & is its independence of the 
initial energy of the neutron. The quantity & can also 
be interpreted as the average change in lethargy per 
collision, where the lethargy u is defined as 


(2.92) 


and the average change in lethargy per collision is given 


aed (2.93) 


The lethargy is a useful alternate energy variable in the 
description of the slowing-down process for neutrons. 

For most nuclides the isotropic assumption generally 
gives a good description of the elastic-scattering process 
for neutrons in the thermal-to-keV energy range; 
however, it is invalid for high-energy (MeV) neutrons or 
for interactions in which chemical binding energies are 
important. Thus in shielding calculations, which usually 
involve high-energy neutrons, the isotropic assumption 
cannot be made safely except for hydrogen, for which 
the scattering is virtually isotropic in the c.m. system 
for neutron energies up to 14 MeV. 

In a more general treatment of the scattering process, 
which can apply to inelastic as well as to elastic 
scattering, the cross sections for scattering are repre- 
sented by an expansion in a series of Legendre 
polynomials: 


o,(E) y 2m+1 


E,n) = 
0,(E,n) ae ; 


Am(E) Pan), (2.94) 


where 
o,(E,n) = energy-dependent differential cross section 
for the elastic scattering of neutrons of 
energy £& through an angle in the c.m. 
system whose cosine is 7 per unit solid 
angle, 


A, (E) = energy-dependent Legendre coefficients 


(except that Ag = 1 always), 
P,,(n) = Legendre polynomials. 


The energy-dependent elastic-scattering cross section 
0,(£) is obtained by integrating the. differential elastic- 
scattering cross section 0,(£,7) over all possible values 
of 7: 

+1 
o,(E) = 20 I o,(E,n) dn. (2.95) 
Representation of the cross section as an expanded 


series of orthogonal functions is necessary in the 
spherical harmonics and S,, solutions of the Boltzmann 
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equation (see Chapter 3). Other orthogonal functions 
can be used — and have been — but the Legendre poly- 


-nomials generally give best fits and are universally used. 


Note that the first term (zeroth order) of Eq. 2.94 
gives the angular cross section for isotropic scattering 
a,/4n. The coefficients A,,(E) obtained by fittine 
angular cross-section data are reported in both labora- 
tory and c.m. systems. Usually the c.m. coefficients 
must be transformed to the laboratory system for 
transport calculations. As the atomic number increases, 
however, the difference becomes small. 

The elastic-collision process is an important energy 
degradation mechanism and is most important for the 
light elements, particularly hydrogen. Efficient neutron 
shields always contain hydrogenous materials. 


2.5.2. NONELASTIC REACTIONS 


All reactions that are not elastic scatterings are called 
nonelastic reactions, and the sum of their cross sections 
is referred to as the nonelastic cross section. In 
experimental determinations of neutron cross sections 
in the energy range 3 to 15 MeV, the nonelastic cross 
section is the quantity that is actually measured. Or, in 
other words, the quantity measured is the total cross 
section minus the elastic-scattering cross section. 

Many models are available for describing nonelastic 
neutron reactions with nuclei, but the one most 
frequently used is the Bohr model. In this model the 
interaction is envisioned as beginning when the neutron 
comes within the range of the nuclear forces of the 
nucleus and as terminating when the products of the 
interaction depart from this region. In the interim the 
combination of the neutron and target nucleus is called 
simply the compound nucleus. It is convenient to 
consider the interaction as occurring in two distinct 
stages: 

1. A compound nucleus is formed when the incident 
neutron loses its identity and becomes assimilated into 
the unstable new system. This event brings to the newly 
formed compound nucleus an excitation energy that 
consists of the available portion* of the neutron’s 
kinetic energy of motion and its binding energy} with 


*Due to conservation of momentum, some of the kinetic 
energy of the incident neutron must appear as translational 
energy at the compound nucleus and is therefore unavailable as 
excitation energy. 

1The binding energy per nucleon is the average energy 
required per nucleon to separate the nucleus into its constituent 
particles, Refer to subsequent discussions of binding energy 
under “Nuclear Reactions.” 


respect to the compound nucleus. The excitation 
energy is statistically shared by all the nucleons. 

2. The unstable compound nucleus separates into its 
products when one or more of the nucleons acquire 
sufficient energy to escape, that is, when the nucleons 
acquire energy that is equal to or greater than their 
binding energy. In some cases the unstable compound 
nucleus may give off excitation energy in the form of 
gamma radiation. In either case the resulting product 
nucleus may itself be unstable and undergo radioactive 
transformation, but in all cases the end result is a 
system of stable products. A schematic representation 
of the two stages of neutron interaction with a nucleus 
is shown in Fig. 2.11. (The example given in the figure 
is specifically for the inelastic-scattering process de- 
scribed below.) 

The compound nucleus concept is being replaced by 
more elaborate models; however, it will continue to be 
useful in providing an easily understood representation 
of nonelastic neutron interactions with nuclei. 

The nonelastic neutron reactions are divided into 
scattering processes and absorption processes, the ab- 
sorption processes often being referred to as nuclear 
reactions. They include radiative capture, charged- 
particle reactions, (”,2n) reactions, and fission. Insofar 
as shielding attenuation is concerned, the (n,2n) and 
fission reactions are usually unimportant unless ma- 
terials such as beryllium and depleted uranium are used. 
The fission reaction is, of course, very important as a 
neutron source in weapons detonations and as the 
energy source in the production of electric power (see 
Section 2.6). 
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Fig. 2.11. Schematic Representation of a Nonelastic Inter- 
action of a Neutron with a Nucleus (Inelastic Scattering). 
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Inelastic Scattering 


In inelastic scattering the incident neutron and target 
nucleus form a compound nucleus from which a 
neutron is subsequently emitted. If the kinetic energy 
of the incident neutron is greater than the lowest 
energy level of the target nucleus,* the nucleus may be 
left in an excited state following the neutron emission. 
This results in the emitted neutron having substantially 
less energy than the incident neutron, thereby con- 
tributing strongly to the energy degradation process and 
therefore to the attenuation of neutrons in shields. The 
excited nucleus then returns to its ground state by the 
emission of gamma rays. An energy level diagram of this 
process is shown in Fig. 2.12. 

Determining the detailed cross sections for inelastic 
scattering is a complicated and difficult task, whether 
accomplished by experiment or by theory. Conse- 
quently, the inelastic-scattering cross sections for many 
elements are unknown and a large fraction of those 
used in shielding calculations are little more than a best 
guess. However, there are some generalizations that can 
be made: 

1. The cross section usually increases with increasing 
neutron energy. 

2. The cross section usually increases with increasing 
atomic weight. 

3. The threshold for the interaction in light elements 
is generally several MeV, and the gamma rays emitted 
are correspondingly harder than those produced in 
interactions with lower energy thresholds. 

4. The threshold for the interaction in heavy elements 
is lower, sometimes much lower, than for that in light 
elements. 


Nuclear Reactions 


The term nuclear reaction is often reserved for those 
interactions in which the product nucleus differs either 
in atomic number or atomic weight from the target 
nucleus. The type of reaction that occurs depends on 
several factors, the most important being the nuclear 
structure of the target nucleus and the energy of the 
incident neutrons. Some reactions are possible for 
neutrons of any energy, while others have threshold 
energies. In some cases more than one de-excitation 
process is possible, with the relative probabilities of 
each being given by their respective reaction cross 
sections. 


*If this condition is not satisfied, the capture scattering form 
of elastic scattering may occur. 
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Fig. 2.12. Energy Diagram for the Inelastic Scattering of a Neutron with a Target Nucleus. 

{Given by 931 MeV [(My + M,) — My] ~7 MeV, where My, M,, and My are the masses, in atomic mass units, of the target 
nucleus, the neutron, and the compound nucleus, respectively. 

bGiven by KE, (Mx/(My + M,,)]. Due to the conservation of momentum, some of the kinetic energy of the incident neutron 
must appear as translational energy of the compound nucleus and is therefore unavailable as excitation energy. 

“Given by 931 (My — Mx) ~931 MeV. 

4Given by KE [(Mx+M,)/Mx]. 

Given by KE, [Mx/(Mx + M,)] — KE. [Mx +M,)/Mx}. 
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The equation for a reaction between a target nucleus 
and a neutron with essentially zero kinetic energy is 
given by 


A 1 A+1 A+1—b 
z4* + on > [z\* aaa 18 


+2, Y5+Q, (2.96) 


where 


z,X4 = target nucleus with atomic num- 
ber Z, and mass number A, 


on’ =incident neutron with atomic 
number 0, mass number 1, and 
essentially zero kinetic energy, 


[z,X4*3] =compound nucleus in its excited 
state, 


(zZ;-Z,)¥4*1~° = product nucleus, 
z,Y° = emitted particle, 
Q = nuclear reaction energy. 


Equation 2.96 requires that a balance exist with respect 
to both the atomic number and the mass number and is 
useful in analyzing nuclear reactions. 

The nuclear reaction energy Q is determined by the 
change in mass experienced by the individual com- 
ponents of the reaction. The original equation can be 
written as a total energy equation. (All masses are 
expressed in terms of their energy equivalent, with 1 
amu being equal to 931 MeV.) The explicit expression 
for Q is given by 


Q=931 [My +M, —(My+M,)] , (2.97) 


where 


Q = nuclear reaction energy (MeV), 
My = mass of the target nucleus (amu), 
M,, = mass of the incident neutron (amu), 
My = mass of the product nucleus (amu), 


M, = mass of the emitted particle (amu). 


It is noted that if Q is positive the reaction is exoergic 
and neutrons with essentially zero kinetic energy can 
produce the reaction. If, on the other hand, Q is 
negative, the reaction is endoergic and the neutron must 
possess sufficient kinetic energy to make up the deficit. 
It can be shown that the threshold energy required of 
the neutron is given by 


My +M, 


My (2.98) 


Er 
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This follows from the law of the conservation of 
momentum when applied to the formation of the 
compound nucleus, since it requires that some of the 
neutron’s kinetic energy must appear as translational 
energy of the compound nucleus and is therefore 
unavailable as excitation energy. 

As pointed out above, those processes that fall in the 
category of nuclear reactions all contribute to the 
absorption cross section (see Fig. 2.10) and include 
radiative capture, (n,2n) reactions, charged-particle re- 
actions, and fission. Of these processes, radiative cap- 
ture is the most important experienced by low-energy 
neutrons in shields. Following radiative capture, the 
complete de-excitation of the compound nucleus is 
accomplished by the emission of gamma radiation. The 
binding energy of the neutron in the compound nucleus 
resulting from slow-neutron capture is given by 

E, =931 (My +M,)—-My] . (2.99) 
In this case the total excitation energy is equivalent to 
the binding energy and is about 8 MeV. 

The other nuclear reactions become increasingly more 
important as the energy of the neutron increases, 
because other forms of de-excitation become more 
probable. However, epithermal radiative captures do 
occur and the gamma rays thus produced may con- 
tribute significantly to or even dominate the overall 
shielding problem. 

Because of limitations in basic nuclear theories and 
the scarcity of experimental data, the standard practice 
has been to assume that the energy spectra of the 
epithermal radiative capture gamma rays are the same as 
those due to thermal capture. Yost and Solomito’? 
have shown that this assumption is inappropriate for 
many materials and that the exclusive use of thermal- 
capture spectra can introduce considerable error in the 
subsequent gamma-ray transport problem and may, in 
fact, be the most questionable of the nuclear data used 
in the calculation, Yost’? developed a method for 
determining capture gamma-ray spectra as a function of 
neutron energy (capture states) and wrote a machine 
code to generate capture gamma-ray spectra with a 
desired group structure resulting from capture of 
neutrons of each energy group. The code is based upon 
an appropriate nuclear model involving the spin and 
parity of the nucleus. 

The (7,2n) reaction involves the emission of two 
neutrons and its occurrence is usually limited to 
neutrons with high incident kinetic energies — greater 
than 10 MeV. The product nucleus is an isotope of the 
target nucleus and is frequently radioactive. The thresh- 
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old energy for the (#,2m) reaction is approximately 
equal to the binding energy per nucleon of the 
compound nucleus and is related to the threshold 
energy for the photoneutron (y,n) reaction by the 
relation 


M,+tM, 


Exam = Ereyn) (2.100) 


For certain very light or very heavy nuclides the 
thresholds for the (n,2n) reaction are sometimes much 
less than 10 MeV. For example, 7H, §Li, and 3B have 
thresholds of 3.34 MeV, 6.2 MeV, and 1.85 MeV 
respectively, and 723Bi, 732Th, and 738U have thresh- 
olds of 7.4 MeV, 6.44 MeV, and 6.0 MeV respectively. 
The (n,2n) production of neutrons is usually ignored in 
shield design; however, it is possible that this effect 
could be important for shield systems containing 
significant amounts of certain materials (beryllium, for 
example). 

In charged-particle reactions the de-excitation of a 
compound nucleus formed by neutron capture is 
accomplished by the emission of a charged particle, 
either a proton or an alpha particle.* This type of 
reaction is most probable in the MeV range; however, 
there are several important charged-particle reactions 
for slow neutrons. Consider the following two methods 
for producing tritium: 


3 He + in > [4He] * > 3H + 1H + 0.764 MeVK. E. 
67:41 717)" 53 4 
3 Li + on > [3 Li] > 1H+ 4 He + 4,785 MeV K.E. 


Another important charged-particle reaction involves 
the }°B isotope of boron: 


*An alpha particle is He’. 


OR + dn > [LB] * > TLi* + He + 2.30 MeV K.E.. 
ZLi* > JLi+ 0.478 MeV y 


The !°B(n,q) reaction is useful for the detection of 
slow neutrons because of its large thermal cross section. 
Also, }°B serves as an important constituent for shields 
designed to suppress high-energy capture gamma-ray 
production. 

In fission reactions the compound nucleus separates 
into fission fragments, concomitantly releasing neutrons 
and gamma rays. This process rarely occurs except ina 
few of the very heavy elements, which, with the 
exception of depleted uranium, are not used in shields. 
(See further discussion of fission reaction in Section 
2.6.) 


2.5.3, NEUTRON TRANSPORT CROSS SECTION 


The neutron transport cross section is defined as 
O41, =% - S, 6,(8) cosd dQ , (2.101) 
nT 


where 
6, = total microscopic cross section, 


0,(6) = differential microscopic cross section for elas- 
tic scattering through a given angle @ per unit 
solid angle, 


6 = scattering angle in the laboratory frame of 
reference. 


The reaction rate 2,, ® may be regarded as an estimate 
of the rate at which neutrons lose their forward 
momentum. In this context the transport cross section 
has been used as an approximation to the energy- 
dependent removal cross section (see Chapter 3, Section 
3.8.2). Quite often, o,, appears naturally in the results 
of a derivation’? which has other formal bases. 


2.6. Neutron Production Processes 


Except for the small contributions made by such 
reactions as the (n,2n), (n,3n), and (y,) reactions, the 
neutrons emitted in the detonation of a nuclear weapon 
are produced directly by the fission process or the 
fusion process. In the fission process the resulting 
neutrons cover a wide range of energies, whereas in the 
fusion process the emissions are essentially mono- 
energetic, with the energies being determined by the 
type of fusion reaction. 


2.6.1. FISSION REACTION 


The details of the fission process are very complex 
and only a brief description is presented here.* As in 
other nonelastic neutron interactions (see Section 
2.5.2), fission begins with the formation of a compound 
nucleus in an excited state. If the excitation energy 
exceeds a critical energy, then fission may occur.t 
Calculated values of the fission threshold energy are 
shown in Table 2.2 for a wide range of atomic 
weights.1* It is apparent that only the very heavy 
nuclides (A > 230) have reasonably low threshold 
energies. 

Experimental values of the fission thresholds for the 
nuclear fuel nuclides of practical interest are given in 
Table 2.3 (ref. 14). The significance of a negative 
threshold energy (73°U, ?35U, and ?%°Pu) is that 
neutrons with essentially zero kinetic energy can 
produce the fissioning of these nuclides (thermal 
fission). All the other nuclides (73?Th, 734U, 73°U, 
2381) and 237Np) can undergo fast fission only. 

When fission occurs, the highly unstable compound 
nucleus almost always splits into two fission fragments. 
Splitting into three or more fragments is theoretically 


* Although this discussion is limited to the fission of the heavy 
nuclei by neutrons, it is pointed out that any nucleus except 
hydrogen may undergo fission if its level of excitation is 
sufficiently high and if this excitation can be imparted to the 
target nucleus by its interaction with gamma radiation (photo- 
fission) or with a variety of high-energy particles. _ 

tThe alternative to the fission mode of de-excitation is the 


emission of gamma radiation, which returns the compound 
nucleus to its ground state, i.e., effectively radiative capture. 
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possible, but such an occurrence is rare and is generally 
ignored. The nucleons are more firmly bound within 
their respective fission fragments than they were in the 
original fissile nuclide, thereby liberating the very large 
energy of fission (about 200 MeV per fission for the 
very heavy nuclides). Fissioning usually results in the 
prompt emission of several neutrons, the average 
number released per fission being a function of the 
neutron energy and the fissile material. For example, 
fissioning of 735U and 7%°Pu by thermal neutrons 
produces about 2.5 and 2.9 neutrons/fission respec- 


Table 2.2, Computed Neutron Fission Thresholds 
as a Function of Nuclear Mass* 


Mass No., Fission Threshold 

A (MeV) 
16 18.5 
60 48 

100 47 

140 62 

200 40 

236 ~s 


*From ref. 14, p. 107. 


Table 2.3, Experimental Neutron Fission Thresholds 


of Nuclear Fuel Nuctides* 

Target Compound Fission Threshold 
Nucleus Nucleus (MeV) 
aah 233TH 1.3 
2334 2344 <0 
2344 2351 0.4 
2354 2364 <0 
2361 23714 0.8 
2381, 239y 1.2 
237N»p 238N, 0.4 
239p) 240p, <0 


*From ref, 14, p. 108. 
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tively. In addition, some of the fission fragments 
contain nuclides with neutrons in excess of those 
required for nuclear stability. These unstable nuclides 
have half-lives up to about 1 min, and the neutrons 
emitted are called delayed neutrons.* (It is these 
neutrons that make control of a nuclear reactor 
relatively easy. If all neutrons were prompt, every 
reactivity excursion would be prompt critical, and a 
reactor would be difficult to operate safely.) Also, 
because the fission fragments are neutron-rich, they 
undergo radioactive decay, emitting primarily beta and 
gamma radiations. 

The mass numbers of the fission fragments versus the 
corresponding fission yields are shown in Fig. 2.13 for 
the thermal fission of 735U, 73®U, and 2°? Pu (ref. 15). 
The fission yield is defined as the percentage of the 
total fissions yielding products of a given mass number. 
As the energy of the incident neutron increases, the 
probability of producing two fission fragments of 
nearly the same size increases and the characteristic dip 
of the thermal-fission-yield curves at a mass number of 
about 120 disappears. 

The energy spectrum of fission neutrons has been 
experimentally determined for 79*U, the nuclide used 
most frequently as a nuclear fuel, and is most con- 
veniently represented by formulas known as the fission 
spectra of Watt,’® of Cranberg et al.,'7 and of Los 
Alamos.!® . 

The formula for the fission spectrum as determined 
by Watt is 


N(E) = 0.484 sinh /2E e-£ , (2.102) 


where M(E) is the fraction of neutrons emitted per 
fission with an energy F, in MeV, per unit energy range. 
This formula is regarded as adequate for the energy 
range 0.075 to 17 MeV, which includes the energy 
range of importance to shielding (3 to 17 MeV). Table 
2.4 presents values of N(E), together with two other 
useful quantities, F(Z) and 1 — FE); F(E) is the 
fraction of neutrons with energies greater than the 
energy £ and is given by 


F(E) = Sf *° NE") dE’ . (2.103) 


The quantity 1 — F(E) is the fraction of neutrons with 
energies less than the energy £ and is given by 


*For thermal fissioning of 235) and 239 py the delayed 
neutron fraction is about 0.7% and 0.2% respectively. 


1-Fe)= f ie ME") dE’ . (2.104) 


The formula for the fission spectrum obtained by 
Cranberg et al.1” is 


ME) = 0.453 e-#/9-965 sinh /2.29E . (2.105) 
The measurements on which this formula is based 


involved an energy range of 0.18 to 12.0 MeV and 
represent no substantial deviation from Watt’s earlier 


Table 2.4, Watt Prompt-Neutron Spectrum Due to 
Fission of 735U by Thermal Neutrons@ 


E M(E) 

(MeV) (MeV—!) bo) FE) 
0.1 0.2023 0.014 0.9860 
0.2 0.2676 0.377 0.9623 
0.3 0.3068 0.663 0.9337 
0.4 0.3300 0.987 0.9013 
0.5 0.3450 0.1321 0.8679 
0.6 0.3528 0.1669 0.8331 
0.7 0.3557 0.2027 0.7973 
0.8 0.3542 0.2385 0.7615 
0.9 0.3513 0.2729 0.7271 
1.0 0.3446 0.3082 0.6918 
1.1 0.3368 0.3425 0.6575 
1.2 0.3271 0.3759 0.6241 
1.3 0.3179 0.4074 0.5926 
1.4 0.3073 0.4383 0.5617 
1.6 0.2841 0.4983 0.5017 
1.8 0.2604 0.5631 0.4369 
2.0 0.2371 0.6028 0.3972 
2.5 0.1839 0.7073 0.2927 
3.0 0.1384 0.7876 0.2124 
3.5 0.1026 0.8476 0.1524 
4.0 0.07472 0.8914 0.1086 
4.5 0.05381 0.9192 0.0808 
5.0 0.03847 0.9463 0.0537 
5.5 0.02729 0.9624 0.0376 
6.0 0.01916 0.9739 0.0261 
6.5 0.01336 0.9821 0.0179 
7.0 9.316(—3)? 0.9876 0.0124 
7.5 6.436(—3) 0.9916 0.0084 
8.0 4.423(-3) 0.9943 0.0057 
8.5 3.034(—3) 0.9960 0.0040 
9.0 2.076(—3) 0.9973 0.0027 
9.5 1.418(—3) 0.9982 0.0018 

10.0 9,53(—4) 0.9988 1.211(—3) 
11.0 4,384(—4) 0.9994 5.501(—4) 
12.0 1,994(-4) 0.9998 2.473(—4) 
13.0 8.953(—-5) 0.9999 1.103(—4) 
14.0 3.997(—5) 1.0 4.875(-5) 
15.0 1.771(—5) 1.0 2.148(—5) 


2From ref, 16, 
bRead: 9.316 X 10°, etc. 
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Table 2.5. Cranberg et ai, Prompt-Neutron 
Spectrum Due to Fission of 235U by 
Thermal Neutrons® 


E N() G(E) 
(MeV) (MeV"!) HME) (MeV) 
0.00 0.000 1.000 1.981 
0.25 0.290 0.948 1.973 
0.50 0.347 0.867 1.942 
0.75 0.358 0.779 1.887 
1.00 0.347 0.690 1.810 
1.25 0.325 0.606 1.715 
1.50 0.298 0.528 1.608 
1.75 0.268 0.457 1,493 
2.00 0.239 0.394 1.374 
2.50 0.184 0.2884 1.138 
3.00 0.138 0.2082 0.919 
3.50 0.102 0.1486 0.726 
4.00 0.0738 0.1050 0.563 
4.50 0.0528 0.0736 0.430 
5.00 0.0375 0.0512 0.324 
5.50 0.0263 0.03544 0.241 
6.00 0.0184 0.02439 0.178 
6.50 0.0127 0.01669 0.130 
7.00 8.77(-3)? 0.01138 9.41(—2) 
7.50 6.01(—3) 7.72(-3) 6.76(—2) 
8.00 4.10(-3) 5,22(—3) 4.83(—2) 
8.50 2.79(—3) 3.52(—3) 3.43(—2) 
9.00 1.88(-—3) 2.364(—3) 2.42(—2) 
9.50 1,27(-3) 1.583(—3) 1.70(-2) 
10.00 8.56(—4) 1.058(—3) 1.187(—2) 
10.50 5.74(—4) 7.05(—4) 8.25(—3) 
11.00 3.84(—4) 4.686(—4) 5.72(-3) 
11.50 2.56(—4) 3.108(—4) 3.947(—3) 
12.00 1.70(-4) 2,058(—4) 2.714(-3) 
13.00 7.48(—5) 8.97(—5) 1.271(-3) 
14.00 3.26(—S) 3.88(—S) 5.88(—4) 
15,00 1.41(-5) 1.67(—5) 2.69(—4) 
16.00 6.07(-6) 7,10(—6) 1.22(-4) 
17.00 2.59(-6) 3.00(—6) 5.48(-5) 
18.00 1.10(-6) 1.26(—6) 2.44(-5) 


“From ref. 18, p. 48. 
bRead: 8.77 X 107°, ete. 


results, Values of M(Z) and FY£) as calculated from 
Cranberg’s spectrum are presented in Table 2.5, along 
with values of G(E), which is the energy per fission 
carried by neutrons having energies greater than F and 
is given by 


G(E) = J. * E’ NE’) dE’ . (2.106) 
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The formula developed by Los Alamos!® for the 
fission spectrum in the energy range below 9 MeV is 


N(E) = 0.77E3/2 e—0-776E | (2.107) 


This equation is a simpler approximation to the 
experimental data and is used primarily for reactor core 
calculations. The error associated with the formula is 
less than 12% for the specified energy range. For the 
energy range 4 to 14 MeV the Los Alamos formula has 
been reduced to the simple exponential 


N(E) = 1.75 e—9-766E (2.108) 


The error associated with this formula is always less 
than 15% and is less than 7% for energies from 5 to 13 
MeV, 


2.6.2, FUSION REACTION 


The fusion reaction involves the nuclear combination 
of two light nuclei having sufficient kinetic energy to 
overcome the repulsive force that exists between nuclei. 
Because the mean binding energy per nucleon for 
certain light nuclei is less than that for nuclei of 
intermediate mass (essentially the same change-in-mass 
phenomenon which characterizes the fission process), 
the fusing process of these nuclei is accompanied by 
considerable energy release and by the emission of 
nucleons. 

The fusion processes of practical interest involve the 
light nuclei deuterium, tritium, and *He. Under proper 
conditions several combinations of these nuclei can 
undergo fusion, the reactions being represented in 
equation form as follows: 


2D +?D > 3He(0.82 MeV) + §n(2.45 MeV) 
?D+?D—>3T(1.01 MeV) + }H(3.02 MeV) 
2D + 7T > fHe(3.5 MeV) + 6n(14.1 MeV) 


21 + 3He > 3 He(3.6 MeV) + }H(14.7 MeV) 


It is the 2.45- and 14.1-MeV neutrons produced in the 
D-D and D-T reactions utilized in nuclear detonations 
that are primarily of concern in shield design. 

Unlike in the fission process, no prompt gamma rays 
are emitted in the fusion process and neither are the 
products radioactive. Bremsstrahlung radiation (see 
Section 2.4.2) is produced in conjunction with charged- 
particle deflections by the Coulomb field of other 
charged particles, but such radiation does not con- 
tribute significantly to the shielding problem. 


2:7. 


There are many possible processes by which gamma 
rays may interact with matter; however, the attenuation 
of gamma rays in a shield is completely dominated by 
three processes: the photoelectric effect, pair produc- 
tion, and Compton scattering. The three processes are 
described below, as are a number of processes that are 
usually neglected in shielding calculations. 


2.7.1, PHOTOELECTRIC EFFECT 


The photoelectric effect is characterized by the 
complete absorption of the gamma ray. The gamma ray 
interacts directly with an orbital electron, which is 
initially bound in an atom. The momentum is conserved 
by the recoil of the entire atom, thereby allowing a 


complete loss of energy by the gamma ray. Conse- - 


quently, the more tightly bound electrons have the 
greatest probability of absorbing the gamma ray. It has 
been shown? that about 80% of the photoelectric 
absorptions involve the K-shell electrons. The orbital 
electron (photoelectron) is ejected from the atom 
having a kinetic energy T = Ey — B,, where B, is the 
binding energy of the ejected electron. The remainder 
of the energy, which is equal to B,, appears as very soft 
characteristic X rays, usually referred to as fluorescence 
radiation. 

The photoelectric process does not lend itself to 
theoretical analysis and most cross-section data for this 
interaction are empirical, with theoretical methods 
being employed for interpolation and extrapolation of 
the data. The qualitative dependence of the probability 
for the photoelectric effect is given by the approxi- 
mation 


Zz 
SpE ~ > (2.109) 
7 
where E, is the incident photon energy and Z is the 
atomic number. Equation 2.109 shows that the photo- 
electric effect is important for low-energy gamma rays 
and large-atomic-number materials, such as lead. Table 
2.6 gives, for the elements of greatest importance to 
shielding, the energy at which the photoelectric effect 
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Interactions of Gamma Rays 


Table 2.6. Gamma-Ray Energies at Which 
Photoelectric Effect Provides One-Half 
the Total Absorption Coefficient* 


Gamma-Ray AS net 
Atomic No., Ionization 
Element Energy 

Z (MeV) Energy 

(MeV) 
1 Hydrogen 10+ 1.4x 1075 
4 Beryllium 0.011 2.2 x 10% 
6 Carbon 0.016 2.8 x 10% 
8 Oxygen 0.025 5.2 x 10% 
13 Aluminum 0.046 1.5 X 1073 
20 Calcium 0.079 4.0 x 107? 
26 Iron 0.11 6.9 x 1079 
42 Molybdenum 0.195 2.0 x 10°? 
50 Tin 0.25 2.9 x 107? 
74 Tungsten 0.42 6.06 X 10°? 
82 Lead 0.50 8.8 X 10-2 
92 Uranium 0.62 1.16 X 107? 


*From ref, 18, p. 144. 


provides one-half the total absorption coefficient, along 
with the ionization energy of the K shell. 

The fluorescence radiation associated with the photo- 
electric effect is significant only when the interaction 
occurs in heavy elements. Even then the range of the 
emitted X rays is comparable to that of the photo- 
electrons, and their contribution to the penetrating 
dose is so small that they can be assumed to be 
absorbed at their point of origin. 


2.7.2, PAIR PRODUCTION 


The pair-production process has a threshold energy of 
1.02 MeV and becomes increasingly more probable as 
the gamma-ray energy increases above 1.02 MeV. In this 
process the gamma ray is completely absorbed and an 
electron-positron pair appears in its place. Pair produc- 
tion must occur within the electrostatic field of a 
charged particle. The momentum is conserved in the 
process by the recoil of the charged particle, thereby 
making possible the complete absorption of the gamma 
ray. Either a nucleus or an atomic electron can provide 
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the necessary electrostatic field; however, due to the 
greater charge of a nucleus, the interactions are pri- 
marily nuclear. 

The sum of the kinetic energies of the electron pair is 
equal to the initial gamma-ray energy less the energy 
equivalent of their rest masses and the kinetic energy of 
the recoil nucleus: 


(2.110) 


T,+T_ =F —2mgc? — Teacoil 
The kinetic energy acquired by the recoil nucleus can 
be neglected because of the relatively large mass of the 
nucleus. If, however, the process occurs in the field of 
an atomic electron, considerable energy is transferred to 
the recoil electron and a triplet is produced. The triplet 
consists of the newly created electron pair (one 
positron and one electron) and of the recoil electron. 
Because of momentum considerations, the threshold 
energy for triplet production is 2.04 MeV. 

Every positron created in pair production loses its 
kinetic energy by ionization and then ultimately inter- 
acts with an ordinary electron, converting the electron 
and itself into electromagnetic radiation called 
annihilation radiation. Thus for each positron, two 
0.51-MeV gamma rays are emitted in opposite direc- 
tions. These secondary gamma rays are most important 
for high-energy gamma-ray sources (greater than 6 
MeV) and medium-to-heavy elements. They are usually 
assumed to be absorbed at the point of formation; 
however, in the more refined Monte Carlo and S,, 
calculations, transport of the 0.51-MeV annihilation 
gamma rays can be included. 

Aside from the effects of screening,* the nuclear 
pair-production cross sections Gppy are theoretically 
and experimentally proportional to Z?. Analytical 
expressions are available only for two special cases: 


(1) If screening is neglected and moc? < E, < 137 
moc? Z~1/3, then 


(2.111) 


(2) If complete screening is considered and £,, > 137 
Moc? Z~ 1/3 : 


Gap road? en (183 Z~ 1/3) — z| , (2.112) 


*The positive charge of the nucleus can, in effect, be reduced 
by the cloud of negatively charged orbital electrons. 


where 


1 e? f 
9 =>" 
° 137 oc? 
=5.80 X 10°78 cm?/nucleus , (2.113) 
moc? =energy equivalent of the rest mass of an 
electron 
= 0.5110 MeV, 


Mo = electron rest mass 
= 0.5488 X 107° amu =9.108 X 10°78 g, 


c = speed of light in a vacuum 
= 2.99793 X 101° cm/sec. 


Table 2.7 gives the gamma-ray energies at which pair 
production provides one-half the total absorption coef- 
ficient for elements of greatest interest in shielding 
calculations. 

Electronic pair production can occur when the 
gamma-ray energy exceeds 2.04 MeV. The analysis of 
this effect is difficult because of screening by other 
electrons and by the nucleus. However, the ratio of the 
total cross section for Z electrons of a given atom, 
ZOpp x, to the cross section for nuclear pair production, 
Op pn; can be written as 


ZOppE ee 


— 2.114 
Sean C2" (2.114) 


Table 2.7, Gamma-Ray Energies at Which Pair 
Production Provides One-Half the Total 
Absorption Coefficient* 


. Gamma-Ray 
nase NG Element Energy 
(MeV) 
1 Hydrogen 78 
4 Beryllium 35 
6 Carbon 28 
8 Oxygen 20 
13 Aluminum 18 
20 Calcium 12 
26 Iron 9.5 
42 Molybdenum _ 75 
50 Tin 6.5 
74 Tungsten $.2 
82 Lead 5.0 
92 Uranium 4.8 


*From ref, 18, p. 145. 


where the constant C depends upon E,, and varies from 
about 2.6 at EB, = 6.5 MeV to about 1.2 at £ = 100 
MeV. 
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The effect of the electronic pair production can be 


included in a total pair-production cross section as 


1 
Opp = Oppn t+ ZOppr = Oppn (1 a 4) » (2.115) 


where 1/CZ =0 if B, <2.04 MeV. 


‘2.7.3, COMPTON SCATTERING 


The Compton scattering process, named for the 
‘experimentalist A. H. Compton, can be described as an 
elastic collision between the incident gamma ray and a 
free electron.* ‘Under these conditions the kinetic 
energy, as well as the momentum, is conserved, which 
results in a precise relation between the loss of 
gamma-ray energy and its change in direction. The 
functional representation of this relation is given by 


E 
fx3—__—-——__, 


(2.116) 
Eo 
1+ (1 — cos@) 


Moc? 
where 
Eo = initial gamma-ray energy, 
F = final gamma-ray energy, 
§=angular change experienced by the gamma 
ray, 
moc? = 0.5110 MeV, 
Mo = 0.5488 X 1073 amu = 9.108 X 10728 g, 
c¢ = 2.99793 X 10'° cm/sec. 


A plot of Eq. 2.116, which is referred to as Compton’s 
equation, is shown in Fig. 2.14. 

Another useful form of Compton’s equation is ob- 
tained if the gamma-ray energy & and wavelength X are 
expressed in units of the electron rest-mass energy (mc? 
= 0.5110 MeV) and the Compton wavelength (h/mc = 
0.02426 A) respectively. In these units the wavelength 
and energy variables are simply reciprocals of each 
other and Compton’s equation assumes the form 


A —Ao = 1 —cosé , (2.117) 


where ) is the final wavelength. 


*As the energy of the incident gamma ray increases, even the 
tightly bound orbital electron behave as free electrons. 
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Fig. 2.14. Reduction in Gamma-Ray Energy by Compton 
Scattering. (Eo = initial energy; EF = reduced energy.) 


The differential cross section o(@) for Compton 
scattering is defined as the cross section per electron for 
scattering through a given angle @ per unit solid angle. 
Expressed in Thomson units, this cross section is given 
by the Klein-Nishina formulat?°® 


BD SRE Po eh 2s 
E Eo 


= re) , (2.118) 
dQ. 16m £2 


TThis form of the equation does not include polarization 
effects that are due to scattering. However, this effect is of no 
consequence to practical problems in radiation transport. 


Cc 


c 


where 1 Thomson unit = (87/3) (e?/mce?)*? = 0.665 
barn. The dependence of the differential cross section 
on the scattering angle 6 for various initial energies is 
illustrated in Fig. 2.15. It is noted that o(@) is 
symmetric only for the nonrelativistic limit E > 0. As 
the initial energy increases, the cross section becomes 
more and more sharply peaked in the forward direction. 
Consider an initial energy E = 2.5 MeV: the cross 
section decreases to one-half its initial value at a 
scattering angle of about 20° and then drops to about 
3% of its initial value for @ = 180°. This characteris- 
tically forward scattering of the high-energy gamma 
rays strongly influences the nature of the deep- 
penetration problem. 

The total cross section for Compton scattering per 
electron is obtained by integrating the differential cross 
section over 47 solid angle: 


o¢(E) = tos o(0) dQ, (2.119) 
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Fig. 2.15. Differential Klein-Nishina Cross Sections. 
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where dQ. = 27 sin@ d@. Substituting Eq. 2.118 into Eq. 
2.119 and integrating provides the following expression 
for the total cross section: 


o(E) = Hee poe +2) —In(1+ 2»)| 


1 1+3E 
+2 jae ony St 0190 
2g nt rane Ca 


where the units of the gamma-ray energy £ are given in 
terms of the electron rest-mass energy (moc? = 0.5110 
MeV). For small gamma-ray energies (E < 1 moc?) the 
value of o.(£) approaches one Thomson unit and for 
large gamma-ray energies it approaches (*4£) (In 2E + 
*h). 

The cross section related to the average fractional 
energy carried off by the scattered gamma ray is 
defined as 


o¢ (6) =f" a o(6) da, (2.121) 


which on integration yields 


_3 |In(1+2£) , 2(1 + £) (QE? — 2F — 1) 
0c A)=3 EB + F*(1 + 2E)? 


8° | . (2.122) 


+ ——— 
3(1 + 2£)° 


The other cross section of importance relates to the 
average fractional energy lost by a gamma ray suffering 
a Compton collision; the fraction increases mono- 
tonically as the gamma-ray energy increases. This cross 
section represents the energy absorption or energy 
deposition* cross section and is defined by 


oc =f" = aaa, 


(2.123) 


*The term energy deposition was suggested by A. B. Chilton 
of. the University of Illinois as a replacement for the term 
energy absorption since the ICRU has designated the energy 
absorption coefficient to mean an essentially similar quantity 
but with correction terms to account for subsequent transport 
of the energy after having made an absorption collision. A 
discussion of these considerations is given in ORNL-RSIC-16 
(ref. 21). One must be careful in the use of such correction 
terms, which should be consistent with the transport model 
assumed or energy may not be conserved. 


which is also simply expressed as 


oc) = 0-(E) - oo (4) : (2.124) 

It is noted that the various Compton cross sections 
are usually given on a per electron basis since the 
scattering occurs with free electrons. The cross sections 
are simply multiplied by the atomic number Z to obtain 
the atomic cross section. More detailed descriptions of 
the Compton cross sections are presented else- 
where.! 819 


2.7.4. NEGLECTED PROCESSES 


As pointed out above, the three processes described in 
the preceding paragraphs are usually the only gamma- 
ray interactions considered in shielding attenuation 
calculations since all other processes make only negli- 
gible contributions. Some of the interactions that are 
ignored are described below. 


Rayleigh Scattering 


Of the neglected processes Rayleigh scattering is the 
most important and of some slight significance at low 
gamma-ray energies, but it is still less probable than 
Compton scattering by 1 or 2 orders of magnitude for 
energies greater than 0.1 MeV. Rayleigh scattering is 
elastic coherent scattering by tightly bound electrons, 
and the process is more probable with heavy elements. 
In the energy range of interest to shielding, Rayleigh 
scattering results in small changes in direction and in 
essentially no energy loss. The net result of coherent 
scattering is to make the angular distribution of 
low-energy scattered gamma rays more nearly isotropic. 
Tabulations of total attenuation coefficients for practi- 
cal design calculations generally do not include the 
Rayleigh scattering cross section, which neglect tends to 
produce a slight overestimate of the gamma-ray flux at 
field points removed from the source. 


Thomson Scattering 


Thomson scattering is quite similar to Rayleigh 
scattering, with the interaction involving the nucleus 
rather than atomic electrons. The effect is very small 
and is truly negligible. 


Nuclear Resonance Scattering 


Nuclear resonance scattering involves the excitation 
of a nuclear level by the complete absorption of an 
incident gamma ray followed by the re-emission of the 
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excitation energy. The probability for nuclear reso- 
nance scattering is extremely small and its effect would 
be to increase the isotropy of the gamma-ray flux. 


Photodisintegration 


Photodisintegration of nuclei is possible if the 
gamma-ray energy exceeds the separation energy of a 
nucleon. Except for the °Be(y,) and ?H(y,n) reac- 
tions, this effect is confined to gamma rays with 
energies in excess of 8 MeV. The cross sections for 
photodisintegration are usually negligible when com- 
pared with those for Compton scattering and pair 
production. 


2.7.5. GAMMA-RAY CROSS SECTIONS 


The total gamma-ray interaction cross section for the 
usual shielding calculation is given by 


O,= Opp t OpptZoc, (2.125) 


where 


o, = total gamma-ray interaction cross section per 
atom, 


Opg = photoelectric microscopic cross section per 
atom, 


Opp = pait-production microscopic cross section per 
atom, 


Oc = total Compton cross section per electron, 


Z = atomic number. 


More useful macroscopic variations of the total gamma- 
ray cross section are the linear attenuation coefficient 
H, given in units of cm™' and defined as 


PNo 


H=0,—-, (2.126) 


where p is the density (g/cm*), No is Avogadro’s 
number = 6.023 X 107? nuclei/g-atom, and the mass 
attenuation coefficient is designated by u/p in units of 
cm?/g. 

The gamma-ray energy deposition coefficient is 
usually taken as 


PNo 


Ho 


(Opp + Opp +Zoc ), (2.127) 
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where gc, is the Compton energy-deposition cross 
section per electron. This relation presumes that the 
photoelectric and pair-production reactions are pure 
absorptions; other variations are sometimes used (see 
Section 2.9.3). These cross sections are often reported 
as mass energy absorption coefficients, u,/p. Extensive 
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tables for the cross sections of these gamma-ray 
interactions, as well as others, have been prepared by 
Plechaty and Terrall?? for elements 1 through 100 and 
gamma-ray energies ranging from 0.001 to 100 MeV, 
and later data for the energies 1 keV to 1 MeV have 
been published by McMaster et al.?* 


2.8. Gamma-Ray Production Processes 


It was pointed out in Section 2.4.2 that the electro- 
magnetic radiations of primary concern in shielding are 
the nuclear transition gamma rays. Gamma rays pro- 
duced by a weapons burst are due to a number of 
sources, with only one, the fission reaction, a direct 
result of the detonation. All other sources are produced 
either by the interaction of fission or fusion neutrons 
with various materials or by the radioactive decay of 
the fission products. The important sources are listed 
below. 


2.8.1. FISSION REACTION 


The fission reaction utilized in weapons detonations is 
always accompanied by the emission of prompt gamma 
rays, whose energies total about 8 MeV per fission. 
These prompt gamma rays, like prompt neutrons, are 
released whenever the compound nucleus formed by an 
incident neutron and a target nucleus splits into two (or 
more) fission fragments (see Section 2.6.1). The energy 
spectrum of such gamma rays is known reasonably well 
only for 775 U, measurements having been made for this 
isotope by Motz?* and by Maienschein et al.?* 


2.8.2. NEUTRON INELASTIC SCATTERING 


In the neutron inelastic-scattering process described in 
Section 2.5.2 the energy of the target nucleus is 
elevated to an excited state, from which it returns to 
the ground state by the emission of gamma rays. In 
weapons detonations the target nuclei can be weapons 
materials or environmental materials, such as those 
comprising the atmosphere and the ground. Investiga- 
tions of the production of such gamma rays have been 
carried out by Dickens and Perey?® and by Orphan and 
Hoot.?7 (See also Chapter 6 of this Handbook.) 


2.8.3. NEUTRON CAPTURE 


-The neutron radiative-capture process, also described 
in Section 2.5.2, results in the de-excitation of the 
compound nucleus via the emission of gamma rays. And 
again, the target nuclei can be weapons materials or 
environmental materials such as those comprising the 
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atmosphere and the ground. The gamma rays emitted 
following radiative capture can be given off immedi- 
ately, in which case they are called capture gamma rays, 
or they can be emitted at later times, in which case they 
are called activation gamma rays. Activated materials 
can contribute to fallout or can prevent access to 
regions in the vicinity of a detonation. 

The radiative capture process is most important for 
low-energy neutrons. Early data on capture gamma rays 
produced by neutron captures in various nuclei were 
published by Bartholomew and Campion,?® Keller et 
al.,?° and Troubetzkoy and Goldstein,?® and later 
information has been presented by Maerker and 
Muckenthaler,?! Bartholomew et al.,?? and Groshev et 
al.33>3* Data on long-term emitters of activation 
gamma rays have been published by Crocker and 
Conners? ® and by Crocker and Wong.?° 


2.8.4. FISSION-PRODUCT DECAY 


The fission fragments produced in the fission process 
undergo radioactive decay through the emission of 
gamma (and beta) radiations. These gamma rays consti- 
tute an important source from the moment of detona- 
tion to long periods following the detonation. Studies 
of the early-time (< 1 min) fission-product gamma rays 
have been made by Engle and Fisher.?’ Data on the 
total radiative decay of fission products for times up to 
70 years have been reported by Crocker and Tumer.?® 
The total gamma-ray energy released through fission- 
product decay is about 7 MeV per fission. 


2.8.5. BREMSSTRAHLUNG 


Electromagnetic radiation, called bremsstrahlung, is 
emitted when a free charged particle experiences an 
acceleration or deceleration by the electrostatic field of 
a nucleus. The total bremsstrahlung produced varies as 
the square of the atomic number of the recoil nucleus 
and inversely with the square of the mass of the 
incident particle. Bremsstrahlung is characterized by a 
continuous spectrum, with the energy emitted by a 
particular particle ranging from zero up to its total 
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kinetic energy. Because of the strong dependence on 
the incident particle’s mass, charged particles other than 
electrons can be completely neglected insofar as the 
production of bremsstrahlung is concerned. 

Large numbers of electrons are produced as second- 
aries during the attenuation of gamma rays; for ex- 
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ample, photoelectrons, Compton recoils, and the 
electron-positron pairs. All these electrons will produce 
numerous bremsstrahlung radiations; however, their 
energies will be relatively low and for all present 
shielding calculations they are assumed to be absorbed 
at the point of formation. 


2.9. Quantities Used to Describe Responses to Radiation 


Although calculations of the interactions of radiation 
with matter deal with the fluxes and currents described 
in Section 2.1, such quantities are not easily measured 
in the laboratory, especially as a function of energy and 
direction. However, the property of a radiation field 
called dose, which is related to the energy deposited in 
a medium by the radiation interacting within it, is 
determinable, and this is the quantity that is. used 
to relate a radiation field to the biological damage it 
causes in a human body. 

In an ideal situation a given dose would always 
produce the same biological effect irrespective of the 
nature and energy of the radiation or of the body organ 
being irradiated. But studies by radiobiologists in which 
doses and biological damage to specific organs are 
correlated for known radiation fields have established 
that nature is not that simple. In addition to being a 
function of absorbed energy, biological responses are 
functions of the irradiated organ and of the type, rate, 
and energy of the radiation. From this discovery has 
evolved the concept of a relative biological effect 
(RBE), which is a weighting factor that can be used to 
compare the biological effects produced by the same 
physical dose (same amount of energy deposited) due 
to radiation of a different type and/or energy. When the 
physical dose is multiplied by the RBE it becomes a 
biological dose. 

The foregoing gives some indication of the problems 
associated with establishing doses for shield design that 
are related to biological hazards. The shield designer’s 
task is further complicated by the perturbation of the 
radiation field by the human body, an effect not 
usually included in a typical shielding calculation. 
Obviously the dose received at a particular location 
within the body is not the same as the dose in a small 
detector at the same location in space with the body 
absent. It is apparent that shielding studies should 
include anthropomorphic phantoms as part of the 
shield configuration, a theoretically possible but usually 
impractical consideration. 

As an alternative, slab or cylindrical phantoms with 
compositions resembling that of the human body have 
been used, and the doses have been calculated or 
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measured as a function of depth in the phantom for a 
given incident radiation field. The results are then used 
as response functions to relate an unperturbed radiation 
field (usually called a free field) to the dose in a human 
body if it had been present. 

It will be apparent in the following discussion on the 
various quantities used to define physical and biological 
doses that confusion has developed in the terminology 
and in the definition of units, both because of the 
burgeoning nuclear science and technology and because 
radiobiologists and shield designers basically differ in 
their viewpoints. The International Commission on 
Radiation Units and Measurements* recognized this 
confusion and in an effort to mitigate it recom- 
mended’?! a consistent set of definitions and units. 
While shield designers have accepted most of their 
recommendations, they have continued to use terms 
not included in the ICRU recommendations because 
these terms are so ingrained in the shielding field and 
because the ICRU did not include all the concepts 
needed in shield design. In the descriptions given below 
the ICRU recommendations and the traditional view- 
points are contrasted. 

In addition to descriptions of the quantities them- 
selves, this section includes data for some of the 
quantities considered most pertinent for shield design. 


2.9.1. ABSORBED DOSE 


Absorbed dose is the energy imparted by radiation to 
a unit mass of matter and as such is a physical quantity 
as opposed to a biological effect. A formal definition 
for the energy imparted, £,,, is: The energy imparted 


*The International Commission on Radiation Units and 
Measurements has as its principal objective the development of 
internationally acceptable recommendations regarding quanti- 
ties and units of radiation and radioactivity, procedures suitable 
for the measurement and application of these quantities, and 
physical data needed in the application of these procedures. 
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by ionizing radiation* to the matter in a volume is the 
difference between the sum of the energies of all the 
directly and indirectly ionizing particles which have 
entered the volume and the sum of the energies of all 
those which have left it, plus the energy equivalent (Q) 
of any decrease in rest mass that took place in nuclear 
or elementary particle reactions within the volume. For 
all practical purposes the energy imparted is usually 
equal to the heating effect, but in some cases part of 
the energy deposited may result in changes in inter- 
atomic bond energies. 

As proposed by the ICRU,! the energy imparted can 
be expressed as 


Ep = DE in - Least ye, 


where 


(2.128) 


» Ein = the sum of the energies (excluding rest 
energies) of all those directly and indirectly 
ionizing particles which have entered the 
volume, 


y Ex =the sum of the energies (excluding rest 
energies) of all those directly and indirectly 
ionizing particles which have left the volume, 


x Q=the sum of all the energies released, minus 
the sum of all the energies expended, in any 
nuclear reactions, transformations and ele- 
mentary particle processes which have 
occurred within the volume. 


The absorbed dose (D) is the quotient of AE, by 
AM, where AE! is the energy imparted by ionizing 
radiation to the matter in a volume element and AM is 
the mass of the matter in that volume element: 


AE 
Hea. 


AM (2.129) 


*lonizing radiation is a radiation consisting of directly or 
indirectly ionizing particles or a mixture of both, ‘Directly 
ionizing particles are charged particles (electrons, protons, alpha 
particles, etc.) having sufficient kinetic energy to produce 
ionization by collision. Indirectly ionizing particles are un 
charged particles (neutrons, gamma rays, etc.) which can 
liberate directly ionizing particles or can initiate a nuclear 
transformation. ; 

tThe notation AE 'p implies that the volume element AV 
associated with the element of mass AM be of an appropriate 
size such that the limiting process D = AE ‘pl 4M yield a 
meaningful estimate of the absorbed dose. The difficulties arise 
because the interactions between radiations and atoms are 
macroscopic and the limiting process cannot proceed to the 
usual mathematical limit AM > 0. 
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The special unit of absorbed dose is the rad, which is 


1 rad ea J/kg = 100 ergs/g . 


aa (2.130) 


The absorbed dose rate is the quotient of AD divided 
by At, where AD is the increment in absorbed dose in 
time Af; i.e., 


Absorbed dose rate = ~. (2.131) 


A special unit of absorbed dose rate is any quotient of 
the rad or its multiple or submultiple by a suitable unit 
of time (rads/sec, mrads/hr, etc.). 

Energy can be imparted to a volume of matter by 
many different particle reactions, and particles having 
the same initial energy do not necessarily deposit the 
same amount of energy. This, of course, is because the 
energy deposition depends not only on the initial 
energy but also on the type of radiation and the kinds 
of interactions that occur. Schematic representations of 
energy imparted to the matter contained within a 
volume element AV for two particles having initial 
energies E§ and £3 are shown in Fig. 2.16 (superscripts 
1 and 2 identify the reactions produced by particles | 
and 2 respectively). The corresponding equations for 
the energy imparted are 


OB, =£}-E}-E!-@!, (2.132a) 


AE, = EG -~E+Q? : (2.1325) 
where AF’) is the energy imparted to the matter within 
the volume element AV, Q is the energy equivalent 
of any change in the rest mass due to nuclear or 
elementary particle reactions within the volume, and 
subscripts 6 and y refer to the type of particle leaving 
the volume. In reaction 1, an incident gamma ray 
undergoes a Compton scattering within the volume 
element AV. Q’ is the binding energy of the Compton 
electron which is usually of negligible magnitude. The 
Compton electron loses some of its energy through 
ionization within the volume and then departs with the 
energy ER. The scattered (degraded) gamma ray also 
leaves the volume. In the second reaction, a neutron 
undergoes radiative capture and a gamma ray is pro- 
duced which leaves the volume. Q? is the binding 
energy associated with the neutron capture which 
appears as excitation energy. Note that no secondary 
collisions within the volume element were included in 
this schematic representation. This was purposely done 
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Fig. 2.16. Schematic Representation of Relationship Between 
Absorbed Dose (Energy Deposition) and Kerma. 


to illustrate that the volume element AV should be 
large enough for many primary interactions to occur 
but sufficiently small so that a primary particle and/or 
secondary neutrons or gamma rays usually do not suffer 
collisions subsequent to the initial collision by the 
incident particle within the volume element AV. 

The average energy imparted to the matter within the 
volume element AV by all reactions of particles of a 
particular type and energy can be written as 


2 Ri (E) EE) 


Ino (2.133) 


E,(E) = 


where 


E 'p(E) = the average energy imparted in ergs or MeV, 


R,(E)= rate of the ith type of reaction for particles 
of energy £ within the volume element AV, 

E,(E) = energy imparted, in ergs or MeV, associated 
with a particle of energy EF undergoing 
reaction i within the volume element AV. 


2.9.2. FIRST-COLLISION DOSE AND KERMA 


The quantity absorbed dose discussed in the previous 
section is a physical variable that is closely related to 


the potential hazard of a radiation field but for a real 
situation is sometimes very difficult to calculate accu- 
rately or to relate to the response in a human being. 
Consequently, other concepts have often been used. 

One of the concepts used is the first-collision dose 
(also called the single-collision dose). This quantity has 
been subject to several interpretations that differ from 
one another in subtle ways. One of the most compre- 
hensive reviews of the meanings of first-collision dose 
has been written by Attix,>® who says that the 
first-collision dose due to a radiation field (neutrons 
and/or gamma rays) can be interpreted in each of the 
following three ways: 

1. The first-collision dose is the absorbed dose 
contributed by all particles incident on an isolated small 
mass (unless the term is explicitly modified to include 
only one component, e.g., neutrons) in which charged- 
particle equilibrium* exists. 

2. The first-collision dose is the absorbed dose 
contributed by primary radiation only in a small mass 
that can be located anywhere and in which charged 
particle equilibrium exists. 

3. The first-collision dose is the energy transferred 
from uncharged incident radiation, regardless of origin, 
to the charged particles formed in a limitingly small 
probe that may be located anywhere. Charged-particle 
equilibrium in the probe is not a requirement. 

The small mass specified in interpretations 1 and 2 
means that the probe is small enough to leave the 
radiation field unperturbed and the probability is 
negligible that the incident particles will interact with 
the probe more than once or that secondary gamma 
rayst produced within the mass will be absorbed 
within it. The small probe specified in interpretation 3 
has the same requirements but it can be smaller than 
that required for the other interpretations because 
charged-particle equilibrium is not necessary. 

It can be stated at this point that interpretation 2 is a 
special and seldom used definition that will not be 
elaborated on here. Interpretation 1, on the other hand, 
is a widely used definition — particularly by experi- 
mentalists — and is often referred to as the free-field 
dose, the air dose, the free-air dose, etc. But in an 


*Charged-particle equilibrium may be viewed as that condi- 
tion when on the average as many charged particles such as 
electrons enter the volume element AV as leak out, thereby 
resulting in an essentially zero net transfer of energy by the 
flow of electrons. 

¥ As used here, secondary gamma rays refer to those gamma 
rays produced by the interaction of the incident radiation (e.g., 
capture gamma rays and inelastic-scattering gamma rays). 
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attempt to alleviate confusion, the ICRU chose inter- 
pretation 3 to use as a basis for an exact definition and 
called the quantity so defined the kerma (kinetic energy 
released in material). In fact, the ICRU encourages the 
exclusive use of the term kerma’ with the implication 
that the term first-collision dose and its other interpre- 
tations should be relegated to the archives. 

The definition of kerma, K, as recommended by the 
ICRU? is 

K=AE,/AM, (2.134) 
where AE, = sum of the initial kinetic energies (ergs) 
of all the charged particles liberated by indirectly 
ionizing particles in a volume element of material AV, 
and AM = mass (g) of the material contained in the 
volume element AV. 

The kerma rate is the quotient of AK by Af, where 
AK is the increment in kerma in time At: 

Kerma rate = AK/At. (2.135) 

Since AF, is the sum of the initial kinetic energies of 
all charged particles liberated by indirectly ionizing 
particles, it includes not only the kinetic energy these 
charged particles expend in collisions but also the 
bremsstrahlung they radiate and the energy of any 
charged particles that are produced in secondary 
processes. 

The kerma or kerma rate for a specified material in 
free space or at a point inside a different material will 
be that which would be obtained if a small quantity of 
the specified material were placed at the point of 
interest. 

With the definition of kerma in mind, the reactions 
shown in Fig. 2.16 are interpreted in terms of the initial 
kinetic energy imparted to charged particles: 


AEX =E} - E,-Q', (2.136a) 


AE? = £3 — £2 +9". (2.136b) 
It is noted that the expression for AE? and the ex- 
pression given in Eq. 2.132b for AE}, are the same, 
while the expression for AE; differs from that given in 
Eq. 2.132a for AE}. This is because the Compton 
electron produced by reaction 1 did not deposit all of 
its energy within the volume element and there was no 
compensating in-leakage of an electron; i.e., charged- 
particle equilibrium did not exist. 

Auxier*® strongly criticizes the ICRU for not for- 
mally recognizing first-collision dose as given by inter- 
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pretation 1. His objections are based on the fact that 
the first-collision doses obtained from measurements 
with ionization chambers, proportional counters, etc., 
are absorbed (imparted) energy and not transferred 
kinetic energy. Therefore, since measured doses cannot 
be interpreted as kerma, except by deduction, some 
dosimetrists are insisting on retaining the term first- 
collision dose. On the other hand, if charged-particle 
equilibrium exists, which it may or may not, and if 
bremsstrahlung radiation is ignored, then first-collision 
dose and kerma are equivalent. It can be seen that, in 
bulk shields or large tissue masses, first-collision dose 
and kerma have nearly the same magnitude but that in 
thin layers, such as clothing or skin, they can be quite 
different. 


The early calculations of first-collision dose per unit 
fluence were performed by Snyder*’ for mono- 
energetic neutrons incident on a four-element tissue 
model. Although the term kerma had not yet been 
introduced, the doses obtained by Snyder are those 
described by interpretation 3 above and thus are kerma 
per unit fluence. His calculations considered only 
neutron captures and elastic scatterings, with the latter 
assumed to take place isotropically in the center-of- 
mass system. Later, Henderson*? made similar calcula- 


‘tions for neutrons for a four-element tissue model and 


for ethylene. He reported his results as absorbed dose 
(in rads/hr per unit flux), but they are actually kerma 
rate per unit flux. 


The most recent and most comprehensive calculations 
of neutron kerma factors (kerma per unit fluence) were 
made by Ritts et al.4* They included essentially all 
reactions, considered anisotropic scattering in those 
reactions for which angular cross-section data were 
available, and in all cases used the most recent 
cross-section data. They used several models, one of 
which was an 1l-element “standard man” model in 
which all tissue organs and bone were homogenized and 
the other models were for specific body organs. The 
elemental compositions used in the calculations are 
shown in Table 2.8, and all the reactions considered are 
listed in Table 2.9. 


The calculations by Ritts et al. were made for 816 
discrete neutron energies corresponding to the super- 
group-subgroup structure of the O5R Monte Carlo 
code** over an energy range 0.023 eV to 19.2 MeV. 
The kerma factor for a particular irradiated material 
was found by summing the averages of the kinetic 
energies imparted to the struck nuclei and the energies 
associated with charged particles that were emitted. 


Table 2.8. Elemental Percentages and Number Densities Used in Neutron Kerma Calculations® 


Percent Composition Number Density (1024 atoms/g) 

nen ates Lung Muscle Bone Brain od, ae Lung Muscle Bone Brain ee. 
H 10 10 10 71 10 10 5.977(—2)©  6.062(—2) 5.977(—2) 4.50(—2) 5.977(—2) 5.977(—2) 
O 60 75 75 39.5 71 74° 2.259(—2) 2.864(—2) 2.824(—2) 1.49(—2) 2.674(—2) 2.7866(—2) 
Cc 24 10 11 22.4 15 12.5 1.204(—2) 5.09(-3) 5.518(—2) 1.13(—2) 7,52(—3) 6.270(—3) 
N 2.9 2.6 2.6 4.7 1.9 3.2 1.25(-3) 1.13(—3) 1.118(—3) 2.05(—3) 8.17(—4) 1.377(-3) 
Ca 1.2 0.011 0.0031 12,37 0.0086 1.77(-4) 1.68(—6) 4.66(—7) 1.87(—3) 1.29(-6) & 
P 1.1 0.11 0.18 13.16 0.34 2.14(-4) 2.17(-5) 3.50(—5) 2.57(—3) 6.61(—5) 
Ny} 0.24 0.24 24 0.41 0.17 4.51(-4) 4.57(-5) 4.51(-5) 7.69(—5) 3.19(—-5) 

0.20 0.20 0.30 0.30 3.08(-5) 3.13(-5) 4.62(-5) 4.62(—5S) 
Na 0.20 0.18 0.16 0.08 0.18 524(-5) 4.78(-5) 4.19(—S) 2.08(—5) 4.72(-5) 
Cl 0.20 0.25 0.18 0.23 0.12 3.40(-5) 4.31(-S) 3.06(—5) 3.91(-5) 2.039(—5) 
Mg 0.03 0.0096 0.019 0.28 0.015 7.43(-6) 4.78(—6) 4.71(-6) 6.87(—5) 3.72(-6) 


2From ref. 43. 
>The values for oxygen and carbon were increased by 1% and 0.5% respectively. 
Read: 5.977 X 10 “, etc. 
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Table 2.9, Reactions Considered in Calculations of Neutron Kerma Calculations® 


Type of Reaction H oO Cc N Ca P A K Na Cl Mg 
Radiative capture x x x x x x 
Elastic: scattering x x x x x x x x 
Inelastic scattering x xX x x x x x x 
(n,p) reaction x? x x x x x 
(n,Q) reaction x x xb x x? x 
(n,n'30) reaction 

(n,2n) reaction xe xe xe 

(n,d) reaction x x 

(n,t) reaction x 

(n, 20) reaction x 

(n,n'p) reaction x x x 
(n,n'0) reaction x 


“From ref. 43, 
includes beta decay after reaction, 
“Includes positron decay. 


Therefore the kerma per unit fluence may be expressed - 


as 
Mg = Ly CNAE) 0;;(E) Ey 5 (2.137) 


where 


N; = number density of element i (atoms/g), 


Gj; = microscopic cross section of reaction j for 
element i (cm?/atom), 


Ej; = average energy deposited by reaction j for 
element i (MeV), 


C= 1.602 X 107° (ergs/MeV). 


Since the kerma values are given per unit fluence, 
they are fluence-to-kerma conversion factors (usually 
referred to simply as kerma factors). Such kerma 
factors for the standard man model are plotted in Fig. 
2.17 and tabulated in Table 2A.! in Appendix 2A, 
along with kerma factors for the individual elements 
making the greatest contributions to the total factor. 
Kerma factors for specific body components (lung, 
muscle, bone, brain, and red marrow) are also tabulated 
in Appendix 2A, as well as unweighted kerma factors 
for the important elements in the standard man. With 
these values a total kerma factor for any arbitrary 
composition of these materials can be easily calculated. 

Kerma rate factors for monoenergetic gamma rays 
incident on carbon, air, and a four-element tissue model 
were calculated by Henderson*? and are plotted in Fig. 


2.18. (Henderson’s gamma-ray results, like his neutron 
results, were actually reported as absorbed dose, the 
term kerma not having been introduced yet.) Prior to 
Henderson’s work, exposure factors for air published by 
Rockwell*® were generally used as the equivalent 
flux-to-dose rate conversion factors for tissue. 


2.9.3 EXPOSURE 


Exposure is a term that should be used only for 
gamma rays. Formerly referred to as exposure dose, it 
came into common usage after the problem of speci- 
fying biological dose associated with X rays was first 
encountered. As recommended by the ICRU, exposure 
(X) describes the deposition of energy in air by 
electromagnetic radiation and is defined by 


X=AQ/MM, (2.138) 


where AQ = sum of the electrical charges of all the ions 
of one sign produced in air when all electrons liberated 
by photons in a volume of air AV are completely 
stopped,* and AM = mass of the material contained in 
the volume element AV. Exposure rate is given by 


*The ionization that would be produced by the brems- 
strahlung associated with secondary electrons is not included in 
AQ. Except for this small difference, which is insignificant for 
photon energies less than 15 MeV, exposure is equivalent to 
kerma. 
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AX/At, where AX is the increment ‘of exposure. during 
the increment of time Af. 

The unit of exposure is the roentgen (r). It is the 
quantity of X or gamma radiation that produces, in air, 
ions carrying 1 electrostatic unit of charge per 
0.001293 g of dry air (or 2.58 X 10~* coulombs/kg). In 
terms of ergs, the roentgen is equal to 87.7 ergs per 
gram of air or 96 ergs per gram of tissue. For gamma 
rays above 3 MeV in energy, the range of secondary 
electrons in low-atomic-number materials becomes com- 
parable with the relaxation length of the gamma rays. 
Consequently, the ionization produced in a small 
volume is no longer a sole measure of the intensity of 
the radiation at that point. The ICRU does not 
recommend the use of the roentgen above 3 MeV, but 
in practice, it is still used above 3 MeV with instruments 
calibrated in terms of energy absorption in air. 

Since 1 rad is equal to 100 ergs per gram of irradiated 
material (see Section 2.9.1), the roentgen and the rad 
are frequently interchanged when tissue exposure is 
referred to, the difference of 4% being no real conse- 
quence in shielding calculations. Strictly speaking, 
however, the rad should be reserved for describing an 
absorbed dose. 

Away from boundaries, the distinction between 
kerma and exposure loses relevance if only the three 
principal processes (Compton scattering, photoelectric 
effect, and pair production) are considered and if the 
photoelectric effect and pair production are assumed to 
be absorptions. For this simple and widely used model* 
the fluence-to-exposure conversion factor is propor- 
tional to the product of the photon energy and the 
energy-deposition coefficient, in which case exposure, 
absorbed dose, first-collision dose, and kerma in air are 
all equal in magnitude. 

In some calculations, however, a slightly more com- 
plex model is used in which only the photoelectric 
process is treated as an absorption. As in the above 
model, the total kinetic energy of the electron pair 
created in the pair-production process is assumed to be 
absorbed at the point of their emission. However, the 
two 0.51-MeV gamma rays produced by the annihila- 
tion of the positron are treated as scattered gamma 
rays, and pair production is viewed as a scattering 
process. This is described in terms of a modified 
pair-production cross section which is given by 


1.022 
lige: Wes (1 2 =) _ E> 1.022 MeV , (2.139) 


*For this and the subsequent model it is assumed that the 
bremsstrahlung radiation that would be produced by the 
secondary electrons (and positrons) is immediately absorbed, 
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where £ is the energy of the incident photon (MeV). 


2.9.4 RBE DOSE; DOSE EQUIVALENT 


Whenever the biological response of a human organ to 
radiation exposure is of concem, merely knowing the 
absorbed dose (energy deposited) is not enough. Bio- 
logical responses vary both with the nature and energy 
of the radiation and with the composition and function 
of the irradiated organs. Thus when an organ is exposed 
to a mixed radiation field or to a field of radiation 
comprised of one type of particle with various energies, 
the energy deposited by each type of particle of a given 
energy must be weighted by some factor before the 
total biological dose received by the organ can be de- 
termined by a summing of the individual contributions. 

For a specific biological effect in a particular mass the 
weighting factor is mainly a function of the Jinear rate 
of energy transfer (LET) to the system by charged 
particles set in motion by the interactions of incident 
radiation. The LET of a material is related to the linear 
stopping power but.the concepts are different. The 
linear stopping power is the average energy loss per unit . 
path length by a charged particle in traversing a medium 
regardless of where the energy is absorbed. The LET, 
however, refers to energy imparted within a limited 
volume. These two quantities are equal only in the 
special case when the LET includes the absorption of all 
secondary particles, in which case it is called LET,. 
These quantities are described more completely in ref. 
21 (p. 51). 

The LET-dependent weighting factor is called the 
relative biological effectiveness (RBE) and is defined as 


rads of standard radiation producing 
a given biological effect 
rads of another type of radiation * 
producing the same effect 


RBE = (2.140) 


The standard radiation referred to is X rays of 250-keV 
energy, and thus by definition the RBE for 250-keV X 
rays is 1. 

When a dose given in rads is weighted by an RBE 
value, the resulting dose, sometimes referred to as the 
RBE dose, is given in rems, a unit derived from the term 
roentgen equivalent man: 

RBE dose in rems = dose in rads X RBE. (2.1.41) 
It follows that for 250-keV X rays the rad dose and the 
rem dose are numerically equal. 


The rem is quantitatively defined as: the absorbed 
dose due to any ionizing radiation that has the same 
biological effectiveness as 1 rad of X rays with an 
average specific ionization of 100 ion pairs per micron 
of water in terms of its air equivalent. It turns out that 
X rays and gamma rays generally do not exceed this 
specific ionization (LET < 3.5 keV/u), which is 
considered to be the boundary between low-LET and 
high-LET radiation and thus to be the limiting condi- 
tion for RBE= 1. 

In the past RBE has been used in the fields of both 
radiobiology and radiation protection, but this gener- 
ated concern because of the differences in application 
and, to some extent, in concept. Consequently ICRU 
recommends that RBE be used only for correlating 
radiobiological experiments and that a new term, the 
quality factor (QF), be used in the field of radiation 
protection. Quality factors are actually those values of 
RBE that are intended to embrace all effects that are 
hazardous to human beings. In other words, QF values 
are not related to specific organs of the body as are 
some of the RBE values. When a QF value is used to 
weight absorbed dose, the resulting dose in rems is 
identified as the dose equivalent (DE). It is assumed 
that dose equivalents, like RBE doses, are additive. 

Although QF is defined in purely physical terms (that 
is, as a function of LET only, which in tum is a 
function of the number of ion pairs produced per 
centimeter of travel by charged particles), the basis of 
the legislated value is biological. Table 2.10 shows the 
recommendations*°**” of the International Committee 
on Radiological Protection (ICRP) for the relation 
between QF and both LET and ion pair formation. 

The ICRP recommends a QF of 1 for X rays, gamma 
rays, electrons, and positrons of any specific ioniza- 
tion:*® QF’s for neutrons, protons, and heavy recoil 
nuclei vary with energy. The ICRP recommendations 
for QF for neutrons between 0.01 and 1000 MeV are 
plotted in Fig. 2.19. Tabulated values for lower energy 
neutrons are also shown. These values are based on the 
calculated results of Snyder and Neufeld*® and of 
Neary and Mulvey.*? As a rough and conservative 
approximation, a QF of 10 is applicable to protons and 
a QF of 20 can be used for heavy-recoil nuclei. 

The preceding QF’s are assumed to be applicable for 
whole-body irradiation. When specific organs are being 
considered, additional modifying factors may be 
needed. For example, ICRP specifies that when the lens 
of the eye is being irradiated a modifying factor of 3 
should be used if the QF is equal to or greater than 10 
and a modifying factor of 1 should be used when QF = 
1. Strangely enough no recommendations are given for 
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1 < QF < 10. Presumably, linear interpolation is to be 
used between these limits. 

When the dose in an organ is nonuniform, a dose 
distribution factor (DF) should also be applied. These 
factors cannot be established with certainty and only 
one is in current use. It is the relative damage factor 
applied in calculating the dose equivalent in bone from 
internal radiation. 

In general the shielding community has accepted the 
use of QF and DE, as is reflected by the latest shielding 
literature. From a practical viewpoint the change 
merely amounts to replacing RBE by QF and calling the 
product of the quality factor and absorbed dose the 
dose equivalent (DE) instead of the RBE dose. Conse- 
quently no confusion should exist when the older 


Table 2,10. QF as a Function of 
LET and Ion Pair Formation? 


F Average Specific Ionization Average LET 
;: (ion pairs/u) (keV/u?) 
1 <100 3.5 
2 200 7.0 
5 650 23 
10 1500 53 
20 2000 175 


"From refs. 46 and 47, 
bof water. 
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literature is consulted. It should be kept in mind, 
however, that QF values are chosen by such groups as 
the ICRP and thus are more subject to change than 
would be an RBE value based on experimental data. 

It should also be pointed out that all QF’s and other 
modifying factors are intended solely for chronic 
exposure to low-level radiation fields, with genetic 
damage being the hazard considered. High accidental 
exposures must be assessed on the basis of particular 
circumstances. The ICRP has no recommendations for 
such exposures, but in many cases the absorbed dose 
will give a better indication of the biological risks than 
will the dose equivalent. In other words, for acute 
effects due to massive exposures, the QF should be 
taken as unity, which partially results because energy 
may be wasted from the standpoint of the production 
of biological effects. Consequently for military applica- 
tions during war, when early death or incapacitation is 
the effect to be considered, the dose in rems may be 
taken as numerically equal to the absorbed or physical 
dose until better information becomes available. 


2.9.5 MAXIMUM ABSORBED DOSE; 
MAXIMUM DOSE EQUIVALENT 


As is pointed out in the introduction to this section, 
the usual shielding calculation provides a detailed 
description of the unperturbed radiation field (free 
field). That is, perturbations of the field by a human 
body are not considered. This greatly simplifies the 
analysis without seriously compromising the overall 
accuracy. But the results must somehow be related to 
the hazards to a human being if he were exposed to 
such a radiation field. 

A correlation between the unperturbed radiation field 
and the dose in a human body has been accomplished 
through the use of slab and cylindrical phantoms, which 
have dimensions and compositions resembling those of 
the human body. Calculations (or measurements) of the 
absorbed dose or the biological dose were made as a 
function of depth in the phantom for given mono- 
energetic neutrons or gamma rays incident on the 
phantom. The incident angular distribution was a 
monodirectional beam or an isotropic flux. (Calcula- 
tions for neutrons include the contribution to the dose 
by secondary gamma rays.) The maximum values in-the 
depth-dose distributions were identified as the max- 
imum absorbed dose and maximum dose equivalent for 
the absorbed (physical) and biological doses respec- 
tively. Since these values were obtained and reported on 
the basis of a unit particle entering the phantom, they 
can be used as response functions to relate an unper- 
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turbed field of mixed radiation to the maximum dose 
that would be received by some part of the body. Use 
of these maximum values is dictated by a conservative 
design philosophy that does not allow the permissible 
dose to be exceeded at any point in the body. 


Neutron Doses in Phantoms 


The first dose calculations in a phantom were 
performed by Snyder and Neufeld*® for a beam of 
monoenergetic neutrons normally incident on a phan- 
tom represented by a slab of tissue. The slab, infinite in 
the transverse directions, was assumed to be 30 cm 
thick. Snyder and Neufeld determined the distribution 
of the absorbed dose through the slab, finding that the 
maximum dose occurred at nearly the surface or within 
a few centimeters of the surface. In these calculations, 
only neutron captures and elastic scatterings were 
considered, which limited the contribution by second- 
ary gamma rays to capture gamma rays. The resulting 
maximum absorbed doses and biological doses as a 
function of the incident neutron energy are shown in 
Table 2.11. These values have been widely used in 
reactor shield design to convert neutron fluences to 
doses and have become a virtual standard since these 
dose equivalents have been stipulated for use by the 
Federal Register.°° 


Table 2,11. Maximum Absorbed Dose and Maximum 
Dose Equivalent for Monoenergetic Neutrons 
Incident on a Slab Tissue Phantom? 


Maximum Maximum 
Neutron Absorbed Dose Effective 
Energy Dose? Equivalent? QF 
(MeV) (mrad (mrem 
neutron! cm~?) neutron! cm~?) 

Thermal 3.2(-7)° 1,04(-6) 3.25 
0.0001 6.9(—7) 1.39(—6) 2.01 
0.005 5.7(—7) 1.22(-6) 2.14 
0.02 5.7(-7) 2.35(—6) 4.12 
0.1 1.10(—6) 8.3(-6) 7.55 
0.5 2.4(-6) 2.30(—5) 9.61 
1,0 3.8(—6) 3.80(—5) 10.0 
2.5 4.3(-6) 3.41(-5) 7.93 
5.0 5.8(-—6) 3.80(—5) 6.55 
15 7.1(-6) 4.16(—5) 5.85 

10 7.0(-6) 4.16(-5) 5.94 


°From ref, 48, 
>Multiply by 3600 to convert to (mrad/hr)/(neutron/cm?-sec) 
or to (mrem/hr)/(neutron/cm?'sec). 


Read: 3.2 X 1077, etc. 


In similar calculations, performed later by Auxier et 
al,°! the infinite slab was replaced by a finite 
cylindrical phantom that more nearly mocked up a 
human body. A schematic of the phantom, 60 cm high 
and 30 cm in diameter, is given in Fig. 2.20. The 
maximum doses (averaged for a volume element) 
occurred in the outermost volume element at the 
midplane on the side facing the beam. These maximum 
doses are given in Table 2.12, and the doses for each 
volume element are presented in Appendix 2B. 

The effective quality factors shown in Tables 2.11 
and 2.12 are the ratios of the maximum dose equivalent 
to the maximum absorbed dose for a given incident 
neutron energy and phantom configuration. The magni- 
tude of the effective quality factor is very close to but 
usually slightly less than the QF’s recommended by the 
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Table 2,12. Maximum Absorbed Dose and Maximum 
Dose Equivalent for Monoenergetic Neutrons 
Incident on a Cylindrical Tissue Phantom? 
Maximum Maximum 
Neutron Absorbed Dose Effective 
Energy Dose? Equivalent? F 
(MeV) (mrad (mrem Q 
neutron™! cm7*) _neutron™! cm~?) 

Thermal 4.68(—7)° 1.15(—6) 2.46 
0.000001 5.89(-7) 1.34(-6) 2.27 
0.00001 5.18(—7) 1.21(-6) 2.34 
0.0001 4.45(-7) 1.01(-6) 2.27 
0.001 4.32(-7) 8.85(—6) 2.04 
0.01 4.34(-7) 9.92(-6) 2.29 
0.1 8.02(—7) 4,86(-6) 6.06 
0.5 1.81(—6) 1.89(—5) 10.4 
1 3.01(-6) 3.26(—5) 10.8 
2.5 3.99(—6) 3.50(—5) 8.77 
5 5.72(—6) 4.41(-—5S) 7.70 
7 5.70(—6) 4.03(-S) 7.07 

10 7.25(—6) 4.31(-S) 5.94 

14 8,.31(-6) 6.15(—5) 7.40 


“From ref, 51. 

>Multiply by 3600 to convert to (mrad/hr)/(neutron/cm?'sec) 
or to (mrem/hr)/(neutron/cm?-sec). 

Read: 4.68 X 1077, etc. 


ICRP (see Fig. 2.19). The differences in the values 
would be expected since the QF’s recommended by the 
ICRP consider only the initial collision of the neutron 
whereas the effective QF’s represent an average of all 
collisions experienced by the neutron within the phan- 
tom. 

The Snyder-Neufeld calculations of maximum ab- 
sorbed dose have recently been improved upon by Ritts 
et al.** They increased the number of elements in the 
infinite-slab phantom from 4 to 11, considered a much 
larger number of neutron source energies, allowed for 
neutron inelastic scattering (and the associated ine- 
lastic-scattering gamma rays) and for a number of other 
nonelastic events, corrected the elastic scattering to 
linear anisotropy, used more accurate secondary- 
gamma-ray spectra, used the latest cross sections avail- 
able, and handled the neutrons with energies below 1 
eV more accurately. The percentages of the 11 elements 
in the slab were the same as those used for the standard 
man kerma calculations (see Table 2.8), and the 
reactions considered were also the same as for the 
kerma calculations (see Table 2.9). The resulting maxi- 
mum absorbed doses are presented in Fig. 2.21 (solid 
line) and in the second column of Table 2.13. A 
comparison with the Snyder-Neufeld results of Table 
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Table 2,13. Neutron Maximum Absorbed Doses Within 
Standard-Man Slab Phantoms? 


Dose [(rads/g)/(neutron/c m?)] b 


Source 


Beam Source Isotropic Flux 
Energy a 
(eV) Infinite Finite Infinite Finite 
Slab Slab Slab Slab 


1.4208(+7)4 7.9598(-—9) 7.4605(—9) 1.4616(—8) 1,3638(—8) 
1.2856(+7) 7.3683(—9) 6.8683(—9) 1.3436(—8) 1.2553(—8) 
1.1633(+7) 7,2273(—9) 6.7712(-9) 1.3077(—8) 1.2274(—8) 
1.0526(+7) 6.8462(—9) 6.4060(—9) 1.2342(—8) 1,1592(-—8) 
1.0000(+7) 6.7111(-9) 6.2789(—9) 1.2070(—8) 1.1341(-8) 


9.5242(+6) 6.5745(—9) 6.1445(—9) 1.1798(—8) 1.1091(-8) 
8.6178(+6) 6.3415(—9) 5.9153(—9) 1.1308(-8) 1.0641(-8) 
7.7977(46) 6.3362(—9) 5.9281(—9) 1.1207(-8) 1.0574(-8) 
7.0557(+6) 6.0811(—9) 5.6837(-9) 1.0736(—8) 1.0136(-8) 
6.3843(+6) 5.8372(—9) 5.4736(—9) 1.0319(-8) 9.7712(—9) 


5.7767(+6) 5.7003(-—9) 5.3445(-9) 1.0002(—8) 9.2060(—9) 
5.2270(+6) 5.4689(—9) 5.1380(-9) 9.6108(—9) 8.8640(—9) 
5.0000(+6) 5.4174(-9) 9.4514(—9) 

4.7296(+6) 5.3667(—9) 5.0331(—9) 9.0074(-9) 8.5870(—9) 
4.2795(+6) 5.2731(-9) 4.9584(-9) 8.8140(—9) 8.4324(—9) 


3.5038(+6) 5.2205(-9) 4.9916(-9) 8.5412(-9) 8.2110(—9) 
3.1703(+6) 4.7108(—9) 4.4959(-9) 7.7798(—9) 7,.4682(—9) 
2.5956(+6) 4.1653(—9) 3.9630(—9) 6.8302(—9) 6.5576(—9) 
2.0000(+6) 3.8389(—9) 6.1462(—9) 

1.9229(+6) 3.7840(—9) 3.6060(—9) 6.0154(—9) 5.7972(—9) 


1.5743(+6) 3.4752(—9) 3.3079(-9) 5.4326(-9) 5.2392(—9) 
1.0553(+6) 3.1381(—9) 2.9991(—9) 4.6764(—9) 4.5326(—9) 
5.2405(+5) 2.1666(—9) 2.0390(—9) 2.9434(—9) 2.8266(—9) 
5.0000(+5) 2.1724(—9) 2.0494(—9) 2.9096(—9) 2.7980(—9) 
4.7418(+5) 2.1766(—9) 2.0577(-9) 2.8652(—9) 2.7590(-9) 


4.2906(+5) 2.1845(—9) 2.6712(-9) 

1.9279(+5) 1.3014(—9) 1.1931(-9) 1.5449(-9) 1.4506(—9) 
9.8803(+4) 9.2925(-10) 8.2811(—10) 1.0393(—9) 9,.5232(—10) 
4.667,1(+4) 6.5317(-10) = 5.5685(-10)  6.6372(-10) = 5.7892(—10) 
2.8308(+4) 5.3582(-10) 4.4115(-10)  5.1846(-10) 4.3284(—10) 


1.0414(+4) 4.5827(-10)  3.5333(-10) 4.2532(-10) —_3.2792(—-10) 
4.9191(+3) 4.6004(-10) 3.5522(-10)  4.2712(-10) —3.2960(—10) 
1.0976 (+3) 5.2192(-10)  4.0294(-10) 4.3888(-10) 3.4036(—10) 
3.1447(+2) 5.3477(~10)  4.1505(-10)  4.4894(-10) —3.4994(-10) 
9.0097(+1) 5.4829(~10) 4.2781(-10) 4.5880(-10) 3.5930(-10) 


3.3145(+1) 5.5774(—10) 4.3691(—10) 4.6500(—10) 3.6674(—10) 
9.4962(0) 5.6353(-10) 4.4345(-10) 4.7006(-10) = 3.7620(—10) 


4.4857(0) - 5.6082(~10) 4.4543(-10) 4.7066(-10)  3.7836(—10) 
2.1189(0) §.5112(-10) 4.4173(-10)  4.6484(-10) = 3.7528(-10) 
1.0009(0) 5.2998(—10)  4.2684(-10)  4.4830(—10) 3.6328(—10) 


6.0707(-1) 5.0478(-10) 4.0796(-10)  4.2706(-10) — 3.4692(—10) 
4.7279(-1) 4.8934(-10)  3.9628(-10)  4.1224(—10) 


2.1800(—1) 3.5443(—10) 2.9505(—10) 1.9709(-10) + 1.6567(—10) 


2From ref. 43. 
>These values correspond to one neutron per cm? entering the slab. 


©The calculations for the finite slab include a buckling term to compensate for 
transverse leakage; the transverse dimensions for the finite slab phantom were 
taken as 40 by 61 cm. 

Read: 1.4208 X 107, etc. 
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Fig. 2.21. Neutron Maximum Absorbed Doses in Standard-Man Slab Phantom Due to a Beam Source. (From Ritts et al, ref. 43.) 


2.11 reveals that the new values are consistently lower, 
as would be expected since the Snyder-Neufeld calcula- 
tions were purposely conservative. 

In order to determine the effect that transverse 
leakage would have on the calculated doses, Ritts et al. 
also performed one-dimensional calculations for a 
30-cm-thick slab with a geometric buckling* value that 
was equivalent to a length of 61 cm and a width of 40 
cm. In general they found that a 5 to 10% reduction in 
the dose could be attributed to transverse leakage. The 
results for the finite slab are included. in Fig. 2.21 
(dashed line) and in the third column of Table 2.13. 
These values for the finite slab are 10 to 20% lower 
than the doses calculated by Snyder and Neufeld for 
the slab with no transverse leakage. 

Ritts et al. also made calculations with and without a 
buckling value for an isotropic neutron flux incident on 
a slab. The results are presented in Fig, 2.22 and in the 
last two columns of Table 2.13. 


*The geometric buckling, B. is the lowest eigenvalue that 
results from solving the wave equation for a system of particular 
size and shape subject to the condition that the flux vanish at 
the extrapolated boundary. In one-dimensional calculations, 
leakage in the infinite directions can be approximated by using 
DB? (D = diffusion coefficient) as a fictitious macroscopic 
absorption cross section. 


Gamma-Ray Doses in Phantoms 


The doses delivered to various parts of a man-like 
phantom by 0.027- to 1.24-MeV gamma rays were 
measured by Jones,°? who later reported his results in 
terms of calibration factors®*? for dosimeters. The 
calibration factors were recommended for use in the 
design and calibration of gamma-ray exposure meters to 
permit a direct reading of the maximum dose (in rads) 
delivered to the critical organ for the particular environ- 
mental condition. The composition of the phantom 
matched soft tissue with respect to density and effec- 
tive atomic number and included a full set of bones for 
the head and trunk. 

For each of 11 photon energies used in the experi- 
ment (see Table 2.14), measurements were made of the 
free-field dose’ (the dose with the phantom removed) 
at a distance from the source corresponding to a 
particular dose position in the phantom. Measurements 
within the phantom were made at positions corre- 
sponding to the location of the brain (front), eye (lens), 
skin (front and back of chest and forehead), gut 
mucosa, ovaries, testes, and bone marrow (seven repre- 
sentative sites). For each source energy the measure- 


This is effectively a measurement of the first-collision dose 
as described in Section 2.9.2. 
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Fig. 2.22, Neutron Maximum Absorbed Doses in Standard-Man Slab Phantom Due to an Isotropic Source. (From Ritts et al, ref. 


43.) 


ments were made with the phantom being irradiated 
from in front, from behind, and while being rotated 
about the vertical axis. Table 2.14 shows the location at 
which the maximum dose was received within the 
phantom for each exposure to gamma rays of different 
energies. The results are given in terms of the ratio of 
the absorbed dose to the free-field exposure (rads/r). 

Recent Monte Carlo calculations by Sidwell et al.°* 
for monoenergetic gamma rays incident on a phantom 
assumed to have the characteristics of water and to be 
shaped like an elliptic cylinder yielded dose results that 
are reported to be in good agreement with Jones’ 
measurements. 

Other recent calculations of the dose delivered to a 
phantom by monoenergetic gamma rays were per- 
formed by Claiborne and Tmubey.** Using the discrete 
ordinates method, with some checks by the Monte 
Carlo method, they calculated the dose distributed in a 
30-cm-thick infinite-slab phantom having the standard 
man composition shown in Table 2.8. Their maximum 
dose rates, which occurred at nearly the surface or 
within the first 2 cm of the surface, are compared in 
Fig. 2.23 with Jones’ experimental results (converted to 
the same units). First-collision dose rates for the 
standard-man composition based on the use of flux-to- 
kerma-rate conversion factors are also plotted in Fig. 
2.23. These values show that the customary practice of 
using kerma-rate factors underestimates the maximum 
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Fig. 2.23. Maximum Gamma-Ray Dose Rate in Slab Phantom 
of Standard-Man Composition. Comparison with kerma calcula- 
tions and with phantom measurements. (From Claiborne and 
Trubey, ref. 55.) 
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Table 2,14. Gamma-Ray Maximum Absorbed Doses at Various 
Locations in a Man-Like Phantom? 


Ratio of 
Gamma-Ray : 
bana Maximum Dose 
Source Energy Location in Phantom 
(MeV) to Exposure 
(rads/r) 
60 . b 
Co 1,24 Skin on back’ 1.06 
1370. 0.660 Skin on chest 1.11 
19215 0.400 Skin on forehead and testes 1,19 
Bremsstrahhing 
Narrow band 0.240 Skin on forehead 1.30 
0.165 Gut mucosa 1.47 
0.125 Testes 1.57 
Filtered by U 0.100 Skin on forehead 1.65 
Filtered by Pb 0.075 Gut mucosa 2.03 
Filtered by Ta 0.060 Skin on chest 1.75 
Filtered by Nd 0.038 Skin on chest 1.34 
Filtered by Sn 0.027 Skin on chest 1.28 


From ref. 52. 


Maximum occurred when phantom was irradiated from behind; all other 
maxima occurred when phantom was irradiated from in front. 


dose delivered to a position within a body, especially 
for the lower gamma-ray energies. This is because the 
energy deposition due to collisions subsequent to the 
first collisions is not accounted for in the kerma-rate 
factors,* 

The excellent agreement between the Claiborne- 
Trubey calculations and the Jones experiment for 
energies greater than 0.07 MeV seems surprising in view 
of the large differences in the configurations. However, 
since most of the energy deposition apparently comes 
from the first collision, the effect due to the geometry 
difference will be small. 


2.9.6. MULTICOLLISION DOSE 


The term multicollision dose, which until recently has 
been used solely for neutron exposure, is not recog- 
nized by the ICRU. The multicollision dose concept 
was developed by health physicists in an effort to relate 
the neutron flux or fluence incident upon a configura- 
tion of tissue representing the human body to the 
absorbed dose or to the dose equivalent. The multi- 
collision dose could, in turn, be related to a maximum 


*The dose rate in a slab-tissue phantom due to high-energy 
photons has been calculated by Alsmiller and Moran,°® but 
except for the values for 10-MeV gamma rays their results are 
for photon energies not encountered in weapons radiation or 
reactor shielding. 


permissible dose (MPD) established by shield design 
criteria. The term has tacitly been defined by common 
usage in shield design as the sum of the maximum dose 
equivalents (or maximum absorbed doses) delivered by 
each neutron energy group to a slab of tissue of a 
particular composition.* This includes the dose due to 
secondary gamma rays produced by the neutrons in the 
slab. This dose is calculated by 


Multicollision dose = J) Wp Dmax(Eg), (2-142) 
& 


where y, is the group fluence as defined in Section 2.1 
and DyaxEg) is the maximum dose equivalent (or 
maximum absorbed dose) in the slab for an incident 
neutron of the energy corresponding to the gth neutron 
group. 

Gamma-ray multicollision doses are similarly calcu- 
lated with Eq. 2.142, with, of course, We and Dmax(Eg) 
being the appropriate values for gamma rays. Because 
the gamma-ray maximum absorbed doses reported in 
Section 2.9.5 have only recently become available, most 
multicollision doses for gamma rays have in the past 
been determined by using the product of the gamma- 
ray energy, the mass energy deposition coefficient for 
tissue or air, and a factor to convert to rads for the 
Dyrnax (Eg) term. 


*See slab tissue phantom described in Section 2.9.5 for 
Snydes-Neufeld caiculations.*® 
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The total dose for a mixed radiation field incident on 
a tissue slab is determined by adding the neutron and 
gamma-ray multicollision doses obtained separately for 
the slab. 

Although the concept of multicollision dose, or dose 
rate, is not recognized by the ICRU, it seems clear that, 
to fully implement the trend to preciseness in the fields 
of radiobiology and radiation protection, the dose 
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calculated in shield design should be given recognition 
by some official body. Perhaps the term multicollision 
dose can be retained for the neutron-induced contri- 
bution and the term total multicollision dose can be 
used when the incident gamma-ray contribution is also 
included. Or perhaps an entirely new term could be 
introduced, such as the term man equivalent dose 
(MED) proposed by Claiborne and Trubey.°* 


2.10. Radiation Exposure Limits * 


Any material will undergo some alteration on ex- 
posure to directly or indirectly ionizing radiation, but 
organic materials are most susceptible to damage, 
particularly living organisms. The earth and its inhabit- 
ants are continuously exposed to radiation from 
cosmic-ray bombardment and naturally occurring radio- 
isotopes. In addition to these natural sources, the 
human body may be exposed to X rays used for 
medical diagnostics and treatments, lumindus wrist 
watch dials, television, etc. To date the effects of these 
exposures, if any, are essentially unknown. 

The International Commission on Radiological Pro- 
tection (ICRP) has taken the position that radiation in 
any amount absorbed by the gonads can cause genetic 
damage and increase the natural mutation rate. Mostly 
as a result of the fear that small chronic dosages will 
significantly increase the mutation rate in humans or 
induce unknown deleterious long-range effects (such as 
leukemia), the ICRP has recommended that workers in 
the nuclear industry be limited to a maximum dose of 5 
rems/year, which amounts to 2.5 mrems/hr based on a 
40-hr work week for 50 weeks/year. This so-called 
maximum permissible dose (MPD) is about 25 times the 
dose that the average American will receive each year 
from environmental radiation. The recommended MPD 
appears to represent an acceptable risk since con- 
servative estimates indicate that 5 rems/year may 
increase the natural gene mutation rate by only 10 to 
20%, which when related to the probable increase in the 
number of potentially defective children becomes a 
statistical exercise. It is interesting to note that no 
measurable defects attributable to radiation were found 
in studies of the children conceived after the parents 
were exposed to intense radiation fields from the 
atomic weapons exploded over Hiroshima and Nagasaki 
and to the Marshall Island natives exposed to fallout 
from the Bikini test shot. 

Much higher doses are required to produce significant 
somatic damage or direct deleterious effects in the 


*This section is largely quoted from Sections 2-8.4 and 2-8.5 
of Volume I of the Office of Civil Defense’s manual for shelter 
design and analysis (see ref. 57). 
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body. And since the body can repair itself, dose rate is a 
factor, a larger total radiation dose being required to 
produce a given degree of injury if the dose is spread 
over a long period of time. But the rate, or fraction- 
ation, effect is of no consequence for acute doses that 
could be received in a nuclear war or from an accident 
in the nuclear industry. Military considerations during 
war require emphasis on minimizing acute exposure and 
the attendant probability of an incapacitation that 
would prevent a necessary mission. Consequently the 
dose delivered to mucosal layers of the gut becomes the 
most important consideration. 

Although data on the biological effects that might be 
expected for human beings subjected to acute doses of 
varying amounts are rather limited, a sufficient amount 
of knowledge is available for drawing useful conclusions 
for planning purposes. Such data have been derived 
from the study of Japanese casualties, a few nuclear 
accidents, and large animal experiments. The following 
paragraphs give a general description of the effects that 
might be expected as the result of whole-body exposure 
to an acute dose received by healthy adults over a 
period of up to two to four days. For a more complete 
and detailed discussion see ref. 58. 

Doses up to 100 rems. — Exposures resulting in doses 
up to 100 rems would give no apparent immediate 
effects. In some cases nausea or vomiting may occur, 
but no other symptom of radiation sickness would be 
apparent. No clinical surveillance or therapy other than 
reassurance would be required. 

Doses from 100 to 200 rems. — Doses between 100 
and 200 rems would cause nausea and vomiting within 
about 3 hr in from 5 to 50% of those exposed. A 
moderate decrease in the number of white blood cells 
would result, making clinical and hematologic surveil- 
lance desirable. Therapy would consist mainly of 
reassurance. Prognosis would be excellent with a con- 
valescent period of several weeks. There should be no 
deaths. 

Doses from 200 to 600 rems. — Doses in the range 
200 to 600 rems would result in a 100% incidence of 
vomiting, which would begin within about 2 hr for 
doses greater than 300 rems. Severe loss of white blood 
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cells (leukopenia) would occur. Purple or livid spots 
(purpura) would appear on the skin. There would be 
some hemorrhaging and an incidence of infection. Loss 
of hair (epilation) would result for exposures above 300 
rems. Therapy would include blood transfusions and 
the administration of antibiotics. A convalescent period 
of from 1 to 12 months would be required. Deaths 
from hemorrhage or infection within two months 
would probably begin at about 250 rems and, for the 
dose range of 200 to 600 rems, would range up to 
about 80%. A value of around 300 rems delivered to the 
mid-line of the body is usually taken as LD. 9/39 (lethal 
dose for 50% within 30 days) for man. 

Doses from 600 to 1000 rems. — Doses approaching 
1000 rems would result in all the symptoms of 
radiation sickness described in the preceding paragraph 
but at a more severe intensity. Therapy would be only 


promising and would involve not only blood transfusion’ 


and antibiotics but also the consideration of bone 
marrow transplants. The incidence of deaths would 
range between 80 and 100%. 

Doses Greater Than 1000 rems. — Doses greater than 
1000 rems would almost certainly be lethal. Symptoms 
of radiation sickness would be very severe and would 
begin within 30 min. Those symptoms described in the 
preceding paragraphs would be augmented by diarrhea, 
fever, and disturbance of electrolyte balance as the 
result of changes in the gastrointestinal tract. Doses 
over 5000 rems would affect the central nervous system 
with attendant symptoms of convulsion, tremors, loss 
of muscular coordination, and lethargy. Prognosis 
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would be hopeless and therapy would consist mainly of 
the administration of sedatives. For doses of 1000 to 
5000 rems, death in almost 100% of the cases would 
result within 2 weeks from circulatory collapse. For 
doses above 5000 rems death in 100% of the cases 
would occur within 48 hr as a respiratory failure or 
edema of the brain, or both. Doses of 10,000 rems and 
up directly affect the nerve cells to the point that 
incapacitation could occur almost immediately. 

To summarize, if a person is exposed to a dose of less 
than about 200 rems, he should not become incapaci- 
tated nor should his ability to work be seriously 
affected. Persons exposed to doses exceeding about 200 
rems will suffer increasing radiation sickness with 
increasing dosage, and the probability of death is 
extremely high if the dose exceeds 600 rems. For doses 
greater than 10,000 rems immediate incapacitation can 
be expected. 

The MPD for military purposes during war will 
depend on the assigned mission and the acceptable 
casualty rate. Obviously the low MPD set for peacetime 
occupational work cannot be used as the criterion in 
military situations. Based on the available information, 
an MPD = 100 rems (absorbed over a few days) seems 
reasonable since the effect of such a dose would be 
hardly noticeable and would not interfere with military 
operations. Since radiation shelter design will require an 
estimate of the free-field dose, the maximum probable 
exposure in a specified time (which includes multiple 
strikes) should be used as the design criterion for each 
shelter in a given area. 


Appendix 2A. Neutron Fluence-to-Kerma Conversion Factors 
for Standard Man and Some Specific Body Components 


This appendix presents neutron fluence-to-kerma con- 
version factors calculated by Ritts et al.*? for a 
standard-man model and for specific body components 
(see Section 2.9.2). The neutron energies range from 
approximately 0.022 eV to approximately 15 MeV, and 
the compositions for the models are those shown in 
Table 2.8. Table 2A.1 gives the conversion factors for 
the standard man, together with those for the elements 
making the greatest contributions in the model. Con- 
version factors for specific body components are 
presented in Table 2A.2. 

Table 2A.3 gives the conversion factors on an atom 
density basis for the most common elements contained 
in man. These values can be used to obtain the kerma 
factors for any arbitrary combination of these elements 
by multiplying the individual value of each component 


by the number of atoms of each element per gram of 
mixture. For example, the kerma factor for H,O (11.9 
wt % H and 88.81 wt % 0) in a field of 15-MeV 
neutrons is 


0.119 7 . | 
{| To0g * (0-6023 x 10 ) (7.75 X 107°) 


0.881 
+|——— x (0.602 24) (6. 6 10724 
[ess (0.6023 X 1024) (6.67 X 107°) Je 0 


=7.41 X 1077 (ergs/g of H,O)/(neutron/em”) , (Al) 


where the numbers 1.008 and 16 are the atomic weights 
of hydrogen and oxygen respectively and 0.6023 X 
107* is Avogadro’s number. 


Table 2A.1. Neutron Fluence-to-Kerma Factors for Standard Man®? 


Energy Fluence-to-Kerma Factor (ergs/g)/(neutron/cm?) 
(MeV) Total Hydrogen Oxygen Carbon Nitrogen Phosphorus Sodium Chlorine 
15.14 7.079(—7)© 4.633(—7) 1.53(-7) 8.08(—8) 9.35(-9) 
14.84 6.948(—7) 4.629(—7) 1.46(—7) 7.53(-8) 9.13(-9) 
14.54 6.811(—7) 4.623(-7) 1.40(—7) 6.89(—8) 8.87(-9) 
14.24 6.652(—7) 4.617(-7) 1.29(-7) 6.47(-8) 8.67(—9) 
13.93 6.501(-7) 4.611(-7) 1.18(-7) 6.16(—8) 8.41(—9) 
13.63 6.276(-7) 4.604(-7) 9.90(—8) 5.93(—8) 7.99(—9) 
13.33 6.119(-7) 4.598(-7) 8.63(—8) 5.72(-8) 7.74(-9) 
13.03 6.007(-7) 4.591(-7) 7.87(-8) 5.44(—8) 7.63(—9) 
12.73 5.966(-7) 4.584(—7) 8.03(—8) 4.96(—8) 744(-9) 
12.43 5.949(-7) 4.577(-7) 8.30(—8) 4.61(—8) 7.29(-9) 
12.13 5.936(-7) 4.568(-7) 8.45(—8) 4.45(-8) 7.04(—9) 
11.83 5.964(-7) 4.559(-7) 9.28(-8) 4.01(-8) 6.77(-9) 
11.52 5.877(—7) 4.550(-—7) 8.90(-8) 3.64(—8) 6.47(—9) 
11.22 5.818(-7) 4.540(-7) 8.48(-8) 3.59(-8) 6.23(-9) 
10.92 5.691(-7) 4.531(-7) 7.49(-8) 3.43(—8) 5.95(-9) 
10.62 5.554(-7) 4.522(-7) 6.75({—8) 2.91(—8) 5.85(-9) 
10.32 5.483(—7) 4.512(—7) 6.81(—8) 2.26(—8) 5.73(-9) 
10.02 §.421(-7) 4.501(-7) 6.46(—8) 2.12(—8) 5.52(-9) 
9.72 §.353(-7) 4.478(—7) 5.98(—8) 2.18(—8) 5.28(-9) 
9.60 5.313(-7) 4.468(-7) 5.70(—8) 2.16(—8) 5.21(-9) 
9.45 §.271(—7) 4.456(-7) 5.40(-—8) 2.19(—8) 4.97(-9) 
9.30 §.222(-7) 4.443(-7) 5.00(—8) 2.23(—8) 4.91(-9) 
9.15 5.173(-7) 4.430(-7) 4.74(-8) 2.14(-8) 4.85(-9) 
9.00 5.123(-7) 4.417(-7) 4.78(-8) 1.75(-8) 4.70(-9) 
8.85 5.078(—7) 4.410(-7) 4.54(-8) 1.61(-—8) 4.64(-9) 
8.70 5.048(—7) 4.402(-7) 4.37(—8) 1.56(—8) 4.60(-9) 
8.62 5.007(-7) 4.399(-7) 4.01(—8) 1.56(—8) 4.53(-9) 
8.47 4.988(-7) 4.391(-7) 3.90(—8) 1.57(—8) 4.40(-9) 
8.32 5.036(-7) 4.383(-7) 4.38(—8) 1.66(-8) 4.28(-9) 
8.17 5.062(—7) 4.376(-7) 4.24(-8) 2.15(—8) 4.30(-9) 
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Table 2A.1. (continued) 


Fluence-to-Kerma Factor (ergs/g)/(neutron/cm?) 


Energy 

(MeV) Total Hydrogen Oxygen Carbon Nitrogen Phosphorus Sodium Chlorine 
8.02 4.971(-7) 4.368(-—-7) 2.94(-8) 2.63(-8) 4.74(-9) 
7.87 5.038(—7) 4,345(—7) 3.81(-8) 2.62(—8) 4.38(—9) 
7.72 5.060(—7) 4.319(-7) 3.89(—8) 3.09(—8) 3.88(-9) 
7.57 4.935(-7) 4.292(-7) 3.65(-8) 2.34(-8) 3.84(-9) 
7.42 4.970(-7) 4.266(-7) 4.11(-8) 2.49(-8) 3.98(-9) 
7.27 4.995(-7) 4.239(-7) 5.32(-8) 1.80(—8) 3.87(-9) 
7.12 4.773(-7) 4.211(-7) 4.13(-8) 1.08(-8) 3,.52(-9) 
6.97 4.650(-7) 4.187(-7) 3.13(-8) 1.10(-8) 3.49(-9) 
6.82 4.745(-7) 4.173(-7) 4.14(—8) 1.21(-8) 3.27(-9) 
6.67 4.707(-7) 4.159(—7) 3.94(-8) 1.17(-8) 3.27(-9) 
6.52 4.605(—7) 4.145(-7) 2.97(—8) 1.25(-8) 3.38(-9) 
6.37 4.731(-7) 4.130(-7) 3.89(-8) 1.71(-8) 3.57(—9) 
6.21 4.611(-7) 4.115(-7) 2.46(-8) 2.09(-8) 3.65(—9) 
6.06 4.572(-7) 4.100(-7) 2.77(-8) 1.53(-8) 3.77(-9) 
5.91 4.633(-7) 4.077(—7) 3.73(-8) 1.42(-8) 3.67(-9) 
5.76 4.457(-7) 4.048(-7) 2.32(—8) 1.38(-8) 3.51(-9) 
5.61 4.502(-7) 4.018(—7) 3.08(-8) 1.37(-8) 3.55(-9) 
5.46 4.312(-7) 3.987(-7) 1.47(-8) 1.42(-8) 3.25(—9) 
5.31 4.294(-7) 3.956(—7) 1.70(—-8) 1.30(—8) 3.39(-9) 
5.16 4.449(—7) 3.925(-7) 3.46(—8) 1.33(-8) 4.17(-9) 
5.01 4.344(-7) 3.892(-7) 2.76(—8) 1.32(-8) 4.01(-9) 
4.69 4.180(—7) 3.810(—7) 1.79(-8) 1.50(-8) 3.73(-9) 
4,50 4.116(—7) 3.760(-7) 1.66(--8) 1.42(-8) 4.34(-9) 
4.31 4.181(-7) 3.709(—7) 2.91(-8) 1.33(-8) 4.48(-9) 
4.16 4.245(-7) 3.666(—7) 3.37(-8) 1.89(-8) 4.98(-9) 
4.01 4.102(-7) 3.623(-7) 2.24(—8) 2.05(—8) 4.78(-9) 
3.82 4.121(-7) 3.562(-7) 2.91(-8) 2.26(—8) 3.88(-9) 
3.67 4.099(-7) 3.511(-7) 3.05(-8) 2.42(--8) 3.67(-9) 
3.52 4.060(—7) 3.459(-7) 3.09(-8) 2.45(—8) 4.48(-9) 
3.33 . 3.981(-7) 3.392(-7) 3.46(-8) 2.07(-8) 3.32(-9) 
3.18 3.727(-7) 3.337(-7) 2.06(-8) 1.54(-8) 2.83(-9) 
3.03 3.555(-7) 3.280(-7) 1.53(-8) 9.44(-9) 2.47(-9) 
2.84 3.499(-7) 3.203(—7) 1.14(-8) 1.60(—8) 2.01(-9) 
2.69 3.396(-7) 3.138(--7) 1.10(-8) 1.26(-8) 1.89(—9) 
2.50 3.279(-7) 3.054(—7) 1.01(-8) 1.09(-8) 1.34(—9) 
2.40 3.150(—7) 3.000(—7) 3.40(—9) 1.03(—8) 1.07(—9) 
2.31 3.136(—7) 2.950(—7) 7.36(-—9) 9.85(-9) 1.21(-9) 
2.21 3.119(-7) 2.898(—7) 1.08(-8) 9.51(-9) 1.61(—9) 
2.10 3.057(-7) 2.835(—-7) 1.14(-8) 9.68(~-9) 9.62(—10) 
2.00 2.995(-7) 2.780(-7) 1.13(-8) 9.16(—9) 9.20(-10) 
1.91 2.952(-7) 2.721(-7) 1.31¢—8) 8.82(-9) 9.82(—10) 
1.80 2.883(-7) 2.648(—7) 1.31(-8) 8.56(-—9) 1.72(-9) 
1.70 2.800(—7) 2.585(—7) 1.21(-—8) 8.36(—9) 8.76(—10) 
1.59 2.717(-7) 2.507(-7) 1.18(—8) 8.12(-9) 9.84(—10) 
1.50 2.644(—7) 2.440(-7) 1.15(—8) 7.93(-9) 9.06(—10) 
1.40 2.577(-7) 2.370(—7) 1.13(-8) 7.73(-9) 1.66(—9) 
1.31 2.602(—7) 2.297(—7) 2.23(-8) 7.52(—9) 6.44(—10) 
1.20 2.426(—7) 2.209(-7) 1,39(-8) 7.25(-9), $.11(-10) 
1.10 2.371(-7) 2.121(—7) 1.63(—8) 7.00(-9) 7.45(—10) 
1.00 2.432(-7) 2.037(-7) 3.21(-8) 6.70(-9) 6.33(-—10) 
9.0(-1) 2.122(~7) 1.939(-7) 1.15(-8) 6.35(-9) 3.39(-10) 
8.0(—-1) 1.980(-7) 1.834(-7) 7.95(-9) 5.97(—9) 3.90(—10) 
7.0(-1) 1.857(-7) 1.730(-7) 6.61(-9) 5.56(—9) 4.97(-10) 
6.0(-1) 1.719(-7) 1.609(-—7) 5.61(-9) 5.07(-9) 2.80(—10) 
5.0(-1) 1.592(-7) 1.483(-7) 5.74(-9) 4,52(-9) §.32(-10) 
4.5(-1) 1.639(-7) 1.403(-7) 1.90¢-8) 4.20(-9) 2.90(—10) 


Energy 
(MeV) 


4.0(-1) 
3.5(-1) 
3.0(-1) 
2.5(-1) 
2.0(-)) 
1.5(-1) 
1.0(-1) 
7.0(—2) 
5.0(—2) 
3.0(—2) 
2.0(—2) 
1.5(-2) 
1.0(-2) 
7.0(—3) 
5.0(-3) 
4.03(~-3) 
3.00(—3) 
2.01(—3) 
1.50(—3) 
1.03(-—-3) 
7.68(—4) 
5,13(—-4) 
4.02(—4) 
3.10(—4) 
2.01(—4) 
1.42(—4) 
1.006(—-4) 
7.1(-5) 
5.03(—5) 
4.11(-5) 
3.43(—5) 
2.74(—S) 
2.05(—S) 
1.60(—5) 
1.26(—5) 
8.0(—6) 
4.0(—6) 
2.0(—-6) 
1.0(--6) 
5.0(—7) 
2.5(-7) 
1.7(-7) 
7.2(-8) 
4.5(—8) 
2.23(-8) 


Total 


1.486(—7) 
1.339(-7) 
1,223(-7) 
1.104(—7) 
9.649(—8) 
8.165(—8) 
6.335(-8) 
4.973(-8) 
3.830(-8) 
2.560(—8) 
1.807(-8) 
1.361(—8) 
9.550(-9) 
6.776(—9) 
4.880(-9) 
3.951(-9) 
2.964(—9) 
2.002(—9) 
1.495(-—9) 
1.030(—9) 
7.774(-10) 
5.299(—10) 
4.244(-10) 
3.347(—10) 
2.299(—10) 
1.764(—10) 
1.417(-10) 
1.192(—10) 
1.071¢—10) 
1.036(—10) 
1.025(—10) 
1.033(—10) 
1.077(-10) 
1.143(—10) 
1.230(—10) 
1.454(—10) 
1.965(—10) 
2.723(—-10) 
3.811¢—10) 
5.354(—10) 
7.529(-10) 
1.144(-9) 
1.398(-9) 
1.766(—9) 
2.489(—9) 


2From ref. 43. 
>The standard-man model also includes the elements calcium, sulfur, potassium, and magnesium but their contributions to the 
kerma are negligible. 
Read: 7.079 X 1077, etc. 


Hydrogen 


1.318(—7) 
1.236(—7) 
1.142(—7) 
1.037(—7) 
9,116(—8) 
7.158(—8) 
6.056(—8) 
4.774(—8) 
3.705(—8) 
2.472(—8) 
1.747(-8) 
1.316(—8) 
9.287(—9) 
6.558(—-9) 
4.832(—9) 
3.820(—9) 
2.855(—9) 
1.928(—9) 
1.436(—9) 
9.832(—10) 
7.368(—10) 
4.925(—10) 
3.867(-10) 
2.986(—10) 
1.939(-10) 
1.372(—10) 
9.753(-11) 
6.920(—11) 
4.94(-11) 
4.07(-11) 
3.43(-11) 
2.78(-11) 
2.13(-11) 
1.71(-11) 
1.40(-11) 
1.01(-11) 
7.20(—12) 
6.63(—12) 
7.59(-12) 
9.87(-12) 
1.35(-11) 
2.03(—11) 
2.48(—11) 
3.14(-1) 
4.18(-11) 
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Table 2A.1. (continued) 


Fluence-to-Kerma Factor (ergs/g)/(neutron/cm?) 


Oxygen 


1.18(—8) 
6.52(—9) 
4.73(-9) 
3.74(-9) 
2.89(—9) 
2.16(—9) 
1.45(-9) 
1.02(-9) 
7.30(—9) 
4.43(-9) 
2.98(—9) 
2.18(—9) 
1.47(-9) 
1.03(-9) 
7.53(-9) 
5.94(-11) 
4.43(—1)) 
2.98(-11) 
2,.22(-11) 
1.52(—11) 
1.14(-11) 
7.60(-12) 
5.96(—12) 
4.60(-12) 
2.98(-12) 
2,.11(-12) 
1.49(-12) 
1.05(--12) 
7.46(-13) 
6.09(—13) 
5.09(—13) 
4.07(-13) 
3.05(—13) 
2.38(—13) 
1,86(—13) 
1.19(-13) 
5.94(-14) 
2.97(-14) 
1.49(-14) 


Carbon 


3.89(—9) 

3.51(-9) 

3.12(—9) 

2.70(—9) 

2.23(—9) 

1,75(-9) 

1.22(-9) 

8.76(—10) 
6.37(—10) 
3.92(—10) 
2.66(—10) 
1.95(-10) 
1.34(—10) 
9.24(-11) 
6.74(-11) 
5.32(-11) 
3.96(—11) 
2.67(—11) 
1,99(-11) 
1,36(-11) 
1.02(-—11) 
6.80(—12) 
$.34(-12) 
4,12(-12) 
2.67(—12) 
1.88(-12) 
1.33(—12) 
9.42(-13) 
6.67(—13) 
$.45(-13) 
4.56(—13) 
3.64(—13) 
2.73(—13) 
2.13(-13) 
1.67(—13) 
1.06(—13) 
$.32(-14) 
2.66(—14) 
1,33(-14) 


Nitrogen 


2.96(—10) 
2.46(—10) 
2.13(—10) 
1.93(—10) 
1.82(—10) 
1.47(-10) 
1.08(—10) 
9,36(-11) 
7.62(-11) 
4.9111) 
3.73(—11) 
3.14(-11) 
2.29(--11) 
1.87(-11) 
1.56(—11) 
1.44(-11) 
1.32(-11) 
1.27(-11) 
1.30(-11) 
1.40(-11) 
1.52(-11) 
1.77(-11) 
1.95(-11) 
2.19(-11) 
2.67(-11) 
3.15(-11) 
3.73(-11) 
4.42(-11) 
5.237(-11) 
5.786(—11) 
6.323(-11) 
7.065(—11) 
8.155(—11) 
9,226(-11) 
1.041(-10) 
1.302(—10) 
1.836(—10) 
2.591(—10) 
3,656(—10) 
5.159(-10) 
7.280(—10) 
1.109(—9) 
1,356(—9) 
1.715(-9) 
2.424(-9) 


Phosphorus 


2.88(—12) 
3.02(-12) 
3.12(-12) 
3.20(—12) 
3.25(-12) 
3.27(-12) 
3.30(—12) 
3.31(-12) 
3.32(-12) 
3.33(-12) 
3.33(-12) 
3.34(—12) 
3.34(-12) 
3.34(—12) 
3.34(—12) 
3.34(-12) 
3.35(—12) 
3.35(-12) 
3.35(-12) 
3.35(-12) 
3.35(-12) 
3.35(—12) 
3.35(—12) 
3.35(—12) 
3.35(—12) 
3.35(-12) 
3.35(—12) 
3.35(-12) 
3.35(-12) 
3.35(-12) 


Sodium 


3.50(—13) 
4.11(-13) 
4.54(-13) 
4.99(-13) 
5.58(—13) 
6.44(—13) 
7.34(—13) 
8.34(—13) 
1.04(-12) 
1.46(-—-12) 
2.09(—12) 
2.97(-12) 
4,15(-12) 
5.89(—12) 
9.00(-12) 
1.10¢—11) 
1.39¢(-11) 
1.83(—11) 


Chlorine 


2.72(—13) 
1.84(-12) 
3.28(—12) 
1.97(-12) 
5§.91(-14) 
1.12(-13) 
3.87(—13) 
9.84(-14) 
1.10(-13) 
1.31(-13) 
1.57(-13) 
1.95(-13) 
2.59(—13) 
3.14(-13) 
3.79(-13) 
5.08(—13) 
7.24(-13) 
1.03(—12) 
1.45(-12) 
2.07(—12) 
2.12(—12) 
2.12(-12) 
2.12(-12) 
2.12(-12) 
2.12(-12) 


© 
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Table 2A.2. Neutron Fluence-to-Kerma Factors fomSpecific Body Components? 


Fluence-to-Kerma Factors (ergs/g)/(neutron/ cm?) 


Energy 
(MeV) Lung Muscle Bone Brain Red Marrow 
15.00 7.04(-7)8 6.97(-7) 5.60{—7) 6.98(—7) 7.01(—7) 
14.54 6.85(-7) 6.78(—7) 5.41(—7) 6.78(—7) 6.81(-7) 
13.94 6.52(-7) 6.46(-7) 5.17(~7) 6.46(—7) 6.48(-7) 
13.48 6.14(-7) 6.08(-7) 4.94(-7) 6.10(—7) 6.11(—7) 
13.03 5.96(-7) 5.90(-7) 4.81(-7) 5.92(-7) 5.93(-7) 
12.43 5.96(—7) 5.90(—7) 4,.73(-7) 5.91(-7) 5.92(—7) 
11.98 6.00(—7) 5.94(—7) 4.72(-7) 5.94(-7) 5.96(-7) 
11.52 5.96(-7) 5.90(—7) 4.63(—7) 5.88(—7) 5.91(-7) 
11.07 5.85(—7) 5.79(-7) 4.55(-7) 5.78(—7) 5.80(-7) 
10.47 5.61(—7) 5.54(-7) 4.34(-7) 5.53(—7) 5.55(-7) 
10.00 5.53(—7) 5.46(-7) 4.25(-7) 5.44(-7) 5.47(-7) 
9.53 5.38(-7) 5.31(-7) 4.15(-7) 5.30(-7) 5.32(-7) 
9.00 5.21(-7) 5.14(-7) 4.01(-7) 5.13(-7) 5.15(-7) 
8.55 5.05(—7) 4.99(-7) 3.90(—7) 4.98(—7) 5.00(-7) 
8.02 4.96(-7) 4.90(-7) 3.91(-7) 4.92(-7) 4.92(-7) 
TAI 4.94(-7) 4.88(-7) 3.85(—7) 4.89(—7) 4.98(—7) 
6.97 4.73(-7) 4.67(-7) 3.60(—7) 4.66(-—7) 4.67(—7) 
6.52 4.67(—7) 4.61(-7) 3.57(-7) 4.60(-7) 4.61(—7) 
6.00 4.58(-7) 4.53(—7) 3.52(—7) 4.52(-7) 4.53(—7) 
5.54 4.41(-7) 4.36(-7) 3.40(-7) 4.36(—7) 4.37(-7) 
5.00 4.39(-7) 4.34(-7) 3.36(—7) 4.33(—7) 4.35(-7) 
4.50 4.13(—7) 4.08(—7) 3.20(-7) 4.08(—7) 4.09(—7) 
4.00 4.09(-7) 4.04(-7) 3.19{—7) 4.05(-7) 4.06(—7) 
3.52 4.05(-7) 4.00(-7) 3.15(-7) 4.01(-7) 4.02(-7) 
3.00 3.59(-7) 3.54(-7) 2.77(-7) 3.55(—7) 3.55(-7) 
2.50 3.28(-7) 3.24(-7) 2.51(-7) 3.25(-7) 3.25(-7) 
2.00 3.01(-7) 2.91(-7) 2.29(-7) 2.98(—7) 2.98(—7) 
1.80 2.90(-7) 2.87(-7) 2.20(-7) 2.87(-7) 2.87(-7) 
1.60 2.74(-7) 2.71(-7) 2.08(—7) 2.71(-7) 2.71(-7) 
1.40 2.59(-7) 2.56(—7) 1.97(-7) 2.56(—7) 2.57(—7) 
1.20 2.45(—7) 2.42(-7) 1.84(-7) 2.42(—7) 2.42(-7) 
1.00 2.51(-7) 2.47(-7) 1.83(-7) 2.46(-7) 2.47(-7) 
8.05(-1) 2.00(—7) 1.97(-7) 1.51(—7) 1.98(—7) 1.98(-7) 
6.0(-1) 1.73(-7) 1.71(-7) 1.31(-7) 1.71(-7) 1.71(—7) 
5.0(—1) 1.60(—7) 1.58(—7) 1.21(-7) 1.58(—7) 1.58(-7) 
4.0(-1) 1.51(-7) 1.48(-7) 1.11(-7) 1.48(—7) 1,50(—7) 
3.0(-—1) 1.23(-7) 1.22(-7) 9.27(-8) 1.22(—7) 1.22(-7) 
2.0(—1) 9.73(—8) 9.67(—8) 7.33(—8) 9.61(-8) 9.63(—8) 
1.0(-1) 6.39(-8) 6.30(—8) 4.80(-8) 6.31(-8) 6.31(-8) 
5.03(—2) 3.88(-8) 3.83(—8) 2.92(--8) 3.84(-—8) 3.84(-8) 
2.04(—2) 1.82(-8) 1.80(—8) 1.38(—8) 1.80(—8) 1.80(-8) 
1.02(—2) 9.69(-9) 9.57(-9) 7.54(-9) 9.57(-9) 9.57(-9) 
4.97(—3) 4.90(-9) 4.84(-9) 4.02(-9) 4.84(-9) 4.84(-9) 
3.00(—3) 2.99(-9) 2.95(-9) 2.61(-9) 2.96(-9) 2.94(-9) 
2.014(—3) 2.02(—9) 2.00(—9) 1.90(—9) 2.00(—9) 1.99(~—9) 
1.025(-3) 1.04(--9) 1.03(—9) 1.18(—9) 1.03(-9) 1.025(-9) 
4.94(—4) 5.16(-10) 5.11¢-10) 7.94(-10) §.12(-10) 5.08(—10) 
2.01(—4) 2.30(—10) 2.29(—10) 5.94(-10) 2.27(-10) 2.29(-10) 
1.00(—4) 1.39(-—-10) 1.39(-10) 5.38(—10) 1.36(—10) 1.41(-10) 
5.03(—5) 1.03(—10) 1.04(-10) 5.26(—10) 9.59(-11) 1.09(-10) 
3.43(—5) 9.72(-11) 9.78(-11) 5.32(—10) 8.76(-11) 1.05(-10) 
2.05(—5) 1.01(—10) 1.01(—10) 5.52(—10) 8.65(-11) 1,12(-10) 
1.03(—5) 1.22(-10) 1,22(-10) 6.00(-10) 9.93(-11) 1.40(-10) 
5.14(-6) 1.61(-10) 1.61(-—10) 6.75(-10) 1.27(-10) 1,88(—10) 
2.00(—6) 2.49(-10) 2.47(-10) 8.33(-10) 1,90(-10) 2.93(—10) 
1.00(—6) 3.48(—10) 3.45(-10) 1.01(-9) 2.62(—10) 4.12(-10) 
5.0(~7) 4.89(-10) 4.83(-10) 1.26(-9) 3.65(—10) 5.81(-10) 
2.5(~7) 6.87(—10) 6.79(—10) 1.61(-9) 5.09(—10) 8.18(-10) 
1.07(-7) 1.04(-9) 1.03(-9) 2.24(-9) 7.68(—10) 1.25(-9) 
4.5(~8) 1.61(-9) . 1.59{-9) 3.24(-9) 1,18(-9) 1.93(-9) 
2.23(-8) 2.27(-9) 2.24(-9) 4.42(-9) 1.66(—9) 2.72(—9) 


4From ref. 43. 


bRead: 7.04 X 107, etc. 


Table 2A.3. Unweighted Neutron Fluence-to-Kerma Factors for the Most Common Elements in Man” 


Energy Unweighted Fluence-to-Kerma Factor [(ergs/g)/(neutron/cm?)] /{atoms/g] X 1024 
(MeV) Hydrogen Oxygen Carbon Nitrogen Calcium Phosphorus Sulfur Potassium Sodium Chlorine Magnesium 
15.00 7.75(-6)® _ 6.67(—6) 6.56(—6) 7.41(-6) 5.02(-7) 7.92(-6) 9.67(-7) 1.35(—5) 1.38(-6) 4.75(-6) 1.11(—5) 
14.84 7.74(—6) 6.48(-6) 6.25(-6) 7.32(-6) 5.02(-7) 7.84(—6) 9.51(-7) 1.34(-5) 1.37(-6) 4.70(—6) 1.11(-5) 
14.54 7.73(-6) 6.20(-6) 5.73(-6) 7.11(-6) 4.99(-7) 7.69(—6) 9.17(—7) 1.30(—S) 1.35(-6) 4.61(-6) 1.11(~5) 
14.24 7.72(-6) 5.72(-6) 5.37(-6) 6.95(-6) 4.94(-7) 7.55(—6) 8.83(-7) 1.28(-5) 1.33(-6) 4.52(-6) 1.11(-5) 
13.94 7.71(-6) 5.23(-6) 5.11(-6) 6.75(--6) 4.91(-7) 7.41(-6) 8.46(—7) 1.25(—5) 1.31(-6) 4.36(-6) 1.10(-5) 
13.63 7.70(-6) 4.38(—-6) 4.93(-6) 6.41(-6) 4.95(-7) 7.26(-6) 8.01(-7) 1.23(-5) 1.29(-6) 4.19(-6) 1.08(—5) 
13.33 7.69(-6) 3.82(-6) 4.75(-6) 6.21(-6) 4.99(-7) 7.10(—6) 7.57(-7) 1.20(—S) 1.28(-6)  4.03(—6) 1.07(-5) 
13.03 7.68(-6) 3.48(—6) 4.52(-6) 6.12(-6) 4.97(-17) 6.94(—6) 7.20(—7) 1.18(—5) 1,.25(—6) 3.87(-6) 1.04(-5) 
12.73 7.67(—6) 3.55(-6) 4.12(-6) 5.97(-6) 4.93(-7) 6.79(—6) 6.84(-7) 1.15(—5) 1.21(-6) 3.73(—6) 1.01(-5) 
12.43 -7.66(—6) 3.67(—6) 3.83(-6) 5.84(—6) 4.96(-7) 6.64(—6) 6.53(-7) 1.13(-5) 1.18(-6) 3.61(—6) 9.74(-6) 
12.13 7.64(—6) 3.74(—6) 3.69(-6) 5.65(—6) 5.00(-7) 6.49(-6) 6.23(-7) 1.10(-5) 1.15(-6) 3.50(-6) 9.34(-6) 
11.83 7.63(—6) 4.11(—6) 3.33(-6) 5.43(-6) 4.93(—7) 6.30(—6) 5.96(-7) 1.08(—5) 1.13(—6) 3.40(—6) 8.89(—6) 
11.52 7.61(-6) 3.94(-6) 3.02(-6) 5.19(-6) 4.85(-7) 6.08(—6) 5.71(-7) 1.06(—5) 1.11(-6) 3.30(-6) 8.44(-6) 
11.22 7.60(—6) 3.75({-6) 2.98(-6) 4.99(-6) 4.74(-7) 5.87(—6) 5.50(-7) 1.04(—5) 1.09(-6) 3.21(-6) 7.90(-6) 
10.92 7.58(-6) 3.32(-6) 2.85(—-6) 4.71(-6) 4.63(-7) 5.66(—6) 5.30(-7) 1.01(—5) 1.07(—6) 3.13(-6) 7.38(-6) 
10.62 7.56(—6) 2.99(—6) 2.41(-6) 4.69(-6) 4.63(-7) 5.46(—6) $.12(-7) 9.86(—6) 1.05(—6) 3.04(—6) 6.85(-6) 
10.32 7.55(—6) 3.01(-6) 1.88(—6) 4.59(—6) 4.62(-7) 5.26(—6) 4.95(-7) 9.61(—6) 1.02(-6) 2.97(—6) 6.32(—-6) 
10.02 7.53(—6) 2.86(—6) 1.76(-6) 4.42(-—6) 4.54(-7) 5.06(—6) 4.81(-7) 9.37(-6) 9.91(-7) 2.91(—6) 5.76(—-6) a 
9.72 7.49(-6) 2.65(—6) 1.81(—6) 4.24(-6) 4.44(-7) 4.88(—6) 4.69(—7) 9.13(—6) 9.58(—7) 2.84(—6) 5.22(-6) oO 
9.60 7.48(—6) 2.52(—6) 1.80(—6) 4.17(-6) 4.41(-7) 4.81(-6) 4.64(-7) 9.03(—6) 9.44(—7) 2.81(—-6) 5.02(-6) 
9.45 7.45(-6) 2.39(-6) 1.82(—6) 3.99(-6) 4.36(-7) 4.72(-6) 4.59(-7) 8.91(—6) 9.23(-7) 2.78(—6) 4.76(—6) 
9.30 7.43(—6) 2.21(-6) 1.85(—6) 3.94(-6) 4.31(-7) 4.62(-6) 4.57(-7) 8.78(—6) 9.00(—7) 2.75(—6) 4,.51(-6) 
9.15 7.41(-6) 2.10(—6) 1.78(—6) 3.89(-6) 4.26(-7) 4.53(—6) 4.57(-7) 8.66(—6) 8.77(—7) 2.72(—6) 4.28(-6) 
9.00 7.39(-6) 2.12(-6) 1.45(-6) 3.77(-6) 4.21(-7) 4.44(-6) 4.57(-7) 8.53(—6) 8.57(—7) 2.69(-6) 4.06(-6) 
8.85 7.38(-6) 2.01(-6) 1.34(-6) 3.72(—6) 4.17(-7) 4.27(—6) 4.56(-7) 8.39(-6) 8.43(—7) 2.66(—6) 3.81(-6) 
8.70 7.37(-6) 1.93(-6) 1.30(-6) 3.69(—6) 4.14(-7) 4.37(-6) 4.57(-7) 8.26(—6) 8.32(-7) 2.63(—6) — 3.56(—6) 
8.62 7.36(-6) 1.78(-6) 1.29(-6) 3.63(-6) . 4.12(-7) 4.34(-6) 4.57(-7) 8.20(—6) 8.28(-7) 2.62(—6) 3.45(—6) 
8.47 7.35(-6) 1.73(—6) 1.30(—6) 3.53(-6) 4.09(-7) 4.05(-6) 4.56(-7) 8.07(-6) 8.30(—6) 2.58(—6) 3.23(—6) 
8.32 7.33(-6) 1.94(-6) 1.38(-6) 3.44(—6) 4.03(-7) 4.07(-6) 4.54(-7) 7.94(-6) 8.27(-7) 2.54(—6) 3.05(—6) 
8.17 7.32(-6) 1.87(-6) 1.78(—6) 3.45(-6) 3.97(-7) 3.58(—6) 4.53(-7) 7.80(-6) 8.23(-7) 2.50(—6) 2.89(-6) 
8.02 7.31(-6) 1.30(-~6) 2.18(—6) 3.80(—6) 3.92(-7) 3.74(-6) 4.51(-7) 7.67(—6) 8.16(—7) 2.46(—6) 2.63(—6) 
7.87 7.27(-6) 1.69(-6) 2.18(—6) 3.51(-6) 3.88(—7) 3.65(—6) 4.49(-7) 7.52(-6) 8.06(—7) 2.42(-6) 2.41(-6) 
7.72 7.22(—6) 1.72(-6) 2.56(—6) 3.11(-6) 3.84(-7) 3.45(—6) 4.48(-7) 7.37(—6) 7.95(-7) 2.37(—6) 2.24(—6) 
7.57 7.18(—6) 1.61(-6) 1.95(-6) 3.08(—6) 3.83(-7) 3.45(—6) 4.47(-7) 7.21(—6) 7.87(-7) 2.33(-6) 2.04(—6) 
7.42 7.14(-6) 1.82(-6) 2.07(-6) 3.19(—6) 3.82(-7) 3.28(—6) 4.47(-7) 7.05(—6) 7.79(-7) 2.29(-6) 1.85(-6) 
7.27 7.09(-6) 2.36(—6) 1.50(—6) 3.10(—6) 3.82(-7) 3.15(-6) 4.46(-7) 6.88(—6) 7.72(-7) 2.24(—6) 1.65(-6) 
7.12 7.05(—6) 1.83(—6) 8.98(-7) 2.83(-6) 3.82(—7) 2.90(—6) 4.46(-7) 6.71(—6) 7.64(—7) 2.20(—6) 1.46(—6) 
6.97 7.00(—6) 1.39(-6) 9.13(-7) 2.80(—6) 3.82(-7) 2.93(-6) 4.46(-7) 6.54(—6) 7.51(-7) 2.16(—6) 1.31(-6) 
6.82 6.98(-6) 1.83(-6) 1.00(—6) 2.62(-6) 3.78(-7) 2.91(-6) 4.43(-7) 6.37(-6) 7.61(-7) 2.11(-6) 1.24(-6) 
6.67 6.96(—6) 1.74(-6) 9.76(-7) 2.62(-6) 3.74(-7) 2.77(-6) 4.39(—7) 6.21(—6) 7.69(—6) 2.07(—6) 1.07(-6) 
6.52 6.93(-6) 1.31(—6) 1.04(—6) 2.71(-6) 3.70(-7) 2.65(—2) 4.37(-7) 6.06(—6) 7.81(—6) 2.02(-6) 9.97(-7) 
6.37 6.91(-6) 1.72(-6) 1.42(-6) 2.86(—6) 3.66(—7) 2.58(—6) 4.34(-7) 5.90(—6) 8.06(—6) 1.98(—6) 8.69(—7) 


Energy 
(MeV) 
6.21 
6.06 
5.91 
5.76 
5.61 
5.46 
5.31 
5.16 
5.01 
4.69 


Hydrogen 


~6.89(—6) 


6.86(—6) 
6.82(—6) 
6.77(—6) 
6.72(-6) 
6.67(—6) 
6.62(—6) 
6.57(-6) 
6.51(-6) 
6.37(-6) 
6.29(-6) 
6.20(—6) 
6.13(-6) 
6.06(—6) 
5.96(—6) 
5.87(—6) 
5.79(-6) 
5.68(—6) 
5.58(-6) 
5.49(-6) 
5.36(—6) 
5.25(-6) 
5.11(-6) 
5.02(—6) 
4.93(-6) 
4.85(—6) 
4.74(-6) 
4.64(-6) 
4.53(—6) 
4.43(—6) 
4.32(-6) 
4.19(—6) 
4.08(—6) 
3.96(—6) 
3.84(-6) 
3.70(—6) 
3.55(—6) 
3.41(-6) 
3.24(—6) 
3.07(-6) 
2.89(—6) 
2.69(—6) 
2.48(—6) 


Oxygen 


1.09(-6) 
1.22(-6) 
1.65(-6) 
1.03(—6) 
1.36(-6) 
6.49(-6) 
7.52(—-6) 
1.53(—6) 
1.22(—6) 
7.91(—7) 
7.36(-7) 
1.29(-6) 
1.49(-6) 
9.91(-7) 
1.29(-6) 
1.35(—6) 
1.37(-6) 
1.53(-6) 
9.11(-7) 
6.77(-7) 
5.04(-7) 
4.88(-7) 
4.45(-7) 
1.50(-—7) 
3.26(—7) 
4.78(-7) 
5.05(—7) 
5.00(—7) 
7.33(-7) 
5.78(-7) 
5.37(-7) 
5.20(-7) 
5.07(-7) 
4.98(-7) 
9.87(-7) 
6.13(-7) 
7.55(-7) 
1.42(-7) 
5.10(-7) 
3.45(-7) 
2.93(-7) 
2.48(—7) 
2.54(-7) 


24 
Unweighted Fluence-to-Kerma Factor [(ergs/g)/(neutron/cm?)] /{atoms/g] X 10 


Carbon 


1.74(—6) 
1.27(—6) 
1.18(—6) 
1.14(—6) 
1.14(-6) 
1.18(—6) 
1.08(--6) 
1.10(-6) 
1.10(-6) 
1.25(-6) 
1.18(—6) 
1.11(—-6) 
1.57(—6) 
1.70(—6) 
1.88(—-6) 
2.01(—6) 
2.03(-6) 
1.72(-6) 
1.27(-6) 
7.84(-7) 
1.33(-6) 
1.05(—6) 
9.01(-7) 
8.52(-~7) 
8.18(—7) 
7.90(—7) 
8.04(—7) 
7.60(—7) 
7.29—-T) 
7.11(-7) 
6.94(—7) 
6.75(—7) 
6.58(—7) 
6.42(—7) 
6.24(—7) 
6.02(—7) 
5.79(-7) 
5.56(—7) 
5.28(—7) 
4.95(-7) 
4.62(—7) 
4.21(-7) 
3.76(-—-7) 


Nitrogen 


2.93(—6) 
3.03(—6) 
2.95(—6) 
2.82(—6) 
2.85(—6) 
2.61(—6) 
2.72(-6) 
3.34(-6) 
3.22(—6) 
2.99(-6) 
3.48(-6) 
3.60(—6) 
4.00(—6) 
3.83(—6) 
3.12(-6) 
2.94(—6) 
3.59(-6) 
2.67(—6) 
2.27(-6) 
1.98(—6) 
1.61(-6) 
1.52(-6) 
1.07(—6) 
8.57(-7) 
9.74(-7) 
1.29(--6) 
7.72(—-7) 
7.38(-7) 
7.90(—-7) 
1.38(—6) 
7.03(—7) 
7.89(-7) 
7.26(—7) 
1.33(—6) 
5.17(-7) 
4.10(-7) 
5.10(-7) 
5.08(-7) 
2.72(-7) 
3.08(—7) 
3.99(-7) 
2.24(—7) 
4.26(-7) 


Calcium 


3.63(—-7) 
3.62(—7) 
3.65(—7) 
3.71(-7) 
3.71(-7) 
3.83(-7) 
3.94(-7) 
4.06(—7) 
4.23(-7) 
4.34(-7) 
4.39(-7) 
4.47(-7) 
4.56(—7) 
4.64(-7) 
4.65(-7) 
4.69(—7) 
4.74(—7) 
5.00(—7) 
5.03(—7) 
4.91(-7) 
4.76(-7) 
4.64(—7) 
4.35(-7) 
4.12(-7) 
3.92(-7) 
3.60(—7) 
3.44(-7) 
3.30(-7) 
3.10(—7) 
2.98(-7) 
2.6%(-7) 
2.42(-7) 
2.23(-7) 
2.10(-7) 
2.02(-7) 
1.96(-7) 
1.87(-7) 
1.85(—7) 
1.86(—7) 
1.30(-7) 
1.19(-7) 
1.52(-7) 
7.55(~7) 


Table 2A.3. (continued) 


Phosphorus 


2.50(—6) 
2.39(-6) 
2.34(-6) 
2.26(—6) 
2.19(-6) 
2.11(-6) 
2.04(—6) 
1.98(-6) 
1.92(-6) 
1.87(-6) 
1.64(—6) 
1.61(-6) 
1.40(—6) 
1.40(—6) 
1.38(—6) 
1.37(—6) 
1.06(—6) 
9.83(—7) 
1.09(—6) 
9.22(-7) 
8.27(-7) 
8.06(—7) 
6.26(-7) 
5.57(-7) 
5.33(-7) 
5.05(-7) 
4.57(-7) 
4.23(-7) 
3.79(-7) 
3.771(-7) 
3.58(-7) 
3.43(-7) 
3.05(—7) 
2.85(—7) 
2.66(—7) 
2.49(—7) 
2.16(—7) 
1.76(-—7) 
1.55(-7) 
1.16(—-7) 
1.07(-7) 
1.51(-7) 
1.05(-7) 


Sulfur 


4,35(-7) 
4.38(-7) 
4.29(—7) 
4.12(-7) 
4.07(-7) 
4.12(-7) 
4.14(-7) 
4,15(-7) 
4.18(-7) 
4.29(~7) 
4,25(-7) 
4.21(-7) 
4.21(-7) 
4.22(—7) 
4.20(-7) 
3.95(-7) 
3.71(-7) 
3.71(-7) 
4.10(—7) 
4.35(-7) 
4.43(—7) 
4.30(—7) 
3.70(-7) 
3.47(-7) 
3.23(-7) 
S127} 
2.96(~7) 
2.83(-7) 
2.59(-7) 
2.45(—7) 
2.44(—7) 
2.30(-7) 
2.09(-7) 
2.02(-7) 
2.41(-7) 
1.80(-7) 
1.46(—7) 
1.58(-7) 
1.30(-7) 
1.34(-7) 
1.27(-7) 
7.73(-8) 
7.24(-8) 


Potassium 


5.74(—6) 
5.51(-6) 
5.40(-6) 
5.23(—6) 
5.05(—6) 
4.91(-6) 
4.76(—6) 
4.60(—6) 
4.43(—6) 
4.04(-6) 
3.84(—6) 
3.57(-6) 
3.35(-6) 
3.15(-6) 
2.78(-6) 
2.49(-6) 
2.16(—6) 
1.80(—6) 
1.50(-6) 
1.29(-6) 
1.13(-6) 
1.01(-6) 
8.51(—7) 
7.61(-7) 
7.03(-7) 
6.42(-7) 
5.71(-7) 
5.52(-7) 
4.98(—7) 
4.32(-7) 
3.99(-7) 
3.43(-7) 
3.24(—7) 
2.64(—7) 
2.43(-7) 
2.07(-7) 
1.80(—7) 
1.54(~7) 
1.20(-7) 
1.01(-7) 
8.99(—8) 
7.85(—8) 
7.30(—8) 


Sodium 


8.19(-6) 
7.72(—6) 
7.42(-6) 
7.45(-6) 
7.58(-6) 
7.72(-6) 
7.84(-6) 
7.74(-6) 
6.98(—6) 
7,.22(—6) 
6.35(—6) 
7.53(-6) 
6.55(—6) 
6.88(-6) 
6.21(—6) 
7.04(—6) 
6.81(—-6) 
6.19(-6) 
5.62(-6) 
5.13(-6) 
6.13(—-6) 
5.95(—6) 
6.67(—6) 
6.70(—6) 
6.00(—6) 
5.89(-6) 
6.45(-6) 
5.41(—6) 
3.88(-6) 
5.33(—6) 
4.59(—6) 
4.42(-6) 
4.01(—6) 
4.73(—6) 
4.40(-6) 
3.85(-7) 
4.62(-7) 
3.13(—7) 
3.98(—7) 
4,58(-7) 
6.13(-7) 
3.54(-7) 
1.36(—7) 


Chlorine 


1.92(-6) 
1.87(—6) 
1.82(—6) 
1.77(-6) 
1.71(—6) 
1.65(—6) 
1.59(-6) 
1.53(—6) 
1.47(~6) 
1.34(—6) 
1.26(—6) 
1.17(-6) 
1.10(—6) 
1.03(—6) 
9.42(-7) 
8.67(—7) 
7.95(-7) 
7.26(-7) 
6.74(-7) 
6.40(—7) 
6.01(-7) 
5.79(-7) 
5.51(-7) 
5.36(-7) 
5.21(-7) 
5.04(—7) 
4.84(-7) 
4.66(—7) 
4.43(-7) 
4.23(—7) 
4.02(-7) 
3.75(-7) 
3.51(-7) 
3.25(-7) 
2.99(-7) 
2.70(-7) 
2.42(-7) 
2.16(-7) 
1.88(—7) 
1.62(-7) 
1.34(-7) 
1.16(—7) 
9.81(-8) 


Magnesium 


7.65(—7) 
6.84(—7) 
7.93(—7) 
8.22(-7) 
8.18(—7) 
7.85(—7) 
7.44(—7) 
7.87(-7) 
7.68(—7) 
7.91(-7) 
7.20(-—7) 
7.10(—7) 
7.29(-7) 
5.85(—7) 
5.26(-7) 
5.83(-7) 
6.59(—7) 
5.58(—7) 
4.19(-7) 
6.20(-7) 
5.72(-7) 
6.12(-7) 
4.67(—7) 
4.42(-7) 
5.05(—7) 
5.65(-7) 
5.57(—7) 
5.19(—7) 
4.47(-7) 
4.13(—7) 
4.90(—7) 
3.15(—7) 
3.10(—7) 
3.66(—7) 
4.69(-7) 
3.02(-7) 
1.72(—7) 
1.99(—7) 
2.37(-7) 
2.48(—7) 
2.54(—7) 
1.97(-7) 
1.89(—7) 


L9 


___ Table 2A.3. (continued) 
Snecma 
Energy Unweighted Fluence-to-Kerma Factor {(ergs/g)/(neutron/cm?)] /{atoms/g} X 1024 
(MeV) Hydrogen Oxygen Carbon Nitrogen Calcium Phosphorus Sulfur Potassium Sodium Chlorine Magnesium 


4.5(-1) 2.35(—6) 8.41(—7) 3.49(-7) 2.32(-7) —‘1.51(—7) 1.12(-7) 7.19(—8) 7.05(—8) 2.18(-7) 9.06(-8) _3.71(—7) 
3.98(-1) —2.20(—6) 4.80(-7) 3.20(—7) DCT): 477) 8.09(—8) 6.63(—8) 6.00(—8) 2.11(-7) 9.05(-8) —-2.13(—7) 
3.51(-1) _—-2.07(—6) 2.88(—7) 2.92(-7) 1.97(-7)  1.14(—7) 9.92(—8) 5.64(—8) 5.88(—8) 1.45(-7) —-8.28(-7) —-2.19(~-7) 
3.00(-1) —-1.91(—6) 2.09%(—7) 2.59(-7) 1.71(-7)  5.84(—7) 6.22(-8) 4.84(—8) 7.89(-8) 1.66(-7)  5.31(-8)  3.34(-7) 
2.51(-1) —_‘1.73(—6) 1.65(—7) 2.25(—7) 1.55(-7) —7.38(—7) 5.80(—8) 4.10(—-8) 3.62(-8) 1.51(-7) — 3.86(-8) —_2.90(—7) 
2.00(-1) —1.53(—6) 1.28(—7) 1.89(-7) 1.46(-7) —-1.50(—7) 5.62(—8) 9.06(—8) 2.16(-8) 1.76(-7) 3.9% -8) —_1.33(—7) 
1.50(-1) —_1.30(—6) 9.55(—8) 1.45(-7) 1.18(-7)  6.74(—9) 5.06(—8) 5.18(—8) 4.31(—8) 6.38(-8)  2.82(-8) _7.57(—8) 
1.01(-1) —-1.01(—6) 6.41(—8) 1.01(—7) 8.64(-8) — 1.02(-8) 3.40(—8) 7.10(—8) 3.27(-8) 4.31(-8) —-1.84(-8) —-9.51(—8) 
7.03(-2) —-7.99(—7) 4.50(—8) 7.21(-8) 7.51(-8) _7.04(—9) 2.38(-8) 2.80(—9) 5.91(-8) 3.70(-8)  1.69(-8)  5.44(—8) 
5.03(-2)  6.20(—7) 3.23(—8) 5.29(-8) 6.12(-8)  5.05(—9) 1.93(—8) 2.39(—9) 2.56(—9) 3.51(-8)  1.35(-8) __2.25(—8) 
3.05(-2)  4.13(—7) 1.96(—8) 3.26(—8) 3.99(-8)  3.07(—9) 1.61(—8) 1.98(—9) 4.09(-9) 1.78(-8) 9.12(-9) —_1.40(-—-8) 
2.04(-2)  2.92(-7) 1.32(-8) 2.21(-8) 2.99(-8)  2.06(—9) 6.21(—8) 1.71(—9) 2.72(-9) 1.22(-8) 6.93(-9) _—-1.92¢-8) 
1.50(-2) —_2.20(-7) 9.64(—9) 1,62(-8) 2.52(-8)  1.51(—9) 8.75(—8) 1.48(-9) 2.06(—9) 9.65(-9)  5.83(-9) 6.959) 
1.02(-2) —_1.55(—7) 6.62(—9) 1.11(—8) 1.84(-8) 1.039) 1.09(—7) 1.15(—9) 2.27(-9) 8.45(-9) 2.94(-9) —-4.48(—9) 
7.03(-3)  1.10(—7) 4.51(—-9) 7.61(-9) 1.50(-8) 7.10(-10) —-1.24(—7) 8.04(-10) 7.21(-10) —-7.83(-9) —-2.96(-9) __—-3.06(—-9) 
5.11(-3)  8.08(—8) 3.33(—9) 5.60(—9) 1.26(-8) 5.17(-10) —-1.33(—7) 5.89(-10)  5.69(-10) 9.44(-9) —-3.89(-9) _—-2.22(-9) 
4.03(-3)  6.39(—8) 2.63(—9) 4.42(--9) 1.15(-6) 4.07(-10)—-1.38(--7) 4.67(-10) 5.54(-10) —-1.76(-8) -3.35(-9)_—:1.74(-9) 
3.00(-3)  4.78(—8) 1.96(—9) 3.29(-9) 1.06(-8)  3.03(-10) _—‘1.43(—7) 3.51(-10) 2.46(-10) —-1.71(-7) —-1.92(-9) _——1.29(--9) 
2.01(-3)  3.23(-8) 1.32(—9) 2.21(-9) 1.02(-8) 2.04(-10) —-1.47(-7) 2.40(-10)  2.09(-10) 1.62(-8) —1.65(-9) —_8.66(--10) 
1.50(-3)  2.40(—8) 9.82(-10)  1.65(—9) 1.04(-8) 1.52(-10) ——-1.50(—7) 1.81(-10) 1.68(-10) 6.77(-9)_—-1.97(—9)__-6.43(-10) 
1.03(-3) —-1.64(—8) 6.72(-10) —-1.13(—9) 1.12(-8) 1.04(-10) —-1.52(-7) 1.28(-10) 1.23(-10) 4.29(-9)_ 3.54(-9) —-4.40(~10) 
7.68(-4) —-1.23(-8) 5.04(-10) 8.46(-10) —1.22(-8)—7.77(-11)_——«'1.53(—7) 9.83(-11) 9.78(-11) -3.67(-9) —-8.01(-9)__—-3.30(—-10) 
5.13(-4)  8.24(—9) 3.36(-10)  $.65(-10) —-1.42(-8)-5.19(-11)_—«1.54(—7) 6.98(-11) -7.30(-11) —-3.43(-9) —-5.42(-8) _—-2.21(—10) 
3.10(-4) — $.00(—-9) 2.04(-10) 3.42(-10)—:1.76(-8) —3.14(-11) —-1..55(~7) 4.80(-11)  5.48(-11) —3.64(-9)_—-5.80(-8) _—-1.34(~10) 
1.42(-4) —-2.30(—-9) 9.32(-11) —-1.86(-10) —-2.53(-8)-1.44(-11)—«'1.56(-7) 3.21(-11) 4.74(-11)—4.70(-9)—-3.30(-9)_—6.22(-11) 
1.01(-4) —_1.63(—9) 6.60(-11) —- 1.11(-10) ~—-2.99(-8) —1.02(-11)_—‘1.56(—7) 3.00(-11)  4.90(-11) 5.63(-9)—1.14(-8) 449-11) 
7.10(-5)  1.16(—9) 4.66(-11) —-7.82(-11) —3.54(—8) —7.27(-12)_——‘1..56(—7) 2.93(-11) 5.38(-11) 6.68(-9) 2.899) 3.2711) 
5.03(-5)  8.27(-10) —3.30(-11)—5.54(-11) —4.20(-8) —-8.43(-12) __—‘1.56(—7) 2.98(-11)  6.0%-11) _7.85(-9) —-3.23(-9)_-—-2.44(-11) 
3.88(-5)  6.44(-10) .2.55(-11)—4.28(-11) 4.778) —-6.87(-12) _1.56(—-7) 3.16(-11) 6.82(-11) 8.93(-9) 4.11(-9)_—«1.99(-11) 
3.43(-5)  5.73(-10) -2.25(-11)_—3.79(-11) — 5.07(-8) —-5.61(—12)_—:1.56(—7) 3.28(-11) -7.23(-11) —-9.52(-9)4.61(-9)——-1.82(—11) 
2.74(-5) 4.65(-10) —-1.80(-11)—-3.03(-11)-—5.67(-8)—-3.76(-12)_—'1.56(—7) 3.45(-11) 8.05(-11) —-:1.07(-8) —5.72(-9)—-1.57(-11) 
2.05(-5)  3.57(-10)1.35(-11)_—-2.27(-11)—6.54(-8) —-2.23(-12) 156-7) 3.80(-11) 9.29(-11) —-:1.23(-8) -7.62(-9) ——-1.34(—11) 
1.60(-5) _2.87(-10)—1.05(-11)—-1.77(-11)—7.40(-8) ~—'1.62(--12) __—‘1.56(—7) 4.14(-11) 1.05(-10) -1.40(-8) —-9.23(-9)_——1.20(—-11) 
1.26(-5) —-2.35(-10) —8.25(-12)_1.39(-11)—-8.34(-8)—1.27(-12) _——‘1.56(—7) 4.54(-11) 1.19(-10)—-1.59(-8) —1.11(-8) ——«1.12(-11) 
8.01(-6) 1.69(-10)  $.26(-12) —8.83(-12) —-1.04(-7) —-8.11(-13) __—*1,56(—7) 5.58(-11) 1.51(-10) 1.99(-8) —-1.49(-8) _—*1.08(—11) 
4.00(-6) —1.20(-10)-2.63(-12) 4.42(-12) —1.47(-7)—4.05(-13) __—1.56(—7) 7.82(-11) 2.13(-10) —-2.79(-8)—.2.13(-8)—1.21(-11) 
2.00(—6) —-1.11(—10)—1.31(-12) —-2.21(-12) -2.08(-7)_—2.03(—13)__—*1.56(—7) 1.11(-10) 3.04(-10)-3.99(-8)_ 3.03(-8)_~—«1..56(—11) 
1.00(-6) —-1.27(-10)6.57(-13)—-1.10(-12)—.2.93(-7)—L.01(—13)_——‘1.56(—7) 1.58(-10) 9.99(-10) 5.67(-8) —4.27(-8) —-2.13(-11) 
S.0M-7) — 1.65(-10) —-3.30(-13) —5.52(-13)4.14(-7) 5.07(-14)_—«1.56(~7) 2.21(-10) —1.39(-9) 7.91(-8)  6.08(-8) —_2.97(-11) 
2.50(-7) -2.26(-10)_:1.67(-13) _—-2.76(-13)5.84(-7) —-.2.53(-14)__—‘1.56(~7) 3.14(-10) — 1.39(—9) 1.12(-7) 6.23(-8) —4.18(—11) 
1.07(-7) — 3.40(-10)_—7.26(-14)—-1.18(-13) -8.90(-7)_—-1.08(-14) _—‘1.56(—-7) 4.88(-10)  1.40(-9) 1.72(-7) 6.23(-8) —-6.38(—11) 
7.16(-8) —4.15(-10) 4.93(-14) 7.89(-14) —1.09(-6) _——7.24(—15) —_—1.56(—7) 5.84(-10)  1.60(--9) 2.10(-7)  6.23(-8) —7.81(—11) 
4.46(-8) — 5.25(-10) 3.18(-14) 4.92(-14)1.38(-6) —4.52(-15)—‘'1.56(—7) 7.31(-10) — 2.01(-9) 2.65(-7) 6.23(-8) —-9.89(-11) 
2.23(-8) —7.00(-10):1.67(-14)-2.47(-14) —-1,.94(-6) —-2.26(-15) _—«1.56(—7) 7.99(-10)  2.65(—9) 3.49(-7) 6.23(-8) —-:1.40(-10) 
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Appendix 2B. Spatial Distributions of Absorbed Dose and 
Dose Equivalent in Cylindrical Phantom Due to 
Incident Monoenergetic Neutrons 


This appendix presents calculated absorbed doses 
and dose equivalents in a four-element cylindrical tissue 
phantom due to a broad beam of 0.025-eV to 14-MeV 
monoenergetic neutrons incident on the phantom. The 
calculations, performed by Auxier ef al! for the 
phantom schematic shown in Fig. 2.20, considered 
dose contributions by both the incident neutrons and 
the capture gamma rays they produced. The neutron 
beam was assumed to be incident normal to the axis of 


the cylinder, and the absorbed doses and dose equiva- 
lents given in Table 2B.1 each represent the average 
value over a radial volume element in layer 3, which 
is located at the midplane of the cylinder. These values 
are multicollision doses per unit fluence and thus are 
fluence-to-dose conversion factors for cylindrical phan- 
toms. 

Table 2B.2 gives the absorbed doses due to the sec- 
ondary (capture) gamma rays only. 


Table 2B.1. Absorbed Dose and Dose Equivalent Due to Monoenergetic Neutrons Incident on a Cylindrical Phantom? 


Absorbed Dose [107!° rads/(n/em?)] for E,, of? 


Layer 3 Dose Equivalent [107° rems/(n/cm?)] for E,, of? 
Volume 0.025 1 10 100 1 10 100 0.025 1 10 100 1 10 100 
Element® eV eV eV eV keV keV keV eV eV eV eV keV keV keV 
1 0.608 1.114 1.533 1.600 1.427 1.616 1.794 0.789 1.472 1.752 2.059 1.912 2.085 2.981 
2 0.474 1.006 1.324 1.378 1.036 1.378 1.783 0.618 1.179 1.550 1.725 1.422 1.799 2.559 
3 0.415 0.562 0.870 1.020 0.759 0.834 1.107 0.466 0.651 0.973 1.174 0.978 1.128 1.530 
4 0.336 0.574 0.789 0.655 0.647 0.737 1.016 0.401 0.630 0.846 0.746 0.824 0.872 1.322 
5 1.175 1.816 2.502 2.360 2.404 2.358 2.796 1.427 2.653 3.489 3.653 3.447 3.790 4.775 
6 0.801 1.361 1.587 1.700 1.660 1.743 2.385 0.962 1.695 2.130 2.451 2.405 2.677 3.785 
7 0.358 0.727 0.689 0.766 0.622 1.041 0.784 0.403 0.788 0.780 0.911 0.821 1.368 1.140 
8 0.164 0.377 0.333 0.414 0.394 0.474 0.553 0.199 0.408 0.351 0.452 0.451 0.536 0.652 
9 1.762 3.467 3.964 3.800 3.928 3.926 4.592 2.386 5.580 6.876 6.896 7.237 7.827 9.324 
10 1.161 1.952 2.178 2.364 2.336 2.314 2.961 1.583 2.892 3.560 3.943 3.858 4.103 5.911 
11 0.269 0.549 0.500 0.591 0.636 0.655 0.917 0.310 0.645 0.569 0.742 0.820 0.981 1.347 
12 0.084 0.192 0.251 0.263 0.255 0.342 0.436 0.090 0.206 0.252 0.301 0.267 0.368 0.493 
13 3.141 5.306 5.816 5.482 5.185 4.844 6.333 5.560 10.566 12.416 11.758 11.148 10.766 17.638 
14 1.603 2.876 3.014 3.022 3.025 2.651 3.280 2.640 5.380 5.816 5.849 6.058 5.441 8.515 
15 0.239 0.378 0.420 0.560 0.462 0.262 0.687 0.272 0.536 0.570 0.741 0.598 0.521 1.044 
16 0.081 0.144 0.167 0.184 0.178 0.145 0.212 0.089 0.148 0.167 0.196 0.183 0.167 0.242 
17 4.680 5.890 5.179 4.449 4.322 4.338 8.018 11.530 13.416 12.096 10.053 8.852 9.916 48.559 
18 2.379 3.130 2.989 2.761 2.548 2.538 5.018 5.500 6.871 6.462 5.586 5.191 5.559 29.011 
19 0.205 0.252 0.347 0.387 0.295 0.299 0.412 0.245 0.376 0.459 0.502 0.376 0.422 0.840 
20 0.052 0.090 0.086 0.141 0.092 0.117 0.112 0.056 0.091 0.086 0.145 0.094 0.119 0.118 


OL 


c 


Table 2B.1 (continued) 


Layer 3 
Volume 500 
Element® keV 


1 2.870 
2 2.608 . 
3 2.286 
4 1.882 
5 4.712 
6 3.372 
7 1.744 
8 1.274 
9 6.903 
10 4.332 
11 1.300 
12 0.685 
13 11.631 
14 6.627 
15 1.034 
16 0.312 
17 18.110 
18 13.619 
19 0.991 
20 0.225 


Absorbed Dose 10719 rads/(n/cm?)] for E,, of? 


1 
MeV 


4.79 
2.73 
2.11 
1.94 
5.92 
3.71 
1.80 
1.15 
9.37 
5.31 
1.51 
0.52 
17.30 
9.67 
1.39 
0.37 
30.14 
21.33 
1.55 
0.22 


2.5 
MeV 


13.1 
12.6 


5 
MeV 


39.5 
21.0 
21.9 
18.2 
37.5 
31.0 
22.2 
12.0 
42.3 
36.8 
17.9 
11.4 
48.9 
38.6 
14.7 

6.0 
57.2 
46.7 
15.2 

6.7 


7 
MeV 


52.8 
39.5 
35.5 
27.5 
39.8 
36.5 
25.2 
21.0 
49.7 
46.2 
24.8 
11.7 
54.6 
49.5 
22.7 
11.3 
57.0 
55.3 
24.7 

8.4 


Dose Equivalent [10 


-10 


rems/(n/cm7)} for E,, of b 


10 
MeV 


43.4 
43.5 
42.3 
44.2 
56.5 
45.5 
32.8 
21.7 
54.5 
51.5 
35.9 
28.5 
58.0 
56.6 
32.4 
20.6 
72.5 
64.7 
30.1 


17.4 — 


14 
MeV 


62.1 
47.4 
50.7 
42.3 
50.0 
49.6 
36.1 
38.6 
68.1 
56.7 
43.0 
27.0 
75.6 
66.1 
39.7 
22.4 
83.1 
77.4 
38.4 
18.4 


500 
keV 


6.192 
5.601 
3.627 
2.635 
14.444 
9.020 
4.401 
1.877 
33.346 
18.130 
2.561 
0.894 
87.508 
43.063 
2.814 
0.370 
188.500 
140.369 
5.766 
0.225 


i 


2.5 5 7 10 14 
MeV MeV MeV MeV MeV 
106.6 274.2 353.9 293.7 379.7 
93.2 163.3 277.3 279.3 324.6 
94.9 174.9 242.3 280.6 336.0 
96.0 139.2 229.0 279.9 322.8 
129.8 272.3 295.6 364.0 369.5 
128.9 231.5 272.1 284.3 354.9 
69.0 165.6 182.9 212.9 238.0 
42.5 102.3 151.5 149.1 245.9 
208.1 333.1 362.1 360.8 491.4 
161.3 285.4 324.8 336.7 427.6 
57.5 132.4 182.2 239.9 317.5 
19.3 85.8 89.7 183.9 185.8 
303.9 378.4 394.2 362.3 §54.2 
217.4 298.4 358.7 357.7 482.9 
58.4 110.4 162.7 209.7 286.5 
18.5 43.2 85.4 146.2 144.4 
349.6 440.7 402.9 431.3 614.9 
301.6 361.7 391.3 403.1 592.5 
70.4 117.5 169.2 193.2 278.2 
11.8 $2.3 58.9 116.1 131.0 


*From ref. 51; values include contributions from capture gamma rays. 
bE n is energy of incident neutrons. 


“See Fig. 2.20. 
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Table 2B.2. Absorbed Dose Due to Capture Gamma Rays Produced by Monoenergetic Neutrons Incident on a Cylindrical Phantom? 


Layer 3 Absorbed Dose [107'° rad/(n/cm?)] for E,, of? 
ee 0.025eV leV 10eV. 100eV ikeV 10keV 100keV SO0keV 1MeV  2.5MeV  5MeV 7MeV- 10MeV- 14MeV 
1 0.590 1.079 1.511 1.554 1.379 1.569 1.676 2.561 2.600 4.810 4.24 3.19 3.46 10.1 
2 0.460 0.989 1.302 1.344 0.998 1.336 1.702 2.307 2.133 2.851 3.55 3.53 4.59 6.08 
3 0.410 0.553 0.860 1.005 0.737 0.805 1.065 2.147 1.781 3.040 2.76 3.17 4.21 5.78 
4 0.330 0.569 (0.783 0.646 0.630 0.724 0.986 1.807 1.518 2.238 3.42 2.93 4.43 6.42 
5 1.150 1.733 2.404 2.232 2.300 2.216 2.599 3.744 3.391 3.637 4.13 3.27 4.60 8.16 
6 0.785 1.328 1.533 1.625 1.586 1.650 2.244 2.813 2.607 2.940 3.26 3.18 4.52 6.31 
7 0.354 0.721 0.680 0.752 0.602 1.009 0.749 1.495 1.457 2.899 2.76 2.55 3.71 5.81 
8 0.160 0.374 = 0.331 0.410 0.388 0.468 0.543 1.209 0.914 2.563 2.28 2.38 3.22 6.16 
9 1.700 3.257 3.675 3.493 3.599 3.538 4.093 4.426 4.396 3.478 3.34 2.80 5.02 8.36 os 
10 1.120 1.859 2.041 2.207 2.185 2.136 2.653 3.020 2.707 3.012 2.84 2.68 4.36 6.63 Ss) 
11 0.265 0.540 0.493 0.576 0.618 0.623 0.873 1.170 1.099 2.010 2.34 2.19 3.40 4.77 
12 0.083 0.191 0.251 0.259 0.254 0.339 0.430 0.663 0.443 1.430 1.74 1.91 3.14 3.72 
13 2.900 4.784 5.161 4.859 4.590 4.233 4.988 4.668 3.675 3.194 2.48 2.39 4.98 7.21 
14 1.500 2.628 2.756 2.741 2.723 2.365 2.670 3.253 2.629 2.786 2.51 2.27 4.23 6.09 
15 0.236 0.362 0.405 0.542 0.449 0.236 0.650 0.861 0.777 1.405 1.60 1.72 2.77 4.41 
16 0.080 0.144 0.167 0.183 0.178 0.143 0.209 0.306 0.334 0.762 1.08 1.31 1.71 3.14 
17 4.000 5.143 4.492 3.888 3.827 3.420 3.309 2.800 2.230 1.840 1.48 1.68 3.79 7.21 
18 . 2.070 2.759 2.644 2.478 2.260 2.039 2.221 2.269 1.985 1.744 1.65 1.79 3.52 5.69 
19 0.201 0.240 0.336 0.375 0.287 0.286 0.358 0.558 0.538 1.049 1.15 1.26 2.07 3.78 
20 0.051 0.090 0.086 0.141 0.092 0.117 0.111 0.224 0.186 0.695 0.459 0.76 1.20 1.74 


"From ref. 51; capture-gamma-ray dose largely due to 1H(n,y)?H reaction. 
bE n iS energy of incident neutrons. 
“See Fig. 2.20. 
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